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Z 

Der Spannungszustand im Ubergangsgehiet einer rechtwinkligen 
: Rohrabzweigung* 

- Von W. Reidelbach 


1. Kinleitung. In der Wand eines geraden, kreiszylindrischen Rohres, das durch Innendruck 
_ belastet ist, herrscht in einiger Entfernung von seinen Enden ein gleichférmiger, in erster Naherung 
_ ebener Spannungszustand, dessen Hauptspannungen parallel und senkrecht zu den Mantellinien der 
-Rohrwand gerichtet sind. Ihre Berechnung bereitet auch bei dickwandigen Rohren keine Schwierig- 
_keiten. Ebenfalls bekannt sind die Randstérungen in der Nahe der Rohrenden, die z. B. durch 
_ Flansche oder — bei Behaltern und Kesseln — durch 
z Deckel bewirkt werden und hauptsachlich als sekundre, 
_ relativ schnell abklingende Biegespannungen auftreten. 
_ Theoretisch unerforscht war bisher die Stérung in der 
_ Umgebung einer Rohrabzweigung (Abb. 1). Die Wand 
_ des durchgehenden Rohres ist unterbrochen und mit- 
_ hin der ebene Grundspannungszustand gestért durch 


abzwelgendes 
Rohr (Rohr) 


Teenie 


7 


vee 


; das Loch, in welches das abzweigende Rohr miindet. 
_ Das Loch wirkt ahnlich wie eine Bohrung in gespannten, 
_ ebenen Blechen: Man hat zunachst mit betrachtlich 
8 gesteigerten Hauptspannungen des Membranspannungs- 
-zustandes zu rechnen. Auf®erdem greifen infolge der durehgehendes 
_ achsialen Innendruckbelastung des abzweigenden Rohres Rohr (Rol) 
am Rande des Miindungsloches auch Krafte normal © 
zur durchbohrten Rohrwand und Biegemomente an, die AD 
. sekundare Biegespannungen bewirken. . 
__ Bekanntlich ist in einer nicht durchbrochenen Rohrwand die Umfangsspannung ungefahr (in 
einem diinnwandigen Rohr genau) doppelt so groB wie die achsiale Spannung. Daraus folgt, da 
der zu einem Umfangskreis parallele Durchmesser des Miimdungsloches starker als der achsen- 
_ paraliele Durchmesser gedehnt, das Loch also ,,ovalisiert“‘ wird. Diese Deformation wird einerseits 
durch das abzweigende Rohr behindert, diesem jedoch andererseits aufgezwungen, so daf als Rand- 
-stérung am einmiindenden Ende dieses Rohres radiale und tangentiale Krafte und — infolge der- 
genseitigen EKinspannung der Rohrwande langs des Lochrandes — Biegemomente auftreten. 

Zur theoretischen Behandlung des Problems wird daher die Biegetheorie der Kreiszylinderschale 

rangezogen, und zwar in einer vereinfachten Form, die u. a. von A. E. Green und W. Zernal? 
atwickelt und bereits auf ahnliche Probleme mit Erfolg angewendet wurde®*. Die daraus ent- 
wickelte Lésung des vorliegenden Problems stellt nur eine erste Naherung dar. Auf gréBere Genauig- 
it wurde bewuBt verzichtet, um den Umfang der numerischen Rechenarbeit méglichst klein zu 
Iten. Durch den Einsatz geeigneter, moderner Rechengerate kénnen jedoch genauere Ergebnisse 
ohne weiteres erziclt werden. 
_Bis auf wenige Ausnahmen werden im folgenden alle dem abzweigenden Rohre zugeordneten 
‘6Ben durch kleine, dem durchgehenden Rohre zugeordneten GréBen durch groBe Buchstaben 
yekennzeichnet. Der aus den beiden Rohren gebildete Kérper ist zweifach symmetrisch. Die Sym- 
metr ebene, welche dic beiden Rohrachsen enthalt, heiBe Lingsebene, die dazu senkrechte, durch 
\chse des abzweigenden Rohres gehende Symmetrieebene heiBe Querebene. 


VW 


CER oe 


» 


_ * Habilitationsarbeit Karlsruhe 1960. 
1 W. Zerna, Ing.-Arch. 20 (1952), S. 357. . 
2 4, E.Green und W. Zerna, Theoretical Elasticity, S. 385, Oxford 1954. ...= ; 
rae I. Lurje, Prikladnaja Matematika i Mekhanika, Neue Serie 10, 3, Moskau 1946, S. 397. Auszug in 
her Sprache in G. N. Sawin, Spannungserhdhung am Rande von Léchern, S. 130, Berlin 1956. 
D. Withum, Ing.-Arch. 26 (1958), S. 435. 
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2. Koordinaten, geometrische Beziehungen. Es sei vorausgesetzt, da der mittlere Durchmesser 2 a 
des abzweigenden Rohres I kleiner als oder héchstens ebenso groB wie der mittlere Durchmesser 2 A 
des durchlaufenden Rohres II ist, da®B die Rohrachsen senkrecht auf einander stehen und sich in 

einem Punkte O schneiden. Dieser 

aE Punkt O sei der Ursprung der parallel 
zu den Rohrachsen gemessenen Ko- 

ordinaten x bzw. X und = /a bzw. 
& = X/A (Abb. 2). Die Langsebene 
schneidet die Schalenmittelflachen I 
baw. II in den Mantellinien | bzw. 
L, von denen aus die Umfangs- 
koordinaten y bzw. Y und die Zentri- 
winkel y = y/a baw. O = Y/A gezahlt 


werden. In irgend einem Punkte p 


dort giiltigen Koordinatenrichtungen 
x und y baw. X und Y durch dazu 
senkrechte, also zur Mittelflache nor- 
male Bezugsachsen z bzw. Z jeweils zu 
einem dreiachsigen, cartesischen Sy- 
stem erginzt, wobei z bzw. Z nach 
auBen positiv zu rechnen sind. 

Bezeichnet y = a/A das Verhilt- 
nis der Mittelflachenradien und €,, ¢, 
bzw. &,, ®, die Koordinaten eines 
Punktes P, der Schnittkurve S der 
beiden Mittelflachen, so ist 


sin ©, = v sing, , 
5,=v cos g,, (2.1) 
1 
f= — 00s Dy, 


Die in P, an die Schnittkurve gelegte 
Tangente t schlieBt mit den durch P; 
gehenden Bezugsachsen y baw. Y die 
Winkel y bzw. Y ein. Aus Abb. 3 
folgt, wenn man eine lings t abge- 
messene Langeneinheit auf die y- bzw. 
Y-Achse und auf die gemeinsame Nor- 
male n der x- bzw. X-Achse projiziert, 


cos y cos y, = cos ¥ cos ®, , 
sin y = cos VY sin @, , 
sin Y = cosy sing, , 
und wenn cos y/cos ®, = cos Y/cos @, 
= C gesetzt wird, 
sin y = C cos 9, sin OD, | (2.2) 
cosy = C cos @, , 


sin ¥W = C cos Q, sin g, , 7 
Abb. 3. (2.3) 
cos ¥= C cos q, . 
Die Tangente t ist die Schnittgerade der Tangentialebenen y bzw. I", welche die Mittelflachen I 
bzw. IT im Punkte P, beriihren. Die Spuren k bzw. K der Tangentialebenen in einer auf t senk- 
rechten Ebene ¢ durch P, schlieBen den Winkel y ein (Abb. 4), ebenso die durch P, gehenden 


bzw. P der Mittelflachen werden die _ | 
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Bezugsachsen 2 baw. Z. Wird von P, aus lings Z eine Langeneinheit abgemessen und auf k projiziert, 
so hat die Projektion die Lange cos @,/cos y und das projizierende Lot die Lange sing, sin ® 

Daher ist s “ 
cos®, 1 


sree cosy eC.” 


cos 7 = sin y, sin O, = y sin? g, , (2.4) 


2 1 1 

eter eae nee ) 
Innerhalb des Koordinatensystems (é, gy) auf der 
Mittelflache I hat die Schnittkurve die Gleichung 


0 ¢ 05 1 


SSS 
e => i" /1— sin? d, = Wer sin2 Ys ; (2.6) 


im Koordinatensystem (2, @) auf der Mittel- 
flache II die Gleichung 


Rae y1 ——\sin? @, = yp? — sin? @, 52 (2.7) 


Abb. 4, Abb. 5. 


Denkt man sich die Mittelflache II auf die langs der Mantellinie L beriihrende Tangentialebene 
abgewickelt, dann nehmen die Schnittkurven die in Abb. 5 dargestellten Formen an. Fir vy = 1 
(a = A) ergibt sich eine cos-Linie, fiir y < 1 ellipsenahnliche Kurven, die um so weniger von einem 
Kreis abweichen, je kleiner y ist. Der Schnittpunkt ®; einer solchen Kurve mit der Y-Achse folgt 
aus.2, — 0, dh, y= sin O*, zu 


i 3 5 
* : 2 4 6. cla aeee 
@; are sin y »(1 |- aa iS age + 336” - ): 


Die gréBte radiale Abweichung der Schnittkurve mit dem Parameter y vom Kreis mit dem Halb- 
messer v, bezogen auf y, 

aa D* —y 1 : See 

= Sent y 10 y - tele 


Ay 


v 


bleibt kleiner als 5° (Ay < 0,05) fiir y < 0,5. Man kann daher im Bereich 0 <y <0,5 gegebenen- 
falls statt der Schnittkurve (2.7) einen Kreis 


2 = y— PC oder , = cos g, ; @, = ysing, (2.8) 


verwenden. Aus (2.2) folgt mit (2.3) und (2.5) 


5 sin? Y cos? ®, 1 — 1? sin? y, Fe) 
cos = 5 — > = . 4 - 
¥ sin? y, sin? ¥ 1 — 7? sin‘ 9, 


: 1 1 Ui 
sin” p, = ye = anye 


Dieser Quotient erreicht mit 


21* 
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den kleinsten Wert 
sin? Y| y (2.9) 


(cos? Y) min Te Fe Gai 9 (1 Seth V1 —**) f 


Im Bereich 0 < »y < 0,5 ist demnach 


1= cay aoe 0,97 , 
also 
We > cosprH 1. (2.10) 


Aus (2.2) folgt mit (2.1) 
siny = Ccosgy,y sing, = Cy sin 2 9, » 
so daB im Bereich 0 <» < 0,5 auch die Naherungen 
sin2y = 2sinycosy » 2siny = Cysin2@Q, , cos2ypr 1 (2.11) 
zur Verfiigung stehen. SchlieBlich seien noch zwei spater auftretende Faktoren betrachtet: 


yi+ cos Y = yl + y sin? y, © 1+ sin*g, = 144 (1 — cos 29,) . (2.12) 


Im Hinblick auf die in Abschnitt 6 formulierten Randbedingungen, die langs der Schnittkurve, 
bei hinreichend kleinem y also langs eines Kreises erfiillt sein miissen, ist es zweckmaBig, auf der 
Mittelflache II ein polares Koordinatensystem (r, pm) bzw. (0 = r/a, y) zu 
definieren, dessen Ursprung im DurchstoBpunkt der Rohrachse I durch 
die Mittelflache IT liegen soll (Abb. 6). Es gelten dann die Beziehungen 


= 7 60s @, Y= sus (2.13) 


£=y 0 cosy, P=—voesng. 


Fiir @ = 1 (konstante Koordinate aller Punkte einer nahezu kreisf6r- 
Abb. 6. migen Schnittkurve) geht (2.13) in (2.8) iiber. 


3. Verschiebungen, Schnittgréfen. Die Verschiebung irgend eines Punktes p bzw. P der Schalen- 
mittelflachen resultiert im allgemeinen aus je drei Komponenten u, v, w bzw. U, V, W. In Abb. 2 
sind die Richtungen positiver Verschiebungskomponenten (im folgenden kurz Verschiebungen 


ie 
MN, 
Ngy I= S ee : 
Nop 
= eae Se 
—_ Nyg 
N, 
rs ae 
: m, 
(Ages ; 
PSG we 7 My 
ly = IL ee Dee be Myr 
aa Pa) Mp eee 
Mya ™y M, Mp 
ee Up 
My Ie Qo Ne 
Moy Op 
Abb. 7 
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genannt) angegeben. Es seien partielle Ableitungen der Verschiebungen nach den Koordinaten Y 
bzw. ® durch Punkte, Ableitungen nach & bzw. © durch Striche gekennzeichnet: 


ou . ou ; eU oU 


os 9E =U > 36. = Us? eee usw. 
Die positiven Schnittgré8en zeigen Abb. 7 und 8. Die Vektoren der Schnittmomente seien als 
Rechtsschrauben definiert. Die Richtungen der positiven Momentenvektoren wurden so gewahlt, 
das in beiden Koordinatensystemen (x, y, z) baw. (X, Y, Z) formal identische Grundgleichungen 
gelten. Man erhalt dann aus den in Abschnitt 5 dargestellten Grundgleichungen fiir das Rohr I 
die entsprechenden Gleichungen fiir das Rohr II, indem man die Kleinen Buchstaben durch die 
entsprechenden grofen Buchstaben ersetzt. 
In bezug auf das durch (2.13) definierte polare Koordinatensystem (0, y) treten statt U, V, W 
die Verschiebungen 
U*= Ucosp + Vsing, 
vV* = — Using + V cos , 
ae | 


auf und statt der SchnittgréBen nach Abb. 8 die SchnittgréBen nach Abb. 9. 


(13.1) 


4, Der Grundspannungszustand. In hinreichend groBer Entfernung von der Schnittkurve erzeugt 
der Innendruck p rotationssymmetrische Spannungs- und Formanderungszustande mit folgenden, 
durch ~ gekennzeichneten Gréfen (FE Elastizitatsmodul, bt Querzahl, h baw. H Dicke der Rohr- 
wande, &) eine Konstante): 


Rohr I Rohr II 
a i ~ A2 nz 
[v= Pe, 02m) E—4), U =F (l—2n) &, 
Verschiebungen: 
| a) ° a0 
op a Lt De (ye 
2). F-44(-9) 
: ~ pa ee pa 
Achsialspannung : n= >> x= 
Umfangsspannung: 7,=pa, peal; 


oder fiir Rohr II mit Bezug auf das Koordinatensystem (o, ¢): 


~ se ~ pA 
pees 2H yo (1+ cos29), N, ="= (3 — cos 29), 


~ EH 4, 
ny eT 9 ~ pA 
y*——7 A Pvosin2g, Ca ee as 
iat p A LM a ipa : 
ameal 4): Nes i sin 2» 


Die tatsachlichen Spannungs- und Formanderungszustande ergeben sich dadurch, daf man den 
Grundzustanden eine Randstérung iiberlagert, die durch eine in den Punkten der Schnittkurve 
angreifende Gleichgewichtsgruppe von (inneren) Schnittkraften und -momenten hervorgerufen wird. 
Der zu iiberlagernde Stérungszustand wird im folgenden mit Hilfe der von Green und Zerna ent- 
wickelten Naherungstheorie fiir die biegesteife Kreiszylinderschale bestimmt. 


5. Grundgleichungen der Biegetheorie der Kreiszylinderschale. Aus den Verschiebungen und 
ihren Ableitungen berechnet man zunachst die Verzerrungen der Schalenmittelflache 


e,=4, eg=— t=, Vey = + (5.1) 
und daraus die in der Mittelflache wirkenden Schnittkrafte 
n, = E*h(¢,-+ ey) 
Bee heme), | zen SE 62) 
ie be 


Nyy = My, = E*h 5) Vxy? 
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Die Schnittmomente ergeben sich aus den Ableitungen der Verschiebung w zu 


E*h Ae : 
Hiempsal T Ue | 
Ti hs ee ” 
my = ae wt’), (53) 
eo ; 
mM, y = My» = 12 a (1— yp) w | 


Da bei Einleitung einer reinen Randstérung die innere oder 4uBere Schalenoberflache nicht belastet 
ist, unterliegen die an den vier Kanten eines Schalenelementes angreifenden SchnittgréBen folgen- 
den statischen Gleichgewichtsbedingungen: 


n, + ny =90, (5.4) 
ny + Ny — qy = 9, (5.5) 
qy + a + ny =0, (5.6) 
m, +.m,y—qy,a=0, (5.7) 
m, + mM,y—q,4= 0, (5.8) 
wahrend die Verzerrungen die Vertraglichkeitsbedingung 
e+ & —yzy—— = 0 (5.9) 


erfiillen miissen. 
Im Rahmen der Naherungstheorie kann man die Querkraft q, in der Gleichgewichtsbedingung 
(5.5) unterdriicken, so daB 


ny + ny = 0 (5.10) 
sein muB. Werden die in der Mittelflache wirkenden Schnittkrafte gemaB 
pa pe =f Ney = —f" (5.11) 


aus einer Spannungsfunktion f(&,q) abgeleitet, so sind die Gleichgewichtsbedingungen (5.4) und 
(5.10) identisch erfillt. Die Querkrafte kénnen mittels (5.7) und (5.8) unmittelbar durch die Mo- 
mente und mit (5.3) schlieBlich durch die Verschiebung w dargestellt werden (Adw= w”’ + w’): 


ie Tok (Aw)’ > (5.12) 
E* h a 

t= Tre (0 9 

Damit und mit n, = f” erhalt (5.6) die Form 

12a° 1, 
AAw +5 sf =0- (5.14) 
Aus (5.2) ergeben sich die Verzerrungen 
__ Mx — Ny __ Ny — Ux _ 20 +4) 
S4 eB Sli ee ie Vay ie She Ray 


als Funktionen der Schnittkrafte n,, n,, n,, und damit der Spannungsfunktion f, so daB (5.9) 
schlieBlich iibergeht in 


on ese 
AAf— ~~ w =a) (5.15) 
Durch Kinfithrung der komplexen Funktion 
g=wt+ikf 
werden die beiden simultanen Differentialgleichungen (5.14) und (5.15) zu einer komplexen Dif- 
ferentialgleichung 
AAg—2i% 9” = 0 (5.16) 
vereinigt. Die Koeffizienten 


a? y12 (1 — p22 Nes ee 
i = 4 3—p) (5.17) 
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haben die Bedeutung von Schalenkennzahlen. Fir Rohr II gilt analog 


C= Pek RF, AACZ 27620" = 0. (5.18) 


A2 V2 (1 — 2) Bn ON 
K=#) a ats O = | 3 (l—p). 


Die Differentialgleichung (5.16) wird durch die Exponentialfunktion 


Pes tPg ? 


(eieg ert (5.19) 
befriedigt, falls deren Eigenwerte der charakteristischen Gleichung 
Pi + 2 Pi Pp + Pp —2i 7 pi = 0 
oder 
PE+ Pot (1+ i) Opp =0 
geniigen. Mit 
z 2 = 2 
Peg Pes Pp % Po 
erhalt man schlieBlich 
Pe= + (14+ i) + y2i— pi. (5.20) 


Bezeichnet man Ableitungen nach @ durch Striche, Ableitungen nach y durch Punkte und 
werden die Verschiebungen U*, V* voriibergehend ohne Stern geschrieben, so gelten in bezug 
auf das Koordinatensystem (9, y) folgende Gleichungen: 

ae,=— U'+7Wsin?g, 
Q4e=V'+U+e0rW cog, (5.21) 
04y,,7=U+eV’—V+orWsin2o, 
N; = E*H (e, ae Le Eq) 2 
Ny = E*H (ég +  &,) 5 


(5.22) 
No eopeyp ee! 
rp Q rp? 
aa E* HP 7D Le 1 L a5 
M, ==" (v +4 ta’), 
E* A 1 1 & 
ae A Set: . == 5.23 
M, = ae (HM + +), ee 
E* HP les 1 j 
M, = = age GML W" aw), 
E* HB ? 
0, = (AW) ? 
ae AW=w’tiwsin. (5.24) 
E*He Q @ 
Q» 12 a (Avy, 


Aus den Beziehungen (2.13) zwischen den Koordinaten erhalt man die folgenden allgemeinen 
Beziehungen zwischen den Differentialquotienten: 


e Ce eee Uso ee ee eeeae| ne 
0B? (vo? cos lee 500) ' Gor a seams e. ap vg epee); 
@ pate coe Lever te ina ae ano. 
6% — Boge YP + & Hoe) * Gor a es ie epee) =e) ep i 


die es gestatten, die Differentialgleichung (5.18) zu transformieren in 


e gee Hee ef teal eG. TS “ 
ie + 56 300) * GOP ey laea v0 Bo) * Woy ay 


ae 1 eG]. ery 
2404 rae oo + [56 apa) t GeF aet| ome + lay op + oe Bpava| 27 
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Multipliziert man noch mit 74 und ersetzt man i durch — i (Betrachtung der es komplexen 

1 ; : 
Gleichung), 162? durch 7? und cos? durch a (1+ cos2q), sin?g durch a (1 — cos 2q@), so 
entsteht die im (0, y)-System giiltige Differentialgleichung 


Lane 
AAG + 17 AG +43 (6" — 76 — 76") cos 29 29(3 CG ") sin2y =0. (5.25) 


o 
Hierin ist 
G=W—ikrF, (5.26) 
und die Schnittkrafte in der Schalenmittelflache sind jetzt nach den Formeln 
Noe . F’ + =P, N, 2 PS Ne a Fe F’ (5.27) 


aus der Spannungsfunktion F zu berechnen. Der Lésung (5.19) entspricht hier die Funktion 


Eg Chee eT Le Come (5.28) 


6. Symmetrie- und Randbedingungen. Da die Formanderungen beider Rohre sowohl zur Langs- 
ebene als auch zur Querebene symmetrisch sind, miissen die Verschiebungen u, v, w symmetrisch 
in bezug auf y = 0, 2/2, x, 3/2 und die Verschicbungen U, V, W symmetrisch in bezug auf 
® = 0, a und © = 0 sein. Ferner ist zu beachten, daB in den Symmetrieebenen keine Schubkrafte 
und Drillmomente iibertragen werden. Daraus folgen als Symmetriebedingungen 


ve nN 32 v=0, ee a (6.1) 
fir p= 0, 55% >? 

Ney = 9, gqy=0, my=0, (6.2) 

= ets Wie ia) en rs : 

fur-O = 05-7: és (3) 

Nyy =90, Oy=—90, Myy=0, (6.4) 

far 8 0: OU: '6,.7. V2 032 WO 0 (6.5) 

Nyy = 04 Qy = 0, Myy = . (6.6) 


Die Randbedingungen bringen die folgenden geometrischen oder mechanischen Gegebenheiten 
zum Ausdruck: 


a) Die Randstérungen, das sind die langs der Schnittkurve S auftretenden Abweichungen der 
Formanderungen und Spannungen von den in Abschnitt 4 angegebenen ZustandsgréBen des Grund- 
spannungszustandes, klingen mit zunehmender Entfernung von der Schnittkurve relativ schnell 
ab und verschwinden nach dem Saint-Venantschen Prinzip in einem endlichen Abstand praktisch 
volistandig. Die Funktionen g bzw. G miissen demnach abklingende Funktionen sein. 


Abb, 10. 


b) Der Verschicbungsvektor eines Schnittkurvenpunktes P, mu8 aus den Komponenten 


pee a epee resultieren wie aus den Komponenten U,, V,, W,. Diese Bedingung verlangt 


w, sing, + v, cosy, = W, sin D, + V, cos ®, , | 
u, = W, cos 6, — V, sin 9, , (6.7) 


U, = COE Ys any ee 


I 
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Verwendet man statt u,, v, bzw. U,, V, die zur Tangente t senkrechten und parallelen Komponenten 


u, = u, cosy + v, sin, 


v, = — u, sin y + 2, Cos yp , os) 
bzw. 

U, = U,cos ¥ + V, sin Y, 

V,=— U, sin Y + V, cos ¥ C2) 


und betrachtet man (Abb. 11) die in der Ebene ¢ liegenden Komponenten u,, w, bzw. U,, W,, so 
lauten die entsprechenden Gleichungen 


DLE as 
U, cos 4 —w, sin 4 = — U, cos $+ W, sin fe, 
U, sin 4 + 1, cos 4 — U, sin 4+ W, cos 4, 
A eh dame a bas pS Sr, te LTA ae 

oder wegen sin — =| 5 (} — cos x), cos Cae |) (1 + cos y) 
0=2,—V,, (6.10) 
0 = (u,+ U,) /1 + cos 7 — (w, + W,) y1—cosy, (6.11) 
0 = (u,— U,) (1 — cos x + (w, — W,) V1 + cos yx. (6.12) 


_ Die beiden Tangentialebenen, welche die Mittelflachen I bzw. II im Punkte P, berihren, 
drehen sich um die Tangente ft (d.i. ihre Schnittgerade) um gleiche Winkel, d. h. so, daB der von 
ihnen eingeschlossene Winkel konstant bleibt. Beschreibt man die Drehung der Tangentialebenen 
- durch rechtsschraubende Drehvektoren ,, w,, w, baw. Qy, Qy, Qz auf den Bezugsachsen x, y, z 
baw. X, Y, Z (Abb. 12), so kommt 


w W: 
oie Oa a? 
w’ Ww 
Det riee vee eae 
ae yee 


Der eingeschlossene Winkel bleibt konstant, wenn die auf t fallenden Komponenten 
WO, = Wy Cosy — O, Sin , 
2, = Qy cos ¥ — Dy sin ¥ 

der resultierenden Drehvektoren gleich sind (0 = w,— 2,): 


0 =~ cosp 4+ ~ siny + 4 cos W + sin. (6.13) 

c) Die Schnittkrafte und -momente, die im Ubergangsquerschnitt vom Rohr I auf das Rohr II 

bzw. langs der Schnittkurve von der Mittelflache I auf die Mittelfliche II und umgekehrt aus- 

geubt werden, sind nach dem Wechselwirkungsprinzip entgegengesetzt gleich. Im Punkte P, der 
~ Schnittkurve wirken im Rohr I (Abb. 13) die Schnittkrafte 


1 
n,= 5 (n, + ny) + > (m,— Ny) Cos 2y + Nzy sin 2p 


6.14 
, =— + (m% —1,) sin 2p + my cos 2p, ( ) 


q, = 4.008 y + gy siny 


302 W. Reidelbach: Spannungszustand im Ubergangsgebiet einer Rohrabzweigung Ingenieur-Archiv 


und die Schnittmomente 


m, ae (m,, + my) +e (m, — my) cos 2y + m,, sin2y, 
: (6.15) 
M59 — = (m, — my) sin 2 p =e M, y COS 2 , 


Xe i 
Ng 


Pip se 
| ge r / Mx 


Th 


Ys 


te iby 
ke SE 
Ty 
Msy Mx \Q% 


Abb. 12. Abb, 13. Abb. 14. 


im Rohr II (Abb. 14) die Schnittkrafte 


N, = : (Nx + Ny) + 2 (Nx — Ny) cos2 4+ Nyy sin2 ¥, 
T, = —— (Nx—Ny) sin2 W+ Nyy cos2¥, Ce 


Q, = Qx cos ¥ + Qy sin Y 


und die Schnittmomente 


M, =-> (Mx + My) + + (My — My) cos2¥ + Myy sin2¥, (6.11) 
M,y = — > (Mz — My) sin2 ¥ + Mgy cos 2. 


Die Krafte t, und T, sowie die Momente m, und M, liegen auf der Tangente #, die Krafte n,, q,, N. 


s 


und Q, sowie die Momente m,, und M,y samtlich in der zu t senkrechten Ebene ¢ (Abb. 15). Die 


Querkrafte q, und Q, schlieBen den Winkel y ein. Das Wechselwirkungsprinzip stellt die folgenden 
Bedingungen: 


T,=—t,, M, = m, , 
Q, cos 4 — N, sin % = —q, cos 4 + n, sin 2, 
Q, sin $+ NN, cos 4 = qs sin 4 + n, cos 4 
oder | 
0= M,—m,, (6.18) | 
O= 7, + t,, (6.19) | 
0 = (Q, + 4) V1 + cos y—(N, + n,) /1— cos x, (6.20) 


0 = (Q,—g,) f1— cosy + (NY, —n,) VIF 08 7. , (6.21) _ 
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Die Randbedingung fir die Drillmomente m,, und M,y ergibt sich aus der Betrachtung eines 
Schalenelementes am Ubergang zwischen den Rohren (Abb. 16). Die an den Schnittflachen des 
Elementes angreifenden Schubspannungen sind — nach den iiblichen Voraussetzungen der Schalen- 
Aiegetheorie — lineare Funktionen der Koordinaten z oder £ = 2 z /him Rohr I und Z oder = 2 Z/H 
im Rohr II. Aus den Schubspannungen resultieren die Schubkrafte 


+h/2 +1 +H/2 re 2 

w= [ ude= > [ad, T= { maZ—=F | mat (6.22) 
—h/2 ss] — H/2 —1 
und die Drillmomente 
+h/2 +1 +H/2 +1 

h? ae 

My = [meea FZ [mca Mwy= | mZdZ= 7 fant dt. . 
—hj2 = — H/2 —1 


Abb. 16. 


Jede unendlich diinne ,,Schicht“, die man durch Schnitte langs € = const und € + df = const 
aus dem Element heraustrennt (Abb. 17), ist im Gleichgewicht, wenn die Schubspannungen der 
Bedingung ‘ 3 

t1 > @=—tm = dC = dS 


geniigen. Damit folgt aus (6.22) und (6.23) zunachst 


t=+1 t=41 
j poe ig joe | ds 
==) t=—1 
c=+ t=+1 
AH 
moma | od, Mee | cds 


und schlieBlich die schon oben angegebene Randbedingung (6.19) und die weitere Randbedingung 
0=m,,H+ M,yh. (6.24) 


Ersetzt man hierin die Schubkrafte und Drillmomente mit Hilfe der Gleichungen (5.3), (5.11) und 
(6.14) bis (6.17) durch Ableitungen von w und f bzw. W und F, so kommt 


1 oe vt * Te 
q 0=5 (f° —f") sin 2y +f" cos 2p + 5 (F — F”)sin2 ¥ + F’ cos2¥, 


‘a | 2 


2 c te 
o h2 Lz (o" —#") sin 2 y -+ w” cos 2y| +a [a (W" — W”) sin 2 ¥ + W" cos 2 ¥. 
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Nach Multiplikation der ersten Gleichung mit 


; He, adhd eLae 
7 k = 0 K a he 
und der zweiten mit a2/h? erhalt man durch Addition und mit w+ikf=g, W+ tiKF=G 


1 ee 47 OQ fe H? & 1 fg vt ° Je 
0= = (g°—s")sn2yp +8 cos 2y + 75 FA ) sin 2 Y¥+ G6 cos2 ¥|. 


Die beiden Randbedingungen (6.19) und (6.24) erweisen sich hiermit als Bestandteile einer einzigen |} 
komplexen Bedingungsgleichung und sind also zwangslaufig beide erfiillt, wenn eine durch geeignete | 
Wahl der in g baw. G enthaltenen Integrationskonstanten ¢c bzw. C befriedigt wird. Damit ist das_ | 
vorliegende Randstérungsproblem zuriickgefiihrt auf acht notwendige und zugleich hinreichende _ 
Randbedingungen, die in den Gleichungen (6.10) bis (6.13) und (6.18) bis (6.21) ausgedriickt sind. | 


Grenzfalle. a) Stimmen die mittleren Rohrhalbmesser iiberein (a = A, y= 1, Abb. 18), | 
so ist 


®,= 9,» y=, 


sin y = C cos g, sing, , cos y = Ccosq,, 


cos 7 = sin’ g, , V1— cos 7 = cos @,, /1 + cosy = y1 + sin®g, . 
Die Randbedingungen lauten dann 
0 = (W, — w,) sing, + (V, — 2,) cos 9, » 
0 = W, cos y, — V, sin g, — u, , 
0 = w, cos y, — v, sing, — U, 
oder 
V=%,—V,, 
0 = (u, + U,) V1 + sin? y, — (w, + W,) cos 9, , 
0 = (u, — U,) cos p, + (w,— W,) /1 + sin? g, , 
0 = (w’ + W’) + (w' + W’) sing, , 
0= M, > Mt,» 
Oe ihe + te: 
0=(Q,+ q,) /1 + sin? y, — (N, + n,) cos g, , 
0 = (Q, —q,) cos g, + (N,—n,) /1 + sin? g, . 


Abb, 18. Abb. 19 


: b) Die Rohrhalbmesser seien gleich groB (v= 1) und RohrI zweige nach beiden Seiten von Rohr ite 
jeweils lotrecht zu dessen Achse, ab (Rohrkreuzung, Abb. 19). Die Schnitikurve zerfallt jetzt in 


zwei Ellipsen mit den Halbachsen a und a y2 in den beiden Ebenen, welche die rechten Winkel 
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zwischen den Rohrachsen jeweils halbieren und ebenfalls Symmetrieebenen sind. Daraus folgt, 


daB u = U,v = V, w= W usw. ist. Von den Bedingungen des Grenzfalles a verbleiben infolge- 
dessen nur 


0 = w, cos py, — v, sing, — u, 


oder 
Ou, V1 + sin? y, — W, COS Y, , (6.25) 
ferner 
0=w’ + w'sing,, (6.26) 
ak (6.27) 
= q, V1 + sin? y, — n, COS Q, . (6.28) 


c) Der Halbmesser a des Robhres I sei nicht gréBer als die Halfte des Halbmessers A des Rohres II 
(0<y< 0,5). Dann folgt aus (6.11) und (6.12) durch Addition bzw. Subtraktion und mit Hilfe 
der Naherung (2.12) 


0 = u, + w, = (1—cos 2 ,) + U, = (1—cos29,)— W,, 


0 =u, — (L—cos29,) —w, + U, + W,  (l—cos 29). 


Wird jetzt aus der ersten Gleichung U,, aus der zweiten u, eliminiert und werden dabei Glieder mit 
dem Faktor 72/16 vernachlassigt, so kommt 


0 = u,+ w, > (1—cos 2 9) — W,, 
0 = w,— U,— W, = (l—cos 2Q,) . 
Durch entsprechende Umformung von (6.20) und (6.21) erhalt man 
v 
0=q,+ n, > (Ll — cos 29,) —N 


0=n,—Q,—N E (1 — cos 2 9,) . 


s s ) 
_ Aufgrund der Naherung (2.10) ist im Bereich 0 Sy < 0,5 
Y. ae OW ew _ W’ 
Vi, Co Ue 5x ¢ el geemell sin Y = re a 


Mia... T= Nye» On Or. N, = N,.. 


Benutzt man jetzt noch die Transformationsformeln (6.7), (6.14), (6.15) und die Naherungen 
(2.10) und (2.11), so lauten die Randbedingungen, wenn der Index s als Hinweis auf die Schnitt- 
_kurve nicht mehr geschrieben wird, 


=—uC+sin2g+v—V*, (6.29) 
=utvC+sin2p+w-> (lL—cos29)—W, (6.30) 
0=w—Ut— WwW (1—cos29), ~ (6.31) 
0=w' +wC% sin2g+ W was, (6.32) 
0 =m, + m,, Cy sin2y9—M,, (6.33) 
0=—(n, — my) Cz sin 2p + Mey + Neg» (6.34) 
0= q+ dy 7 sin 2g +m, x (1 —cos2%)—N,, (6.35) 
0 =n, + n,, Crvsin29—Q,—N, > (L—cos 29). (6.36) 


+e imi lial 
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7. Die Lésung des Randwertproblems. Mit Riicksicht auf die Symmetriebedingungen (6.1) 
und (6.2) wahlt man zweckmaBigerweise pg = - un (n= 0,2, 4;...), so daB 5 


z, = c, oP! (en + eine) = (a, + i b,) ef cosng. (7.1) 
Wegen des Abklingens der Randstérungen kommen von den Eigenwerten p; nur diejenigen mit 
negativem Realteil in Betracht: 


per = —A,— 104, Ppe2 = —A, +1. 
Aus (5.20) erhalt man 
oO oO 
1, 2D: ee 
o 3 o = 
= 5 (We set Os A es (7.2) 
‘Ant 2n? 
ws=|/ ee es pes 
Zerlegt man 
Sn = Wy, +ikf, = (a,, + i b,1) (cos wo, & —i sin w, §) e~** cosn@ 
++ (dng + ib, 9) (cos @y& + i sin w, €) e—** cos np 
in Real- und Imaginarteil, so kommt 
Wy, = (A, COS W,& + 8, sin w, &) e~** cosnyp 3 
+ (Gy2 COS Wy & — bay sin Wy &) e—*§ cosng = w, cosn®@ , (7.3) 
k f,, = (6,1 cos w, € — a, sin w, €) e~** cosng 
+ (bys COS &-+ dyy Sin w,€) e~*£ cosng=kf,cosng. (7.4) 


Aus w,, und f,, berechnet man mit Hilfe von (5.1) bis (5.3) und (5.11) bis (5.13) alle ibrigen Ver- 
schiebungen und SchnittgréBen, die samtlich in der Form 


5 [(44 G14 — By yx) cos @, E + (By ayy + Oy 5,1) sin @, »] e~** 
+ [(%g Gre + By Ong) COS gE + (By Gng — Og bag) sin Wy é] e—~** 


erscheinen. Die Koeffizienten 6, 0, 61, %2, $2 entnimmt man Tabelle 1. Alle ZustandsgréBen stellen 
sich dar als Produkte aus sin n y oder cos ny und einem nur von  abhangigen Faktor, der wie in 
(7.3) und (7.4) durch Querstrich gekennzeichnet sei. 

Eine der Lésung (7.1) entsprechende Funktion fiir Rohr II 


Gy = (Ay +i By) 2" cos N® 


mit 
Pz, = —A,—iQ,, Psg=—A, +i, 


r=) 


wirde die Symmetriebedingungen (6.3) und (6.4) automatisch erfiillen. 


Im Hinblick auf die Bedingungen (6.5) und (6.6) beachte man zundchst, daB auch W” und 
Myy verschwinden, wenn W’ = 0 ist fiir alle ®, und daB die Bedingungen U= 0 und N xy = 0 
zugleich V’ = 0 erfassen. Es sind daher statt scheinbar sechs tatsachlich nur vier Bedingungen 
vorgegeben. Sie kénnen jedoch nur mit der trivialen Lésung Ay, = Ay. = By, = By, = 0 
erfiillt werden. Auch der Verzicht auf die Erfiillung unwesentlicher Bedingungen — etwa U = 0 — 
oder Nyy = 0 oder beide — bringt nur geringen Nutzen, denn obwohl jetzt die Konstanten nicht 
samtlich verschwinden und eine oder zwei der Konstanten frei bleiben, ist eine geschlossene An- 
passung der Lésung an die Randbedingungen aus folgenden Griinden nicht méglich: 

Alle GréBen, die den Randbedingungen zu unterwerfen sind, miissen als Funktionen einer 
laufenden Koordinate der Schnittkurve — z. B. y, — dargestellt und dann die Konstanten G13 
Ay, usw. so bestimmt werden, daB die Randbedingungen im Bereich 0 < Y, Sa/2 fir jeden Wert 
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Tabelle 1 
6 Oy By OX B 
2 
ape cosy [Ue @, — Ws — ps, —A, — tM, +o, by om 
saa se 
ge ee HA, — Wo) (A, —o,)—(1+-) n® — B(A4, + @,) —8 (1,4 @,)-+(1+ 2) n? 
w cos np 1 0 1 0 
w cosn p —A, — Oo, —A, —, 
n? : : 
Me | > cosny 0 1 0 1 
1 
ny | = cosng O (A, + @,) — 3 (A, — w,) — n? — B(A, — ws) —@0(A, + w,) + n? 
Ln 
Nz y |, snny = 0; Ay Ws —A; 
1 : 
de | app Cry 28o,—n* —20/1,—n? —2%o, +n, 2 BA, — n? 
n b 
one te — (Ay — ay) — (A, + @,) Ig ++ Wy — (Az — @) 
a 
Mm, Loz paces MPAA —o)+1—p) nF OA, + a) 8 (Ag+ 2) + (1—p) n? 8 (Ag — ws) 
a 
My | az gn 08 2 PH HA—o,)—(1—w) nw 8 (Ay + 1) M O(Ag+ @2)—(1—) n? 1 B (Az — we) 
na. 
M,y| 9], g2 Sin np Ay QO, A, We 


des Winkels ¢, erfiillt sind. Man hat dabei die Kreisfunktionen der Winkel y, Y, y nach (2.2) bis 
_ (2.5) und ferner die Argumente &,, &,, ®, der Kreisfunktionen in den Lésungen g bzw. G mittels 


ae 


(2.1), (2.6) und (2.7) durch Funktionen des Winkels y, zu ersetzen. 


Es entstehen dabei Ausdriicke von der Form 


De oe 1/2 ; 
ae (= ie *) cos | w Eases sin? mw) Cv sin? 9) (1 — sin? co 
1} y Ps y2 Ys SND, 


- (1 — »? sin? p,)1/? 


oder 
A,» cos 
e 1 P 


s cos (Q, v cos y,) cos [N arc sin (y sin g,)| C? [cos? p, — sin? y, (1 —v? sin? ¢,)] (1+ y sin? g,)"””, 


d.h. die Bestimmungsgleichungen fiir die Konstanten enthalten von g, abhangige Koeffizienten, 
die sich nicht herauskiirzen lassen. Um diese Schwierigkeit zu iiberwinden, kann man die Rand- 


_ bedingungen in einzelnen Punkten der Schnittkurve befriedigen, indem man — falls j die Zahl der 
; liber den Bereich 0 < —~, < 7/2 verteilten Punkte ist — so viele Losungen g, und Gy superponiert, 


a 


Pa ie ok aah 


= 


daB 87 freie Konstanten zur Verfiigung stehen. Da wenigstens in drei Punkten (z. B. py, = 0, 
m/4, 2/2) die Randbedingungen erfiillt werden sollten, sind mindestens 24 Konstanten zu be- 
stimmen, wobei allein die Berechnung der 24? Koeffizienten in den 24 Bestimmungsgleichungen 
einen erheblichen Aufwand an numerischer Rechenarbeit darstellt. 


Beschrankt man sich jedoch auf kleine Radienverhaltnisse y, so fallt die Schnittkurve innerhalb 
der Mittelflache I ausreichend genau auf einen Umfangskreis mit der Koordinate x, = A, inner- 
halb der Mittelflache II sehr genau auf einen Kreis mit dem Radius a. Die Randbedingungen sind 
"jetzt einerseits fiir £, = 1/y und andererseits fiir @ = 1 zu erfallen. Hierzu eignet sich fir Rohr I 


die Reihe 
g= do, eC) cosnp i els eet. (7.5) 
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die aus (7.1) durch eine Verschiebung des Nullpunkts der Koordinate € hervorgeht, und fiir Rohr Ii 
in Anlehnung an (7.5) eine Lésung in der Form 


G = S1T,(0) cos ny n=O 2, Ameo (7.6) 


in der J7,(0) eine Funktion nur von g ist und mit wachsendem @ ahnlich einer Exponentialfunktion 
mit komplexem Exponenten abklingen soll. Die strenge Lésung der Differentialgleichung (5.25), 
d.i. die Funktion (5.28), 14Bt sich nicht in dieser Form darstellen, da sie in bezug auf y = + 7/2 
nicht symmetrisch ist. 

Setzt man (7.6) in (5.25) ein, so erhilt man als Differentialgleichung fir if, 


Q 


n 


Batt, + 7 A* IT, + Tal Sih + % Or cos 2 cos a 


2 7? (3 Tf. pHs) sin 2 sin n g} = 0 (7.7) 
AERC rd n? | 
a" de Kode 
oder wegen 
1 
cos 2p cos ny = 5 cos (n—2)p + > cos (n + 2) 9, 
(7.8) 
1 1 
sin 2 @ sin ny =~ cos (n— 2) py —-~> cos (n + 2) p 
auch 
1 Yr 2n— 3 , 22 
X [4° Att, + ot A* at oP (Mea — Ta + * 7 * Hn) 
i. 2 2) eee 2+ 2 \ 
+g (Hes + ~ fhe g ite - bs na) cosng=0. (7.9) 


Damit (7.9) fiir alle @ erfillt ist, muB jedes Glied der Summe > fiir sich verschwinden, d. h. (7.9) 
zerfallt in ebenso viele Gleichungen, wie Glieder in der Reihe (7.6) verwendet werden. Jedes Glied 
IT, ist dabei mit den ,,benachbarten“ Gliedern //,,__. und //,,,2 gekoppelt. Da jedoch der Kopplungs- 
faktor 7? mit dem Quadrat des Radienverhaltnisses abnimmt, kann bei kleinem Radienverh4ltnis 
die Kopplung vernachlassigt werden. 


Die Differentialgleichung fiir /7,, lautet dann 
AS AP TTS A eT (7.10) 
und ihre bereits auf abklingende Funktionen beschrinkte Lésung 


= Gay at Cre Hn @) (7.11) 


hat die hinreichend einfache Form, die im Hinblick auf Anzahl und Gestalt der Randbedingungen 
notwendig erscheint, um die numerische Rechenarbeit in angemessenen Grenzen zu halten. 
Withum' hat bei der Behandlung eines wesentlich einfacheren Randwertproblems mit Hilfe der 
Stérungsrechnung eine strengere Lisung des Systems (7.9) entwickelt, die grundsatzlich auch hier 
angewendet werden kann. Allerdings ergibt sich dann jedes Glied /7,, der Reihe (7.6) selbst als 
Funktionsreihe 


IFe= ye + Hay + Hae + es 
so daB der Umfang numerischer Rechnungen im Vergleich zur einfachen Lésung (7.11) erheblich 
ansteigt. 


Bezeichnet R, bzw. J, den Real- bzw. Imaginarteil der in (7.11) enthaltenen Hankelschen 
Funktion erster Art ; 


HY (Vive) =R, +i Jn, 
1 siehe FuRnote 4, Seite 293. 


BN a a a Ue a A 
WN » sy . x 
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so ergeben sich aus der komplexen Funktion 
KS 1 
GW. iK yr |(4na —i Byx) Ga + (Ang—t Bus) (Ry + i J,)| c08 np 


durch Zerlegen in Real- und Imaginarteil die Verschiebung 


A, 1 — 
Ve= | oe + An, R, + Bi, In| cos ny = W, cosng (7.12) 
und die Spannungsfunktion 
B oe 
ky f= ae + Bro Ra— Ans J,| cosng = Ky? F, cosngy (7.13) 


und hieraus schlieBlich alle tibrigen Verschiebungen und SchnittgréBen mittels (5.21) bis (5.24), 
jeweils in der Form von Produkten aus sin n @ oder cos ny und einem nur von o abhangigen Faktor, 
der wie in (7.12) und (7.13) durch Querstrich gekennzeichnet sei. 

Die Verschiebungen U,, und V,, sind infolge der Koeffizienten sin? y, cos? y und sin 2 ¢ in (5.21) 
nicht allein mit W,, sondern auch mit W,,_. und W,,,2 gekoppelt. Diese Glieder entsprechen den 
Koppelgliedern in (7.9) und werden wie die letzteren vernachlassigt. Es gelten dann die Gleichungen 


62,—=U,+> Wn; a fon = — (Vi J U,) + 


1 / 7 
aVron — Oo (U,— Ve) SE Ve 


oder nach Kirzen durch sin ng baw. cosng 


ee a = Ih e5 otis , i= 
08, = ay (—H Ee ete Un + > Wi.» 


ee et ey litre. tS n i 
lon = pa (Pn — Pat ue Fil= 9 On oe Aarne 


a La, i, ie eS ay 1 = 
*j= pr te)n(> he e on BEAL” Tih gree 


Hieraus berechnet man leicht U,, und V,,, wobei fiir die aus der Hankelschen Funktion abstammen- 


den Glieder zweckmaBigerweise die aus (7.10) abgeleitete Beziehung E H y W,, =2aA*F, benatzt 
wird. Man beachte, daB in (5.21) und den daraus entwickelten Gleichungen U = U* und V = V* 
ist. Die Verschiebungen und SchnittgréBen lassen sich jetzt alle in der Form 


OG, tee ee, + 6, A,, + fs Baa) 


_ darstellen. Die Koeffizienten 6, Gr Be Om ip entnimmt man Tabelle 2. Die Ableitungen der Zylin- 


derfunktionen R,, und J, nach 0 bestimmt man mit Hilfe der Formeln 


R,=—7 R,+ Ov Roi, In = =F In + OY Ini: 


2 A, 2 (2) 
Ro — R, + er j,—— teal Jn = a a eis leer In—1> 


Reg i Jn = Vi (VLOr 0). 


Die Funktionen n-ter und (n — 1)-ter Ordnung ergeben sich nach bekannten Beziehungen aus den 


Funktionen nullter und erster Ordnung. 


lye Darts ; 
Fir n = 0 verwendet man, abweichend von (7.11), statt C,, on die Funktion i By, In g und die 
davon abgeleiteten GréBen 


Bee Bo 

Cas 27 6? 9 ” 
ha B 1 Bu 
Ne= KR ge” Ng Ky? Go 


22 
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Tabelle 2 
Ew Os ya ea z a 
U Ok ee d+) m— 2 n Q se Ks a 
1+e a sare ae —n-1 aL eet ee J Ww R 
Y | ae"? Farne@—? ae Bea 
J | 
W cos np (ont 0 Ky Jn 
W' cos np ners 0 R. Af 
1 a ; et nt , 
- Re oat. 0 —n(n+1)e@ : ie tae Ri, Fake. 
No = cos np 0 n(n + 1) 9" =” eee 
n 5 ‘ —n-1 he 1 Jf R’ 1 R 
Nig cay 0 —(n+l1)eo a8 See: i ae - : 
Q, ey coeny 0 0 In — Ri, 
ies a 
Gis 0 — J R, 
Qo 2Ky @e wees 0 n : 7 
rales ae ye COSND (1—p) n(n+1) 07 0 OP? In+(1—p) Rf —p 0? v? R+C-) Jn 
9 | 9 <5 a cosny |—(l—p)n(n+lor 0 @2 y? J, —(1—p) RY - 6? v? Rya-(1-p) J, 
a(l — pL) Lies ds at & , 1 ¢ 1 
eee n aed Pe == UR — ees 
M,9 2K@vo  * (n+ 1)e 0 at n Sop 


8. Die expliziten Randbedingungen bei kleinem Radienverhiltnis y. Setzt man die aus Grund- 
spannungszustand und Randstérung insgesamt hervorgehenden ZustandsgréBen 


w=Ww-+ > wu, cosng, W=W+ DW, cosnp usw. (8.1) 


in die Randbedingungen (6.29) bis (6.36) ein und verwendet man zur weiteren Umformung die 
goniometrischen Formeln (7.8) sowie 


sin 2 y cosng = 5 sin (n+. 2)9—> sin (n — 2), 
1 (8.2) 


cos2ysinng = =y sin (n + 2)o+> sin (n — 2) 
und die Reihenentwicklung 
1 
52 


Cy =y(1—2? sin’) | Sy +7 Psinty +3 Seino te, 


so erhalt man bei Vernachlassigung von Gliedern, die y? und héhere Potenzen von y als Faktoren 
enthalten, das folgende Gleichungssystem: 


0=—a> sin 29 —7 (ln —2—itn 2) sinng + 30, sinng—V*—SVssinng, (8.3) | 
0=%+ Yu, cosny + — ¥ (Pn42—Tn—2) cos ny + + W (1 —cos 29) 


+3 >) W, cos n@ —+ > (Wn42 + Wr—2) cos ny — Vey W,, cos ng, (8.4) 
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0=6+ > w, cosng— U*— > (ih cosnp — W = (1 —cos 2.9) 
Vv = —— — 
ay = W,, cos ng oe (Was2-+ Wa—2) cos ng, (8.5) 


0 = Jw, cosng—— > [(n + 2) H,42—(n—2) W,_»] cosng+ YWieosng, (8.6) 


0= 2 M,, cos ng + ay 2 (mz y (r42) — Mx y(n—2)) cosng— >’ M,, cosng, (8.7) 
ave por % x ss a 
0 = — = (n,,— n,) sin 2 —= a [(2,, — Ny)n—2 — (n, — Ny)n42] sin np 
+ Diligen sinng + N,» + Sea sinn@, (8.8) 


0 = > dn cosng + > (Gon+2) — Fy (n—2)) cosn.gy + ii, (1 — cos 2 ¢) 


+ > Dd) Ngn COS NY + D> (Mx (n+-2) + Mx (n—2)) Cos n Y — Ness Se, N,,, COS ng, (8.9) 


~ = Vv —. sy — TA 
0=n, +) n,, cosng + ae > (7x y (n £2) — Mx y (n+2)) cosng— »'Q,, cosnp 


v 


ea: N, (1 — cos 2 g) —+ > N,,, cos ng+ - SS (N, (+2) = N, (n—2)) cosng. (8.10) 


Die Form der in Abschnitt 4 angegebenen Grund-ZustandsgréBen — von g unabhangige und 
mit sin 2 » bzw. cos 2 m veranderliche Glieder — verlangt, daB die RandstérungsgréBen mindestens 
die ersten zwei Glieder der Summen 5) in (8.1), also wy + w, cos 2 y, W, + W, cos 2 p usw. ent- 
halten und die Bedingungen (8.3) bis (8.10) mindestens fiir n = 0 und n = 2 erfillen. Um die 
Genauigkeit der Lésung zu erhéhen und die Konvergenz der Reihen beurteilen zu kénnen, wird 
man auch w,cos4qg, W,cos 4g usw. in die Rechnung einbeziehen und die Randbedingungen 
auch fiir n = 4 erfiillen. Da die von y unabhangigen Spannungszustande keine Verschiebungen 


v baw. V* und keine Schubkrafte n,, bzw. N,, liefern, sind die Bedingungen (8.3) und (8.8) fiir 


n = 0 automatisch erfiillt und es verbleiben in diesem Falle noch sechs Bedingungen: 

a, es ae a 2 eee (L—$)— Wo. (8.11) 

6= ae a Ei, 2 EF (2p) —0,— 4 ale 5 2+ W,, (6.12) 

0= w+ Wo — > we, (8.13) 

0S meg ++ mae — M05 (8.14) 

0= Got Gat yet Fo — 4 tea PA— Myo, (8.15) 
0=F- t Meo 4 5 Ny ye =O0—-5 No ; iN pds (8.16) 

- Fur n= 2 gelten acht Bedingungen 

ee 0) et | @.10 

Can 2a See ” i, — 2 i) —W,+F%— FM, = 24G18) 

. 22" 
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0= w+ W,—vw,, (8.19) 
0= Mig, — M+ > Mya, (8.20) 
= . = pa =e Yy> = ; 
0= 5 mE 5 (n, Ny)o a Nye ae ar oe N92 a8 Ah (n, Ny) a> (8 21) 
eS = A = i eS 
0 = deo 2 re 5 Ny ’ Ny, Ie fi — Nig Wyss gee? (8.22) 
= a2 = D ie 
0= Nee r2 . pA 5 Nes 5 N,o 5 Ney | 4 N,4> (8.23) 


=m, — F252) 4+5 Fa (IZ) tg M—g tame Cam 
fiir n= 4. ebenfalls acht Bedingungen 
0=—7u,+%,—VZ, : (8.25) 
0=4,—7%+ 7 u,— 7 w,—W,, (8.26) 
0=w,—Uj—5W,4+—M,, (8.27) 
0=m+>u,+%, (8.28) 
0 = m1— > My,— M4, : (8.29) 
0=— = (,— 7p + Mya t Noa (8.30) 
0 = Ga— Gyn + > Mea — | ea — Nes (8.31) 
0S yO, A Nt (8.32) 


Zur Befriedigung der Bedingungen (8.11) bis (8.32) stehen die 22 Konstanten 


Gin boys Migs aus der Lésung gy 
@g9> bg45 Gao, Bas aus der Lésung g, fir Rohr TI, 
O41, 544, Ago, by, aus der Lésung g, 


Boy, Aga, Bog aus der Lésung G, } 
A,,, Byy, Ags, By. aus der Lésung G, fir Rohr IT 
Ajy;, By,, Ags, By aus der Lésung G, 
zur Verfiigung. do, tritt nur in der Verschiebung uy auf und beschreibt eine formanderungs- und 


spannungsfreie Verschiebung des Rohres I parallel zu seiner Achse, die durch die radiale Auf- 
weitung des Rohres IT herbeigefiihrt wird. 


9. Numerische Rechnungen. Hat eine Rohrabzweigung beispielsweise die in Abb. 20 dar- 
gestellten GréBenverhaltnisse mit den Parametern 


P02 = d. h. a oe 

ge eee d. h. a= 5,45 h, 

a= 16, d. h. A = 9,69 HT, 
so kann wegen des kleinen Radienverhdltnisses zur Bestimmung der Konstanten’ G1, by, usw. die 
in Abschnitt 8 entwickelte Form der Randbedingungen verwendet werden. In Tabelle 3 sind die 


mit obigen Parametern berechneten Konstanten als Vielfache der GréBe P = pal/E angegeben, 
und es sind die Konstanten der dreigliedrigen Ansitze g, -- 82 + g, und G+ G,+ G, den ent- 
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sprechenden Konstanten der zweigliedrigen Ansdtze g, + g, und Gy + G, gegeniibergestellt. Die 
letzteren lésen das Gleichungssystem (8.11) bis (8.24), nachdem man dort alle Glieder mit dem 


Index 4 gestrichen hat. Den Bedingungen (8.25) bis (8.32) 


sind die zweigliedrigen Ansatze nicht zu unterwerfen. 
Die Zustandsgréfen im Ubergangsgebiet zwischen 
den Rohren (9 = 1, §€—€&, = 0) einschlieBlich der 


Spannungen stellen sich simtlich in den Formen 


ky + ky cos 2 ym + ky cos 


oder 


k, sin2y + k, sin 4p 


dar, in der zweiten Form diejenigen, die einen von @ 
unabhangigen Summanden k, nicht enthalten. Die Ko- 
effizienten ky, k,, k, entnimmt man Tabelle 4. Die 
Normalspannungen resultieren aus Membranspannungen 


und Biegespannungen: 


4 p 


eas bf ba 
O;, a h a h2 ? Oy = Th. ae h ? 
y= N, ob 6M, No 6 Mo 
isthe Fe H? ’ Oem =H, Abb. 20. 
Tabelle 3 
3-gliedriger 2-gliedriger 3-gliedriger 2-gliedriger 
Ansatz Ansatz Ansatz Ansatz 
44/P 10,6247 10,6239 p 5,3051 5,2944 
oi/P 13,9131 13,9111 Ayalp 78,8088 78,7900 
Ay /P 62,1408 62,1342 BoolP 11,8401 11,8397 
a, ,/P 7,2436 7,2478 A,,/P —99,1044 | 99,1404 
b, 1/p 4,6872 4,703] lp 13,3323 13,2506 
22/P — 0,00224 — 0,0304 Ba>/p 3,3719 3,3289 
44/6 0,13469 Aq,/p — 3,5144 
bailp | — 0,02985 Bulp | — 2,1368 
a5/P 0,08132 | 13/P — 0,03945 
baolP 0,11537 Byo/B 0,04721 
Tabelle 4 
Pp 63,370 —2,204 0,039 U* |p 8,379 —6,935 —0,055 
ae 0 12,400 —0,013 V*[p- 0 6,059 0,096 
w/p 15,258 —14,07 0,216 Wp 66,216 —5 122 0,144 
w’/p 9,865 0,397 —1,122 W' |p —8,468 —0,437 —0,284 
n,/pa 0,500 —0,191 —0,014 N,/|pa 0,810 0,241 0,008 
ny/pa 2,948 1,908 0,026 Nelpa 7,980 4,603 0,180 
a 0 —0,942 —0,013 rglP a > —U, 
ae 0,750 0,309 —0,002 Q,/pa 0,337 —0,135 0,092 
qy/pa 0 —0,020 —0,004 Qp[pa 0c —0,239 0,082, 
m,/pa —0,1417 —0,0122 0,0024 M,/pa* —0,1419 —0,0119 Pea 
my/pa? | —0,0425 0,0239 | —0,0011 M,|pa —0,1398 0,2265 —0,0298 
mz y/pa?* 0 —0,0005 0,0025 M,/pa 0 —0,0825 i 
Spannungen in der inneren Oberflache 
a 20 0,043 
—22,56 —3,21 0,354 Gis oe | 1,63 | 0, | : 
a6 8,50 | 14,65 | —0,056 OglP 12,30 14,02 —0,320 
Spannungen in der auBeren Oberfliche 
0,74 —0,011 
27,96 1,13 —0,506 o,|p 4,71 | 5 | 
fa | 23,64 | 6,14 | 0,336 OylP 18,60 3,82 1,020 
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Verschiebung U* Verschiebung V* 


Verschiebung 0 


Verschiebung 70 


PaO is 


Verschiebung W 


Verschiebung u 


Mabstab: 


Mabstab: 2op 
203 ee 
——— 
Abb. 21. Verschiebungen der Schnittkurvenpunkte Abb. 22. Verschiebungen der Schnittkurvenpunkte 
in bezug auf das Koordinatensyxtem x, y,2 der Rohres I. in bezug auf das Koordinatensystem r,g,Z des Rohres II. 
Normalkratt Np Schubkraft Nog Normalkratt Tu, Schubkratt riary 


? da ae a 


| 


Normalkratt My Normatkraft roy 


——_-———_-—~ 


Mabstab: 

5pa. 
———— 
Abb. 23. Membrankrafte im Rohr I, 


Abb, 24, Membrankrifte im Rohr II. 
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Das positive (negative) Vorzeichen gilt fiir die Spannungen in der inneren (auBeren) Oberflache 
der Rohre. Der Verlauf der ZustandsgréBen langs der Schnittkurve ist in den Abb. 21 bis 28 dar- 
- gestellt. 
Die Materialanstrengung, die der einfachen Berechnung wegen als Differenz zwischen gréBter 
und kleinster Hauptspannung bestimmt wurde, erreicht ihr Maximum in der inneren Oberflache 
des abzweigenden Rohres, in dem in der Langsebene liegenden Punkte des unteren Randes (€ = 1/y, 


gy =9). Dort ist o, = — 25,4 p, oy = 23,1 p, o, = — p, also die Anstrengung o, = Craton 
= 48,5 p. Der ungestérte Grundspannungszustand in der inneren Oberflache des Rohres I mit 
ee 2:73 p, Gy = 5,45 p und G, = — p bewirkt die Anstrengung 6, = 6,45 p. Die Randstérung 


steigert also die Anstrengung auf das 7,5-fache; dieser Effekt beruht jeweils etwa zur Halfte auf 
der Ovalisierung des Miindungsloches, die zu einer betrachtlichen Steigerung der Umfangsspannung 
Oy fuhrt, und auf der Randeinspannung, welche die Achsialspannung o, stark negativ macht. 

Der Beitrag der Koeffizienten k, zu den Maximalwerten der Spannungen ist durchweg vernach- 
lassigbar klein. Da ferner die ky und k, praktisch unabhangig davon sind, ob mit dreigliedrigen 


oder zweigliedrigen Ansatzen gerechnet wird — siehe die Ubereinstimmung der entsprechenden 
Konstanten in Tabelle 3 — kénnen vorbehaltlich der Bestatigung durch weitere numerische Rech- 
Biegemoment Mp Querkratt 9, Biegemoment M, Querkratt (7, 
> = : - - > - . - = 


te aN, 


p= ° 


od Biegemoment My Querkratt g~0 


Biegemoment Mey guerkraft ky 


MaBstab: 


apa, Mabstab: 
a O2pa, 
Abb, 25. Biegemoment und Querkrafte im Rohr I. Abb. 26. Biegemoment und Querkrafte im Riohr II. 


nungen die zweigliedrigen Ansadtze gy + g, und G, + G, bei hinreichend kleinem Radienverhiltnis 
als befriedigende Lésungen betrachtet werden. Der hierbei zu beachtende obere Grenzwert des 
Radienverhdltnisses kénnte durch Vergleich mit Lésungen aus der erwahnten Stérungsrechnung 
bestimmt werden. Der Einsatz von Rechenautomaten erscheint hierzu unerlaBlich. 


10. Zusammenfassung. Der Spannungszustand im Ubergangsgebiet einer durch Innendruck 
belasteten, rechtwinkligen Rohrabzweigung — eines im Maschinen- und Apparatebau haufig ver- 
wendeten Konstruktionselementes — unterliegt verhaltnismaSig komplizierten Randbedingungen, 
die mit den klassischen Hilfsmitteln der Biegetheorie der Kreiszylinderschale im allgemeinen nur 
in einzelnen Punkten auf der Schnittkurve der beiden Rohrmittelflachen erfiillt werden kénnen. 


Unterschreitet das Radienverhaltnis a/A (a Halbmesser des abzweigenden, A Halbmesser des 
durchlaufenden Rohres) den Wert 0,5, so lassen sich die Randbedingungen soweit vereinfachen, daB 
ihre Erfillung in allen Punkten der Schnittkurve durch eine geschlossene Lésung des Randwertpro- 
blems méglich wird. Dies gelingt hauptsachlich durch Verwendung eines polaren Koordinatensystems 
auf der Mittelflache des durchlaufenden Rohres und einer vereinfachten Schalenbiegetheorie von 
Green und Zerna. Die Lisung hat zunachst die Gestalt unendlicher Reihen; man kann sie jedoch 
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in der inneren Obertltiche in der auBeren Obertliche 


ee - 


Mabstab: 
20p 


Abb. 27. Normalspannungen im Rohr I am Ubergang zum Rohr II. 


in der inneren Oberfliiche Inder ese Obertiche _#| 


Abb. 28: Normalspannungen im Rohr II am Ubergang zum Rohr I. 


bei kleinem Radienverhdltnis auf die ersten zwei oder drei Glieder der Reihen beschranken. Der 
Kern der numerischen Rechnungen besteht dann in der Bestimmung von 14 oder 22 Konstanten, 
mit denen schlieBlich samtliche Verschiebungen, inneren Krafte und Momente und Spannungen 
im Ubergangsgebiet der Rohrabzweigung berechnet werden kénnen. Ein spezielles Zahlenbeispiel 
ergibt eine Steigerung der Materialanstrengung im Ubergangsgebiet im Vergleich zur Anstrengung 
im ungestérten Grundspannungszustand um den Faktor 7,5. 


(Eingegangen am 16. August 1960.) 


Anschrift des Verfassers; Privatdozent Dr.-Ing. Willi Reidelbach, Sindelfingen/Wttbg., Eyachstr. 10, 
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Die Relation zwischen Temperatur und Geschwindigkeit in der Theorie 
kompressibler turbulenter Grenzschichten 


Von W. Szablewski 


1. Einleitung. Die Relation zwischen Temperatur und Geschwindigkeit spielt bekanntlich eine 
fundamentale Rolle in den Berechnungsverfahren kompressibler turbulenter Grenzschichten. Die 
bisher bekannt gewordenen Formulierungen dieser Relation, die die Temperaturverteilung als 
quadratische Funktion der Geschwindigkeitsverteilung ergeben, unterliegen jedoch gewissen Ein- 
schrankungen und Einwendungen. EF. R. van Driest! und jiingst J. Rotta? lassen in ihrer Herleitung 
der Relation den Hinflu8 der Dichteschwankungen auf den turbulenten Austausch unberiicksichtigt. 
Dichteschwankungen treten in kompressiblen turbulenten Grenzschichten notwendigerweise wegen 
der in diesen vorhandenen Dichteverteilungen auf, wie anschaulich sofort erkennbar ist, wenn man 
sich den ProzeB der turbulenten Vermischung als durch Verschieben ganzer Fliissigkeitsteilchen von 
einer Schicht in die andere erzeugt vorstellt. Auf die Bedeutung der Dichteschwankungen fiir die 
Theorie kompressibler turbulenter Grenzschichten haben bereits T. Y. Li? und H. T. Nagamutsu® 
hingewiesen*, jedoch tragen T. Y. Li und H. T. Nagamutsu dem Einflu8 der Dichteschwankungen 
nur in der Formulierung des Impulssatzes Rechnung. Beriicksichtigt werden die Dichteschwan- 
kungen auch in einem Artikel von G. B. Schubawer und C. M. Tchen®. Doch erscheint die dort 
gegebene Ableitung der Temperatur-Geschwindigkeitsrelation nicht zwingend. In der vorliegenden 
Untersuchung geben wir eine neue Formulierung der Temperatur-Geschwindigkeitsrelation unter 
Beriicksichtigung der Dichteschwankungen und unter sehr allgemeinen Annahmen. 
Wir machen im Folgenden die iibliche Voraussetzung der Grenzschichttheorie, da8 sowohl der 
_molekulare und der turbulente Austausch wie auch ihre Gradienten innerhalb der 
Strémungnurin Richtung normal zur Wand von Bedeutung sind. Weiterhin beschranken 
wir uns auf den wandnahen vollturbulenten Bereich der Grenzschicht, der gekennzeichnet ist 

-durch konstante Schubspannung und konstanten Energieaustausch. Die Strémung wird als im 
zeitlichen Mittel stationér und als ebene Strémung vorausgesetzt. 


2. Kontinuitat und Impulssatz. Es bedeuten im Folgenden: x bzw. y baw. z cartesische Koordi- 
naten langs bzw. normal zur Wand bzw. normal zur (x, y)-Ebene; u,v, w die entsprechenden 
Geschwindigkeitskomponenten; p den statischen Druck, 0 die Dichte, den Zahigkeitskoeffizienten, 
y = lo die kinematische Zahigkeit, T die absolute Temperatur, c, bzw. c, die spezifische Warme 
bei konstantem Druck bzw. konstantem Volumen, A den Warmeleitungskoeffizienten. Uberstreichen 
_ bedeutet Mittelwertbildung. 


1. Die Kontinuitatsgleichung lautet dann 


Lue, SNe (1) 


Spalten wir nach O. Reynolds die Momentanwerte der Strémung auf in zeitliche Mittelwerte und 
_ momentane Schwankungen: 


; w=ut+w, e=O0+0', (2) 
so folgt | 
4 au, ee _ A= eV), (3) 
ax Oy is ay 
2. Fiir den Impulssatz haben wir 
dou? , douv  _— dp di a ( po). (4) 
Ox ay dx ay \" ay 


1 E.R. van Driest, Journ. Aeron. Sci. 18 (1951), S. 145. 

2 J. Rotta, Z. Flugwiss. 7 (1959), S. 204. ; 

3 7. Y. Li and H. T. Nagamutsu, Journ. ees ioe (1951), S. 696. 

4 Vol. f W. Szablewski, Ing.-Arch. 20 (1952), 5S. 67. , 

e Ge ‘Be Schubaber we Cc. M. Tchen, High Speed Aerodynamics and Jet Propulsion, Bd. V (1959). 
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Mit 
Quv=v(dut+ou)+oQur +ug vw oun+ eur tug (5)? 
folgt 
dou? douv __ dp a du OT Tra, (Pag 2S 6 
oe et lea) t leer) tee) (6) 


3. Energiesatz. Der Energiesatz lautet bekanntlich 


div 


(ce T J we + \o 0 = div (A grad T) 


0 rr ae ee ou dv\ | ow ou 
+ gle 2 ev ! wolay + ee) woos Oz 
ou 


+ 
7) é 7) Dine ne ov ow 
my Hay | = | wo(2— 5 div) | wro(s. BS om 
fe ee 
me ja a | wolgt wy) Tem? Gs — 3 aiv 0) (7) 


Cp wird dabei als Konstante angesehen. 
Die Mittelwertbildung ergibt dann fiir den turbulenten Austausch normal zur Wand 


¢, Tour ce [o To +0 oe Ep! ef (8) 
und fiir 
(uc? =? 0") a BN aera 
osteo oe” zunachst ~aurov, 
wenn wir P 
vt wu ue 
2 Se 


fiir die Grenzschichtstrsmung annehmen. Mit 


w=wuv+ duu + u?, ev=ov+tov’+ve+ov 


folgt weiter 
wWovw [ev +u?)+ 2uve w)] + 20uu + Wo’ 
und bei Vernachlassigung von Gliedern kleinerer Gré8Benordnung? 
Wovw owv+ 2ouw +ao'v'. (9) 1 


Sehen wir im Folgenden, wie eingangs erwahnt, die Gradienten des Austausches nur in Richtung : 
normal zur Wand als wesentlich an, so erhalten wir aus (7) mit (8) und (9) 


zl T+ =$")o0] +2((,7+ +") a3] 
OTB du, av ercay tee aa 
= ay(¥ ay) + ale (Se +e) +H (2ay — a avo) tae (E+ 


+ ye % T’v') + ¢, T(—@' v) + u(— wv) + Cer | (10) fj 


Dabei haben wir in die Konvektionsglieder die an und fiir sich vernachlassigbaren Mittelwerte v2 |] 
wieder aufgenommen. | 

Wir fiihren jetzt mit der Wandtemperatur T, die Temperaturdifferenz 9 = T— T, ein (der | 
Index 0 kennzeichnet auch weiterhin die Wand). Wir kénnen dann bei Benutzung der Kontinuitats- 
gleichung (3) in der Formel (10) fiir den Energiesatz die Variable T' durch @ ersetzen. 


1 Der Herleitung von Formel (5) und (9) liegt zugrunde, daB wohl |u'| bzw. lo’| als klein gegen u bzw. @ 


anzusehen sind, jedoch nicht |v’| als klein gegen » angenommen werden kann; vel. W. Szablewski, a. a. O, i 
2 Zur Begriindung siehe FuBnote 1, ; =| 
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Nunmehr schranken wir unsere Betrachtung auf den wandnahen Bereich konstanter Schub- 


spannung und konstanten Energieaustausches (durch Molekiile und Turbulenzelemente) ein. Aus 


der Impulsgleichung (6) erhalten wir dann durch Integration 


T= Ma + 0(—w'v’) + u(—o'v’) (11) 
und aus dem Energiesatz (10) 
16($2) = 1-2 + 06, F 0) + 9-072) 
+ u(t Blt ee PE Lawn) tunfe a) 
+ 9uCwv) +5 Cor). (12) 


Die Glieder des molekularen Austausches in (12) lassen sich nach G. B. Schubauer und C. M. Tchen1 
auch auf die Form bringen 


0 (w+ 24+ wy? ov du\ , ov Ow 1 ou ov Ow 
He 2 |tu(use— og) te (ee — oe) a ay ee 
Das dritte Glied hier ist aus Griinden der Symmetrie (beide Richtungen der z-Achse sind vom Stand- 
punkt der Strémung nicht zu unterscheiden) gleich Null zu setzen. Da wir nur den molekularen 


Austausch normal zur Wand als wesentlich ansehen, kénnen wir auch pu = und uv - vernach- 
x x 


lassigen. Sehen wir noch das Glied 1/, v? als klein gegeniiber u2 + v? + w? an, so erhalten wir schlieB- 
lich fiir den Energiesatz 


mi) , OO = Guay. aw Y SD 
ue ) =I +0 (ed v’) +e, 9 (—o'v) 


ey] 
Cee sae ON Ge ee yy 
: + ea) te we) + es (13) 
_ diese Formel entspricht der Ausgangsform 
u? + v2 + w? oT oT 0 /wWtvt w? 
fey (eh ’9(5), Hs 3 ); (14) 


die in ttbersichtlicher Weise turbulenten und molekularen Austausch gegeniiberstellt. Gleichung (14) 
haben G. B. Schubauer und C. M. Tchen? angegeben. 

Wir betrachten weiterhin nur noch den vollturbulenten Teil des wandnahen Bereiches, in dem 
also der molekulare gegeniiber dem turbulenten Austausch vernachlassigbar ist. Aus (13) folgt 
dann unter Beachtung von (11) 


68 See eee pees Pa) = tt Oceans 
Jy ($2) = Bey =O 0) + 6, 8 (— ev) + wa—F CET). a 


_ Die unterstrichenen Glieder ergeben zusammen die Dissipation der Strémung. Das bemerkens- 
_ werte Ergebnis ist, daB bei Beriicksichtigung der Dichteschwankungen im Ausdruck fiir die Dissi- 
_ pation der Strémung zu dem Glied 1, u, das mit Einschlu8 der Dissipation der Hauptbewegung die 


der Hauptbewegung entzogene Energie zur Erhaltung der Nebenbewegung darstellt, das Glied 
(a?/2) ou hinzukommt, das als Diffusion kinetischer Energie mittels des turbulenten Massen- 
transports aufzufassen ist (vgl. hierzu auch Abschn. 4). 


4. Die Temperatur-Geschwindigkeitsrelation kompressibler turbulenter Grenzschichten. Wir legen 
im Folgenden die sehr allgemeinen Boussinesqschen Ansatze zugrunde 


—wy = &(x, y) i : (16a) 

rg Eeuy) a, (16b) 
potas 60 

—o'v = E¢(x, y) ae (16¢)3 


Der Austauschkoeffizient ¢(x, y) hat dabei die Dimension einer kinematischen Zahigkeit. 


1G. B. Schubauer und C. M. Tchen, a. a. O. 
2 G. B. Schubauer und C. M. Tchen, a. a. O. 
3 Vgl. hierzu W, Szablewski, a, a. O. ; 
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Wahrend (16a) und (16b) in der Theorie turbulenter Grenzschichten gelaufige Ansatze sind, 
ist zur Erlauterung des Ansatzes (16c) zu sagen: Die obigen Gleichungen entsprechen in Analogie 
zur kinetischen Gastheorie der Vorstellung, da8 Fliissigkeitseigenschaften von einer Schicht in 
die andere transportiert werden. Gleichung (16a) spricht das fiir die Geschwindigkeit, (16b) fir | 
die Temperatur aus. Da nun in einer kompressiblen Grenzschichtstrémung verschiedene Schichten _ 
auch verschiedene Dichten aufweisen, ist ein Transport der Dichte von einer Schicht in die andere || 
anzunehmen. Gleichung (16c) ist dann nur der (16a), (16b) entsprechende Ansatz. 

Der Faktor E ist ein Ma® dafiir, da8 beim turbulenten Transport fir die substantiellen Kigen- | 
schaften bekanntlich ein anderer Austauschkoeffizient in Erscheinung tritt als fiir die wahrend des | 
Transports Druckschwankungen ausgesetzte Geschwindigkeit. 


Wir erhalten dann nach (11) und (15): 


Impuls . 
/— du — 00 
T= 8 [0-5 + rus). (17) 
Energie | 
ad 20 8 " ui? 00 | 
(22) = Eee a LT) U Des iy (8)q 


Mittels der Gasgleichung fiir konstanten Druck (gema8 der Annahme der Grenzschichttheorie) 
OT QT, (19) | 


kénnen wir in den Gleichungen 9 durch T = i oe @ ersetzen. 


Fiihren wir noch mit der Schubspannungsgeschwindigkeit v, = |/T)/@) dimensionslose Koordi- 
naten 


Stey Rite =e 
7= ae Sige und 4 pe 
und die Parameter 
ad 
ats) 
Apes a ee (20) 


aes M. = 
Cy Qo Ty Vy 0 Vo—De,™ cp To (vy = ¢,/¢,) 


ein, so kommt: 


Impuls 
we Scr le bce Os te ax/én 
mir Boge 1) } 
Energie 
Seppe 2 i ee ON wo 
Ay — (vy — 1) M3 Eat Mes]. (22) 


Aus diesen Gleichungen gewinnt man die Temperatur-Geschwindigkeitsrelation in folgender Weise: 
Nach (21) ist 
€ 0w/8n é dy/én 
BER CONC Ee he Oe 
 1+x Ae % +x 


oder unter Beachtung von (22) 


é€ bo/én 1+ A, o—(y—1) M3 w?/2 


(23) | 


Mm ity — + 1+ (y—1) M3 o*/2 
Mit der Form (22) 
& Ox/én _ Ay—(y—1) Mio 
Mm (+x? 1+ (y~1) M3 w/2 
erhalt man durch Division der beiden letzten Gleichungen 
1 da/an 1A oe G1) Meee 
1 ae 0 oe 5 07/2 
E ayléq ‘t+ %) Ag Qn Mio a’ (24) 


? Der Parameter M), der von R.G. Deissler NACA TN 2138 (1! i ti i 
mit der Schubspannungsgeschwindigkeit v, gebildeten Machzabl, se CN a ley id aay eapy se Wl 
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Integration ergibt dann die neue Relation 


(l+7)®=Cj1+ A,o—(y—1) ma& 


oder ao ; ; ; (25) 
Ss =C)1+ A X— fy yt (4). 


Die Integrationskonstante C, die als abhangig von den Parametern Ay und M, anzusehen ist, ist 
im Rahmen der vorliegenden Untersuchung nicht weiter bestimmbar. Ihre Ermittlung wiirde die 
Integration der durch den molekularen Austausch vervollstandigten Differentialgleichungen von 
der Wand aus unter Beriicksichtigung der Randbedingung erfordern. Dabei ware dann noch eine 
Hypothese iiber den turbulenten Austausch, d. h. iiber die funktionale Abhiangigkeit des Austausch- 
koeffizienten ¢ von den mittleren ZustandsgréBen des Feldes, zu treffen. 


Wir wollen noch anmerken, daB man bei Vernachlassigung der Dichteschwankungen ‘die Glei- 
chungen: 


Impuls 
__ € 0a/én 
“we TEK" () 
Energie 
Ae ay Me os et Hen (b) 
9 Vek % 
: erhalt, aus denen die Relation 
2 
E x = const + Ayo — (y—1) MBS (c) 


folgt!. Die Integrationskonstante tritt also hier additiv auf, wahrend sie in (25) als multiplikativer 
-Faktor erscheint. 

Die Relation (25) enthalt eine bemerkenswerte Folgerung. Fiir den Grenzwert 1 + y + 0, der 
mit gegen Null strebender Tangente bei nicht vernachlassigbarer Dissipation (M, + 0) zu erwarten 
ist (vgl. Abschn. 6), folgt aus (25) als entsprechender Grenzwert der Geschwindigkeit mit gegen 

_Null strebender Tangente der aus der Gleichung 


sich ergebende Wert 
1 
(y—)) My 


Ay + VA? + 2(y — 1) M2}. (26) 


Om 


Diese Formel gibt bereits einen ersten Uberblick tiber die Staffelung der Geschwindigkeitsprofile 
nach den Parameterwerten A) und Mo. 


5. Vergleich mit Messungen. Zum Vergleich ziehen wir eine neuere amerikanische Messung® 
heran. Die MeBpunkte haben wir dabei, da uns die Originalarbeit nicht zur Verfiigung stand, einer 
Arbeit von J. Rotta? entnommen. 
| Der von uns erhaltenen Relation (25) kann man auch die Form geben 


(00) 


2 
Ig (+) = lg [1 + 4po—(y —1) MBS] + Zig C. (27) 


In der entsp.cechenden Auftragung der einzelnen MeBreihen (vgl. Abb. 1) miiBten also die im voll- 
‘turbulenten wandnahen Bereich gelegenen MeSpunkte auf parallelen Geraden gelegen sein. 


? Aus der Neigung dieser Geraden miBte ferner das Ubertragungsverhaltnis E (vgl. Abschn. 4) 
_ablesbar sein! Die in Abb. 1 wiedergegebenen vier MeBreihen verschiedener Machzahlen M,, (Mx 


~Machzahl der Strémung auBerhalb der Grenzschicht) darf man wohl als eine Bestatigung der 
Relation (27) ansehen. Aus der Auftragung ergibt sich weiter, dab das Ubertragungsverhaltnis nur 
‘wenig von E = 1 verschieden ist; der von H. Ludwieg* erhaltene Wert E = 1,1 scheint sich der 


‘Messung am besten anzupassen. 


1 Vel. J. Rotta, a. a. O. ' 
2 R-K. Lobb, E. M. Winkler, J. Persh, Experimental investigation of turbulent boundary layers in hyper- 
‘sonic flow. U.S. Naval Ordnance Laboratory, NAVORD Rep. 3880 (1955). 
_ 8 J, Rotta, a.a. O. 
4 H. Ludwieg, Z. Flugw. 4 (1956), S. 73. 
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—_— 


In der graphischen Darstellung der Abb. 2 wird die Relation (25) durch ein Strahlenbiischel mit 
dem Nullpunkt als Ursprung wiedergegeben. Auch diese Darstellung, der wir E = 1,1 zugrunde- | 


legten, scheint unsere Relation zu bestatigen. 
SchlieBlich wollen wir noch demonstrieren, daB die alte Relation [Gl. (c), Abschn. 4] 


1 wo 
l+y=Zll+ 40—¢~—) Moz + konst. | 
die Verhaltnisse weniger gut wiedergibt. In der entsprechenden Auftragung der Abb. 3 miiBten nach 
der obigen Relation die MeBreihen durch parallele Geraden der Neigung 1/E wiedergegeben werden. 
-lg(*x) 
06 " < 
7 
wie 
05 + Sige 
A 
Uf y 
Ml “s : 
a4 sy 
Ze 
Wl Ya | 
V f EZ, We, 
be Le Oy ae 
i } a L AE i a 
Wie > Va | 
iy pp 
ge Ve x — f= : | 
his J ae ag Theorie 
3 Ae y, o Nr 1=M,=5 | 
ay MEDIAS 6 fey 58 Messung R.K. Lobb, 
Say c L ag 68 E.M, Winkler, IPersh, 
Mf fey ee 1 NAVORD Rep.3880 (1955) 
/ 7 84 
4 ee is) Z Ys 
Za wa YY -\y[74A,o-(y-4)m 072| 
yl Oo G2 Q3 OF O5 be 
Abb. 1. 
i a BO a a ee 8 
1A, o-(y-1) My wife a 
|_| 
gs 
£ 
iF r 7 rr 92 
eres = a3 
8 . ) 
1 RG Tae ES i en 


05 
eo T _ 0, 6 
‘ : | Mae Sale 
we | Theorie (E=41) o 
, Bye oNn4 Moo=5 
F : 7 58 Messung R.K. Lobb, 08 


2 10 68 6. Winkler, XPersh, 
5 3, NAVORD Rep. 3880 (1955) 


13 82 
tedweas 
40 


j (t+) 
Abb, 2. 


“Y | | 
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Legen wir wieder E = 1,1 zugrunde, so ist aus Abb. 3 deutlich ersichtlich, da dies nicht der Fall ist. 
- Erwartungsgem4B treten dabei die Abweichungen vom theoretischen Verlauf bzw. der EinfluB der 

Dichteschwankungen um so starker in Erscheinung, je groBer M,, ist. Wir merken noch an, daB 
“Abb. 2 und Abb. 3 aquivalente Darstellungen der neuen und der alten Relation sind. 


Seah Ls 09 08 G7 06 G5 OF 93 G2 of 0 
Ay o-(y-1)lMyoie ead 
ya 
ae 7s aL) 
TE Cala Za 93 
ibe oalkee 
7 We ie? 
CRN nas ol 
ary Fi AE: OY 
7A 4 se? ? 
we 70 4 
Z 7 Via gv 
Zi a7 7e 7] 05 
Vanes YE 4 
fia 6 WA Be 
Za Ze WA 7 
“itis 7 A At gb 
a Taf Z 
(a) Sa wy 
s YA we 7 Je 
Be ASNT VA 
aaa er, marae L g7 
ae oe me ——— diltere Theorie . 
ae folie 7 o,7 oNr? Mi=5 
eae “_¥ ‘4 ahr gg MessungRK.Lobb, og 
Die ie oe ag 68 E.M.Winkler, J. Persh, 
YN x Speer Zz 9 NVAVORD Rep. 5880 (155) 
py ee de spats 82 
wali [mie ane 2 99 
th 
o ° | 
— — 10 
° E (14x) 
Abb. 3 


6. Geschwindigkeits- und Temperaturverteilungen. Die Gleichungen (21) und (22) erméglichen 

- bei Zugrundelegung einer Hypothese iiber den turbulenten Austausch die Berechnung der Geschwin- 
digkeits- und Temperaturverteilung. Als eine solche Hypothese bietet sich die Prandtlsche Mi- 
schungsweghypothese 


du € dw 
A = 72 vy? —— b . = SOF VE . 28 
eae y By ZW 5 en an (28) 


Einsetzen von (28) in die Gleichung (23) 


joey be yo — (— 1) Me w4/2 (29) 
a S47 1+ Y —1) Mj 0/2 


ergibt dann 


: E+] 
do __ PEL AD+4,0—-—¢—1) Mio] ** 
on a 1 V1 + —D M o?/2 
bzw. 
—Ing+D= a: Presi? eer (30) 
- [1 + A, o— (y—1) Mj w?/2] 28 


6 
mit der Integrationskonstanten D. 
Spezialfalle: 1) M, = 0 bzw. Dissipation vernachlassigbar. Nach (25) haben wir 


(l+ 4)" =C(1+ 40), (31) 
wahrend aus (30) folgt 
2E 
aie fh fal Woon 
Teeter [c# 4, 2 (Zin + D) + | (32) 


2) A, = 0 baw. warmeisolierte Wand. 
23 


setae ih VS 
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Gleichung (25) ergibt die Relation 


wo? / 
(+ 0%= C]1—y—1) M35], (33) 
und (30) ergibt mit = M, @ 
1 : 
ze (1 ie: pee at Jl + 03 
ema 2) Vela me M, gr da, - (34) 


(1 — o{) 2E 


Formel (34) beschreibt (siehe Abb. 4) einen monotonen Verlauf der Kurve w tiber 7 mit gegen Nall 
strebender Tangente und dem Grenzwert 


im allgemeinen Fall siehe Gleichung (26). 


Y, 
10 
| 
Q5 
| Sea a 0-2 doy, 
ae eeu _ 

0 5 10 

Abb. 4. 


Das in Abschn. 4 iiber die Ermittlung-der Integrationskonstanten C Gesagte gilt auch fiir die 
Integrationskonstante D. Wir fiigen noch hinzu: Die bisher formulierten Hypothesen tiber den 
turbulenten Austausch enthalten bekanntlich — man denke an das logarithmische Gesetz — zwei 
empirische Koeffizienten, deren Charakter insbes. bei kompressiblen turbulenten Grenzschichten 
noch weitgehend ungeklart ist. So zeigen z. B. Experimente! ein Anwachsen der laminaren Unter- 
schicht gegeniiber inkompressibler Strémung an. 


7. Zusammenfassung. Beriicksichtigt man die Dichteschwankungen kompressibler turbulenter 
Grenzschichten, so wird der Energiesatz insbes. insofern modifiziert, als im Ausdruck fir die Dissi- 
pation der Strémung zu der der Hauptbewegung durch die Nebenbewegung entzogenen Energie 
eine durch die Dichteschwankungen bedingte Diffusion kinetischer Energie hinzutritt. Bei ent- 
sprechender Beriicksichtigung der Dichteschwankungen auch im Kontinuitats- und Impulssatz 
erhalt man eine neue Formulierung der Relation zwischen Temperatur und Geschwindigkeit, aus | 
der sich u. a. Grenzwerte der Geschwindigkeitsverteilungen ergeben. Der Vergleich mit Messungen 
ergibt eine Bestatigung der Theorie. 


(Eingegangen am 1. September 1960.) 
Anschrift des Verfassers: Prof. Dr. W. Szablewski, Berlin W 8, Mohrenstr. 39. 


1 Vgl. F. K. Hill, Journ. Aeron. Sci. 23 (1956), S. 35. 
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Das Auswuchten elastischer Rotoren, ein Problem der Strukturanalyse 


Von E. Hiibner 


I. Einleitung. Das Auswuchten elastischer Rotoren hat in den letzten Jahren zunehmende 
Bedeutung erlangt und ist Gegenstand zahlreicher Arbeiten geworden!. Die meisten Untersuchungen 
beschrankten sich aber auf homogene Gebilde und leiteten in rein qualitativer Betrachtungsweise 
die Bedingungen her, die beim Auswuchten elastischer Rotoren grundsatzlich zu beachten sind. 
So wertvoll die hierbei gewonnenen Erkenntnisse auch waren, so ermiglichten sie doch keineswegs, 
die zur Kompensation der Unwuchten in einem Rotor allgemeinster Form erforderlichen Tarier- 
massen treffsicher nach GréBe und Lage zu bestimmen. Der Auswuchtvorgang gestaltet sich 
weiterhin recht aufwendig und fiihrt nicht immer zum gewiinschten Ergebnis. 

Von Meldahl? wurde erstmals ein Weg zur quantitativen Lésung der Aufgabe gewiesen. Jedoch 
bereitet die praktische Durchfiihrung des Verfahrens — Auswuchten des Laufers bei seinen kritischen 
Drehzahlen — nicht unbetrachtliche Schwierigkeiten und ist nicht ohne Gefahr fiir Maschine und 
Rotor anwendbar. 

Die nachstehenden Ausfiihrungen sollen zeigen, wie die Aufgabenstellung auf ein wohldefiniertes 
Problem der Strukturanalyse zuriickgefiihrt werden kann, das sich mit Hilfe der klassischen Berech- 
nungsverfahren lésen JaBt. Allerdings erfordert der Lésungsgang, zu dessen Durchfiihrung die 
Higenlésungen des Rotors bis zu einer durch die geforderte Genauigkeit festgelegten Ordnungszahl 
notwendig sind, einen ziemlichen Aufwand an Berechnungsarbeit.. 

Das entwickelte Verfahren ist fiir starre und nachgiebige (isotrope) Lagerung sowie auch fiir 
mehrfach gelagerte Rotoren anwendbar und gewinnt dadurch fiir das betriebliche Auswuchten eine 
gewisse Bedeutung. 

Der meBtechnische Teil des Auswuchtvorganges beschrankt sich hierbei unter AusschluB jeglichen 
Experimentierens mit Zusatzmassen auf das Messen von Lagerkraften (Auswuchten in starren 


~ Lagern), bzw. auf das Messen der Schwingungswege der Lager (Auswuchten in nachgiebigen Lagern) 
- bei einigen diskreten, im Voraus passend zu wahlenden Drehzahlen, die entsprechend den Sicher- 


heitserfordernissen geniigend weit ab von den kritischen Drehzahlen liegen kénnen. Die zum Aus- 


_ gleich der Unwuchten in den frei wahlbaren Ausgleichsebenen anzuordnenden Tariermassen kénnen 


nach vollzogener Messung rechnerisch bestimmt werden. 


2. Grundbeziehungen. a) Ruhender Biegungstrager. Fir den ebenen Belastungsfall lat 
sich die Abhangigkeit zwischen der Durchbiegungsfunktion y(z) eines geraden stabférmigen Bie- 
gungstragers mit kontinuierlich verteilter Masse und Steifigkeit und der Belastungsfunktion p(z) 


_ unter der hier zulassigen Vernachlassigung des Schubspannungseinflusses in der bekannten quellen- 


maBigen Darstellung zu 


y(z) = JG, 6) pC) a, (La) 


mit G(z, ¢) als Greenscher Funktion, anschreiben. 


Wirken auf den Trager zwei verteilte ebene Belastungen p,(z) und p,(z) in zwei zu einander senk- 


rechten Ebenen, so ergeben sich die Durchbiegungen der elastischen Linie aufgrund des Uberlage- 


rungssatzes als geometrische Summe der in den einzelnen Belastungsebenen hervorgerufenen Biege- 


- pfeile y,(z) und y,(z). Da man die beiden verteilten Belastungen als Projektionen einer raumlich 


=) 


verteilten Belastung senkrecht zur Tragerachse auffassen kann, ist damit der Weg gewiesen, wie 


sich zu einer raumlich verteilten Belastung die Biegepfeile, die in ihrer Gesamtheit die raumliche 


j Durchbiegungsfunktion bestimmen, finden 148t. Die fiir die beiden Teilbelastungen eines Biege- 
. tragers mit rotationssymmetrischen Querschnitten geltenden Gleichungen 


yilz) = f Gz, 6) ps(C) ae , Yolz) = f G(z,C) pal) dC 
kénnen mit Hilfe der symbolisch-komplexen Schreibweise durch eine einzige Gleichung 
(2) = f G(s, 0) pd) dc (1b) 


1 Ausfihrliche Schrifttumshinweise finden sich beispielsweise bei H. Becker, ATM, V 8224—9 (1960), S. 76. 
2 A, Meldahl, Z. angew. Math. Mech. 34 (1954), S. 317. 
23* 
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ausgedriickt werden, in der die durch einen dachformigen Akzent gekennzeichneten und durch 


p(z) = p,(z) + ip2(z) als komplexe Belastungsfunktion, 
4¥(z) = ¥;(z) + iy.(z) als komplexe Durchbiegungsfunktion 


(mit i = y¥— 1) definierten GréBen die Belastung und die Durchbiegung des Tragers nach Betrag 
und Phase, bezogen auf eine feste Bezugsebene, bestimmen. 


b) Schwingender Biegungsstab. Die die Schwingungsbewegung eines kontinuierlichen stab- 
férmigen Biegungstragers, der durch eine verteilte, zeitlich harmonisch veranderliche Erregerkraft 


mit der Kreisfrequenz = /2 und der Amplitudenfunktion p(z) zu erzwungenen Schwingungen 
erregt wird, in quellenmaBiger Darstellung beschreibende Gleichung 


y(2) = f Gz, 6) [A o(6) (2) + pS) ae (Le) 


bringt den bekannten Sachverhalt zum Ausdruck, wonach man den Vorgang als quasistationaren 
Belastungsfall behandeln kann, wenn als Belastung des ruhend zu denkenden Biegetragers die 
Summe aus den Erregerkraften und den d’Alembertschen Tragheitskraften 


seared ts P*(2) =4 o(2) y®) + Pl) (1d) 
eingelubrt wird. 


Die Amplitudenfunktion y(z) kann aufgrund des Entwicklungssatzes mit "y(z) als den Eigen- 
funktionen, und k, als den Entwicklungskoeffizienten (auch Normalkoordinaten oder generalisierte 
Koordinaten genannt) als Summe 


ve) = 3 hy "ye 2a) 
angeschrieben werden. Die Entwicklungskoeffizienten sind dabei in bekannter Weise durch 


1 
= | pv) de 
k= 
if Ly (A,/A art 1) 
bestimmt, wobei A, die Eigenwerte und 
UM, = Jf "y?(z) (2) dz 
die generalisierten Massen des Gebildes darstellen. 


c) Umlaufender Biegungstrager. Man denke sich ein rotationssymmetrisches stabférmiges 
Gebilde durch Radialschnitte in infinitesimal diinne Scheiben von der Masse 


dm = 0(z) dz, 


(2b) 


mit 0(z) als Massenbelag, zerlegt. Als Folge von Werkstoff-Inhomogenitaten und kleiner Bearbei- 
tungsfehler werden die Schwerpunkte der Scheiben im allgemeinen nicht in die Rotationsachse 
fallen. Es werde angenommen, da die von der Gesamtheit der Schwerpunkte gebildete Schwerelinie 
des unverformten Gebildes eine zumindest abschnittsweise stetig verlaufende, raumlich gekriimmte 
Linie darstellt, die durch ihre Projektionen auf die beiden zueinander senkrecht stehenden kérper- 
festen Kraftebenen gegeben sei. Der Abstand des Schwerpunktes einer elementaren Scheibe von 
der Rotationsachse, die Exzentrizitat, ist dann durch zwei Projektionen bestimmt und kann in 
symbolisch-komplexer Schreibweise durch eine komplexe Zahl ausgedriickt werden. Diese ist eine 
Funktion der Lagekoordinate z und werde als komplexer Exzentrizitatsbelag ¢(z) bezeichnet. 

Das Gebilde laufe mit der Winkelgeschwindigkeit Q gleichférmig um. Unter der Wirkung der. 
bei der Rotation hervorgerufenen Fliehkrafthelastung werden Verformungen auftreten. Mit Q2— 


und der Belastungsfunktion 
p*(2) = 2 [¥(z) + e(2)] o(2) (3) 
ergibt sich unter Vernachlassigung der Kreiselwirkung die komplexe Durchbiegungsfunktion gemaB 


(1b) zu 
¥(2) =A J G(z,£) [H(S) + e(C)] of) de. 


Tragt das Gebilde zusatzlich m diskrete Massen M,; mit den Exzentrizitaten ¢@,, so kann die Be- 
lastungsfunktion durch : 


p(s) =2 (9G) + OC) 02) + B os(2) MCV) + &) 


Bs 
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bzw. 
Pe) =A(ol) + S046) mh) x6) +2(0e0) of) + Solo a1) () 
Te 7 ait 
ausgedriickt werden. Hierbei bedeuten die GréBen p;(z) sogenannte Nadelfunktionen (im englisch- 
sprachigen Schrifttum Dirac-delta-functions genannt) gemaB Abb. 1, die den Bedingungen 
g,(2)=0 — fiir ae a+, 
p;(2) = A fiir Bl z<2+6 


mit 6—> 0 und lim 46 = 1 geniigen. 
6-0 


(5) 


g(2)—— 


Ss 


7 


Abb. 1. Nadelfunktion 9(z) (Dirac-delta-functions). 


Die Belastungsfunktion p*(z) in (4) stimmt formal mit der in (1d) dargestellten Funktion iiber- 
ein. Aus diesem Grunde kann die gesuchte Durchbiegungsfunktion des umlaufenden Rotors als 
Linearkombination der dem Gebilde zugehérigen Eigenfunktionen analog zu (2a) durch 

He) = Fk, "y(@) | (6a) 
mit den komplexen Entwicklungskoeffizienten 


jf ic (z) e(z) + ~ pj(z) M; é) ry(z) dz 


| = (ApJ — 1) (ob) 
dargestellt werden. Mit den Setzungen 
u(z) = Q(z) e(2) (6c) 
als komplexer Unwuchtbelag und 


als komplexe Unwuchten! der Einzelmassen ergibt sich aus (6b) nach teilweiser Integration 
[1y(2)- BG) de + 3 y(es) iy 
{poner a 6 
fi A —2) Ge 


Das Integral im Zahler des Bruches auf der rechten Seite der Gleichung kann mit Hilfe der bekann- 
ten Quadraturformeln? (Trapezregel, Simpsonsche Formel oderGaufsche Formel) in eineSumme gemaB 


J *y(@) We) ds = S We "y(e) Ae) 


umgewandelt werden, in der die GréSen W; durch die Wahl des Integrationsverfahrens festgelegte 
,,Gewichte“ sind, die beispielsweise bei der Simpson-Formel die bekannte Folge h/3 4/3 h 2/3 h 
4/3 h---h/3, mit h als Intervallbreite, bilden. Damit geht (6e) in 
SW tye) Wei) + Dy(e;) aj 
eee J 6f 
fi, (JAI) © 
iiber. Die Gesamtheit der Entwicklungskoeffizienten laBt sich in Matrizenschreibweise zu einem 
Spaltenvektor zusammenziehen. Mit den Setzungen 
"y(21) | 
"y (22) 
"t= | ys) 


(7a) 


"Y@n) 
1 Entgegen der iiblichen Bezeichnungsweise sei hier unter Unwucht das Produkt Masse mal Achsabstand 


tanden. : 
ES Siche etwa F. A. Willers, Methoden der praktischen Analysis, 3. Auflage, S. 137 und S. 152, Berlin 1957. 
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als approximierender n-dimensionaler Eigenvektor r-ter Ordnung, 


0, W, (2) a uy 
Us W, U(Z.) alas iu, 

t= |e. =| eae) +a eal 
&) write) +6, 


als n-dimensionaler Spaltenvektor der Unwuchten! des Rotors erhalt man aus (6f) mit (7) 


=> 


1 see 
Sores oreo oe 


woraus sich schlieBlich mit der Diagonalmatrix 


F, 
F, 
F= 2 F, 5 (8b) 
deren Elemente durch 
1 
ae sea (pasde2oarsh (8c) 
bestimmt werden, und der Matrix der approximierenden EKigenvektoren 
X= Pere *e...) (8d) 
der Spaltenvektor der Entwicklungskoeffizienten endlich zu 
ke 
ae k, ; By 


errechnet. Die in den durch die gewahlten Teilungspunkte des Integrationsbereiches gelegten 
Radialebenen — sie seien MeSebenen genannt — herrschenden (komplexen) Durchbiegungen 
ergeben sich dann aus (6a): 


Die Gesamtheit dieser GréBen kann zu einem n-dimensionalen Spaltenvektor zusammengezogen 
werden, der sich mit (8d) und (8e) zu 


y (2) 
¥ (22) a : 
H=]ye%)|=XFE=HXRXU (9) 
= 5G) 
erg . 


3. Rotor auf starren Stiitzen. a) Auflagerkrafte. Die zwischen der Belastungsfunktion p(z) 
und der die Durchbiegungen beschreibenden Funktion y(z) gemaB (la) geltende Crundbeziehung 
der Elastizitatslehre kann mit Hilfe der sogenannten Steifigkeitsfunktion C(z,¢) auch durch 


P(2) = f C(z,C) y(C) de 
dargestellt werden. Unter Beachtung von (6a) erhalt man daraus mit den Setzungen 


"P(2) = f C(z,C) (6) dé (r= 1,2,3...) 


‘ : : : ; : : 
Da im allgemeinen nicht in allen MeBebenen Einzelmassen angeordnet sein werden, kénnen einzelne 
(oder alle) Glieder @; in ll verschwinden. 


= es A é 
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als den den Eigenformen zugeordneten Belastungsfunktionen, den Kigenform-Belastungsfunktionen, 


Pe) = > 8, pb). (10) 


Diese Gleichung driickt den bekannten Satz aus, daB cine beliebig vorgegebene Belastungsfunktion 
nach Higenform-Belastungsfunktionen entwickelt werden kann. 

Bezeichnet man die bei der Verformung eines Biegetragers in ciner Eigenform auf die Stiitzen 
wirksam werdenden Krafte — sie mégen, ebenso wie die Kigenwerte und Eigenfunktionen aus 
einer Schwingungsberechnung herriihrend als bekannt vorausgesetzt werden — mit 


aS rk eee = 10s TG 


wobei r die Ordnungszahl der zugehérigen Eigenfunktion und | die Anzahl der Stiitzen sei, so 


ergeben sich die bei einer vorgegebenen Verformung des Gebildes herrschenden Auflagerkrafte 
wegen (10) zu 


4; = 3k, "8; (==, Janae (11a) 
FaBt man die zu einer Eigenfunktion gehérenden Auflagerkrafte zu einer Spaltenmatrix 
"Sy 
iS 
GS) — ¥S5 


('s, 
zusammen, so kann man deren Gesamtheit zu der Matrix der Eigenform-Auflagerkrafte 
S = (SSG...) 


vereinigen. Bildet man noch den Spaltenvektor der zu einer gegebenen Verformung gehérenden 
Auflagerkrafte, der in transponierter Schreibweise durch 


{’ = (A, 4, ...4,) (11b) 
gegeben ist, so kénnen die Gleichungen (lla) zu einer Matrizengleichung 
| Gok 
zusammengezogen werden, die schlieBlich mit (8e) die endgiiltige Form 

T=SFxX i (11c) 


annimmt. % 

Die Auflagerkrafte erscheinen als lineare Funktionen der Unwuchten U;. Bestimmt man daher 
mittels Messung die an den Auflagern eines mit der Winkelgeschwindigkeit 2 umlaufenden Rotors 
wirksam werdenden Auflagerkrafte nach Betrag und Phase (bezogen auf die rotorfeste Bezugs- 


ebene), so laBt sich die gesuchte, durch den Spaltenvektor fl beschriebene (approximierende) Un- 
wuchtverteilung am Rotor durch Auflésen des Gleichungssystems (11c) berechnen. 

b) Lésbarkeitsbedingung. Bei einem kontinuierlichen Gebilde geht die Anzahl der Frei- 
heitsgrade gegen Unendlich. Die Ordnungszahl der Matrizen ©, % und X in (11c) tibersteigt daher 
alle Grenzen. Erfahrungsgema8 ist der EinfluB der Kigenlésungen héherer Ordnungszahl von 
vernachlassigbarer GréBenordnung, so da8 man bei diesen Matrizen nur eine endliche Zahl von 
Eigenlésungen — ihre héchste Ordnungszahl sei mit o bezeichnet — zu beriicksichtigen braucht. 
In Abb. (2a) ist die Matrizengleichung (llc) in einem geometrischen Schema aufgezeichnet, in 


dem die an den Kastchen angeschriebenen Indizes die Anzahl der Spalten und Reihen der Matrizen 


kennzeichnen sollen. Man erkennt, daB sich bei der Ausrechnung des Matrizenproduktes © % X’ 
eine (J, n)-Matrix ergibt. Da n> I ist, bedeutet dies, daB zur Bestimmung der n Unbekannten 
U; zu wenig Gleichungen zur Verfiigung stehen. 

Die noch fehlenden Gleichungen lassen sich einfach dadurch beschaffen, da man mehrere, 
beispielsweise « MeBlaufe bei entsprechend zu wahlenden Winkelgeschwindigkeiten §2,, Q,... 2; 
vornimmt. Fiir die einzelnen MeBlaufe kann man je eine Gleichung gemaB (11c¢) 


i= 6%, X ti foil) 2a Sees) (11d) 
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anschreiben, die in ihrer Gesamtheit mit der in transponierter Darstellung geschriebenen Matrix 
Ke = (HH... WC) (11e) 

durch die Gleichung 
j*=2u, (11 f) 


die in Abb. (2b) schematisch dargestellt ist, ausgedriickt werden kénnen, deren Auflésung 
schlieBlich die gesuchte Unwuchtverteilung 


fi = Q-1 S* (11g) 


7 7u 
: 

tt Z tt 
7 7% 


7 7 
™ ia Tu 


Abb. 2. Schematische Darstellung der Matrizengleichungen (11d) und (11f). — a) fiir eine MeBlaufgruppe; b) fiir e-Meblaufgruppen 


liefert. 


ba) 


TL 


eo 


Das Gleichungssystem (11f) kann nur dann nach tt aufgelést werden, wenn die Produktmatrix 
& nichtsingular ist, wenn also 
el=n 


ist. Bei unganzzahligem Quotienten n/I ergeben sich iiberzahlige Gleichungen, die bei der Auf- 
lésung gestrichen werden miissen. Wurde beispielsweise n = 10 bei 1= 3 gewahlt, dann sind 
wegen € = 10/3 vier MeBlaufe vorzunehmen. Von den dann zur Verfiigung stehenden 12 Gleichun- 
gen diirfen aber nur 10 fiir die Weiterrechnung verwendet werden. In diesem Falle ware es ratsam, 
entweder die Zahl der MeBebenen auf n = 12 zu erhdhen (ergibt bei gleichem MeBaufwand ein 
genaueres Ergebnis), oder n = 9 zu wahlen, um mit ¢ = 9/3 = 3 MeBlaufen auszukommen. 

Es ist bei einem mehrfach gelagerten Rotor nicht unbedingt erforderlich, die Auflagerkrafte 
in allen Stiitzen zu messen, sofern es der konstruktive Aufbau des Rotors gestattet, die dann zu- 
sitzlich notwendig werdenden MeBlaufe bei hohen Drehzahlen vorzunehmen. Die hiermit ver- 
bundene Vereinfachung des meBtechnischen Aufwandes, es sei nur auf das betriebliche Auswuchten 
verwiesen, kann oftmals sehr erwiinscht sein. 


c) MeBlaufdrehzahlen. Die Elemente der Diagonalmatrizen %, sind nicht nur von der 
Ordnungszahl der Eigenlésungen, sondern auch noch von der Umdrehungszahl des Rotors sehr 
stark abhangig. Bei ungiinstig gewahlten Drehzahlen kann unter Umstanden der EinfluB wesent- 
licher Elemente von %, auf die Bildung der Produktmatrix 2 im Stérpegel der MeBanordnung 
untergehen, so daB8 dann im Verein mit den beim numerischen Rechnen unvermeidlich auftretenden 
Abrundungsfehlern die Auflésung des Gleichungssystems (11f) in Frage gestellt werden kénnte. 


mh 


~ dann zu erreichen, wenn die von den Elementen F gebildete Wertefolge bei jeder Drehza 
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Zur Wahrung der linearen Unabhangigkeit dieser Gleichungen miissen deshalb die MeBlaufdreh- 
zahlen so gewahlt werden, daB die Matrix & einen moglichst grofen Betrag annimmt. Dies ist 


: ‘ hl még- 
ichst stark verzerrt wird, ganz besonders aber dann, wenn bei den einzelnen MeBlaufen jeweils 


eine andere Eigenform besonders kraftig herausgehoben wird. 


Wieviele Eigenformen beim Aufbau der Matrizen ©, % und X zu beriicksichtigen sind, d. h., 
wie grofs die Ordnungszahl o gewahlt werden muB, hangt ab von den jeweils vorliegenden Ge- 
gebenheiten und ist von Fall zu Fall zu entscheiden. Die hierbei grundsatzlich zu beachtenden 
Beziehungen seien an Hand eines einfachen Gebildes, eines zweifach gelagerten homogenen Rotors, 
kurz erértert, fiir den in Abb. (3) die Drehzahlabhangigkeit der Matrizenelemente F’, in dimensions- 
loser Form dargestellt ist. Man erkennt, daB der EinfluB der Elemente hdherer Ordnungszahlen 


NIE 


ie 


A, 
89 7 2 3 a 


Abb, 3. Frequenzabhingigkeit der Elemente F, der Diagonalmatrix % fir einen homogenen zweifach gelagerten Biegungstriger. 
(M = Masse des Tragers; 4, = Eigenwert 1-ter Ordnung; a= Q = Erregerkreisfrequenz.) 


desto kleiner wird, je naher die Mefdrehzahlen an die kritischen Drehzahlen zu liegen kommen, 
bzw., daB umsomehr Eigenformen beriicksichtigt werden miissen, je weiter man von den kritischen 
Drehzahlen abbleiben will. Im allgemeinen wird man mit o = ¢ eine ausreichende Genauigkeit 
erlangen. Daraus folgt aber, daB man die Anzahl der MeBebenen in Einklang mit der Stiitzen- 
zahl I und der als gréBtzulassig erachteten MeBlaufdrehzahl wahlen muf. 

Selbstverstandlich wird man Messungen bei itberkritischen Drehzahlen nur mit vorgewuchtetem 
Rotor vornehmen kénnen. Sollte es sich als notwendig erweisen, innerhalb einer MeBreihe zusatz- 
liche Tarierunwuchten zur Erhéhung der Laufruhe anzuordnen, so braucht man den bereits ge- 
fahrenen Teil der Reihe keineswegs zu wiederholen, da der Einflu8 der Zusatzunwuchten rech- 
nerisch beriicksichtigt werden kann. Der hierfiir erforderliche Rechnungsgang gestaltet sich 
analog zu dem in Abschnitt 5 beschriebenen Verfahren, so daB hier auf Einzelheiten verzichtet 
werden kann. 

Zu erwahnen ware noch, daB bei den MeBlaufen keinerlei Auswuchtarbeiten im engeren Sinne 
zu leisten sind, so da8 man die vorgesehenen MeBdrehzahlen nur kurzzeitig wahrend des Ablesens 
der MeBwerte einzuhalten braucht und unter Umstanden ein ganzer Satz von Mefdrehzahlen 
in einem Zuge durchfahren werden kann. 


d) Genauigkeit des Verfahrens. Die Genauigkeit, mit der die durch (11g) festgelegte Un- 


-wuchtverteilung mit dem wirklich vorhandenen Unwuchtverlauf iibereinstimmt, hangt natur- 


gem4B ab von der Anzahl der beriicksichtigten Kigenlésungen, von der Art des angewendeten 
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Integrationsverfahrens und der Zahl der Integrationsintervalle. Gute Ergebnisse liefert die auf 
die Lagrangesche Interpolationsformel gegriindete Gaufsche Quadraturformel, bei der der Inte- 
grand durch eine ganze rationale Funktion n-ten Grades (n = Anzahl der MeBebenen) ersetzt 
wird!, Die Unterteilung kann dabei in Intervalle beliebiger Breite vorgenommen werden (besonders 
wichtig, wenn am Rotor Einzelmassen zu beriicksichtigen sind, oder wenn einige MeBebenen in 
die am Rotor konstruktiv vorgesehenen Ausgleichsebenen gelegt werden sollen), im Gegensatz 
zur Simpson-Formel, bei der eine aquidistante Intervallteilung gefordert wird, und die Annaherung 
an die wirkliche Funktion nur vom zweiten Grade ist. 

Die Anwendung der Quadraturformeln setzt einen stetigen Integranden voraus. Unstetig- 
keiten miissen bekanntlich durch Unterteilung des Integrationsbereiches in Teilabschnitte beriick- 
sichtigt werden. Wahrend man bei Laufern mit diskreten Massen den Sitz méglicher Unwuchten 
von vornherein kennt, lassen sich bei kontinuierlichen Gebilden keinerlei Voraussagen tiber den 
Verlauf der Unwuchtverteilung machen, so daB eventuell auftretende Unstetigkeiten bei der 
Berechnung nicht beriicksichtigt werden kénnen. Die hier grundsatzlich liegende Quelle méglicher 
Fehler kann leider nicht ausgeschaltet werden. 


e) Austarieren der Unwuchten. Mit (11g) ist die approximierende Verteilung der Un- 
wucht langs der Rotorachse gegeben. Durch Hinzufiigen von entsprechenden Zusatzmassen, bzw. 
durch Entfernen von Kérpermasse in den der Berechnung zugrundegelegten MefSebenen lassen 
sich die errechneten Unwuchten kompensieren. In diesem Falle sind die MeBebenen gleichzeitig 
Ausgleichsebenen. Will man den Ausgleich nicht in allen MeBebenen, sondern nur in einigen vor- 
gegebenen Ausgleichsebenen durchfiihren, so kénnen die hierfiir notwendigen Tarierunwuchten 
aus der Forderung hergeleitet werden, da im umlaufenden Rotor, und zwar bei allen Drehzahlen, 
die Eigenformen bis zu einer gewahlten Ordnungszahl o unterdriickt werden. Es muf also 

bk =0 (s=1,2,...,0) 


erfillt sein. 
FaBt man die Gesamtheit der in den Ausgleichsebenen anzuordnenden Tarierunwuchten U7; 
in einen — hier transponiert geschriebenen — Spaltenvektor 


ie = (U2 Us Caro Ug.) 
zusammen, so geht (8a) fiir den ausgeglichenen Rotor iiber in 
Re Ua © (+ Up). 


Die Auswuchtbedingung, das Verschwinden der ersten o-Eigenformen fiihrt dann zu den 0-Glei- 
chungen 

sy’ (I+ Up) = 0 (s = 1,2;,...,0), (12a) 
aus denen sich die Tarierunwuchten berechnen lassen. Ordnet man beispielsweise mit 0 = 3 die 


Ausgleichsebenen in der in Abb. 4a gekennzeichneten Lage an, so nimmt der Spaltenvektor tik 
wicder transponiert geschrieben, die Form 


fir = (00 Oy, 0 U7, 0 0 fiz, 0) (12) 
an, mit der sich aus (12a) ein inhomogenes Gleichungssystem 
*y(é) Ur, + *y(e5) On, + ty(e) Ur, = —¥ 
*y(zs) Up, + *y(e5) Un, + *y(e9) Or, = —*4' tL, leg 
*y(z5) Un, + *y(2s) Orr, + 8y(e6) Op, = —*y' B 
ergibt, das mit (11g) und den Setzungen 


ty (25) ®y(25) ®y(25) 


[ty (25) ®y(23) ®y(zg) 
| 
*y (Z) *y(%s) 2y (2s) 


1 Siehe FuBnote 2 : i i % oes . 
“Sedan ay New York suis 327, sowie Tables of eae Interpolation Coefficients, National Bureau of 
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als Matrix der reduzierten Eigenvektoren 


als verkiirzte Matrix der Kigenvektoren, 


Ur, 
als reduzierten Spaltenvektor der Tarierunwuchten, nach dem gesuchten reduzierten Vektor der 
Tarierunwuchten aufgelést werden kann: 


ft = — (€¥)7 &* f= — &¥)1 QF, (12d) 


Axl, 


Abb. 4. Dreifach gelagerter Rotor mit n = 9 MeBebenen. 


a) Tarierunwuchten in drei Ausgleichsebenen; 
b) Anordnung von Probe-Unwuchten in drei MeBebenen zur experimentellen Ermittlung der Lésungsmatrix Q-1 (siehe Gl. 16c). 


Aus diesen Gleichungen 1a8t sich die bemerkenswerte, von der in der Praxis tiblichen Auffassung 
abweichende Tatsache ablesen, daB ganz unabhangig von der Anzahl der Rotor-Lagerungsstellen 
jede zu beseitigende EKigenform nur je eine Ausgleichsebene erfordert, bzw., daB mit jeder Aus- 
gleichsebene eine Eigenform unterdriickt werden kann}. 


4, Rotor auf nachgiebigen Stiitzen (isotrope Lagerung). Ein mit der Winkelgeschwindigkeit 
a ya umlaufender elastischer Rotor mit isotroper Lagerung verformt sich unter der Wirkung 


- nichtausgeglichener Fliehkrafte und rotiert wie ein starrer Kérper. Es liegt auch hier ein quasi- 


stationdres Problem vor, das in der gleichen Weise wie unter Abschnitt 3 beschrieben, behandelt 
werden kénnte. Dabei miiBten aber die erforderlichen dynamischen Parameter auf den nachgiebig 
gelagerten Rotor bezogen werden, so daB die elastischen Parameter der Stiitzenkonstruktion 
implizit in den Eigenlésungen enthalten waren.Anderungen im Aufbau der Auswucht-An- 
ordnung wiirden eine neuerliche Berechnung der Eigenlésungen bedingen. Aus diesem Grunde 
ist der nachstehend beschriebene, zwar etwas umstandlich erscheinende Lisungsweg empfehlens- 
wert, der den Vorteil besitzt, daB die elastischen Kigenschaften der Lagerung nicht in die Eigen- 
lésungen eingehen. 

Die Gesamtheit der auf den umlaufenden Rotor einwirkenden Fliehkrafte bildet mit den Reak- 
tionen der Auflagerkrafte, den Stiitzkraften, ein Gleichgewichtssystem. Man kann daher ganz 
willkirlich zwei beliebige Kraftangriffspunkte als Auflager wahlen und den verformten Rotor als 
zweifach gelagerten Biegungstrager auffassen, der unter der verformenden Wirkung von zwei 
Kraftegruppen, den Fliehkraften und den von der wirklichen Lagerung herrithrenden Stiitzkraften, 
steht, wie dies beispielsweise in Abb. 5 fiir einen mehrfach gelagerten Rotor dargestellt ist. 


1 Optimierungsfragen mégen hier unberiicksichtigt bleiben, 
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Die Verformungen der elastischen Linie des Rotors lassen sich hierbei in zwei Anteile zerlegen 
und ergeben sich als Summe aus den Durchbiegungen des zweifach gelagert gedachten Ersatz- 
tragers und der durch die Senkung der gewahlten Stiitzpunkte hervorgerufenen Anteile. Fir eine 
ebene Belastung gilt fiir den /-fach gelagerten Rotor gemaB Abb. 5 


t eS 
@) = y¥@) +7 at Eu: 


Handelt es sich aan eine nichtebene Kraftegruppe, dann kann man in symbolisch-komplexer Dar- 


stellun 
i A L—2z x Zan : (13a) 
q2) =¥2+>Z-Ht7TH 


Abb. 5. Mehrfach gelagerter Rotor auf elastischen Stiitzen. 


a) Schemazeichnung des Rotors mit MeBebenen; 
b) verformter zweifach gelagert gedachter Ersatzrotor. 


schreiben. Bezeichnet man die von den Lagern auf den Rotor ausgeiibten Stiitzkrafte in komplexer 


Schreibweise mit P,, so ergibt sich mit den diesen Kraften zugeordneten Nadelfunktionen ¢,(z), 
die den gem4B (5) analogen Bedingungen entsprechen, die komplexe Belastungsfunktion zu 


PC) =A (4) + A) 02) + Sole) My Gls) + &)] + Ze) B 
bzw. mit (13a) : a 
© BG) =A (ele) + Foy(0) Mj) FH) + A (0) e) + BIG) Mya) 
é) J 


+4(F ti + F fi) (06) + Foe) M) + Dele F. (13b) 


Die in dem auch hier wahlbaren Lésungsansatz 
He) = & ke, "y(2) 


enthaltenen Entwicklungskoeffizienten ergeben sich, da (13b) formal mit (1d) iibereinstimmt, 
analog (2b) mit den Setzungen (6c) und (6d) zu 


b= aay || YA) & + Bye) a 


+ |) (o) + Z oy) M,) (* 


Sew a 1 ‘ 
= = Gy +F ai dz + 7% Ye) P;}. (13) 
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Bei der vorausgesetzten Isotropie der Lagerung des Rotors besteht Proportionalitat zwischen den 
Stiitzkraften und den Stiitzenverschiebungen. Fiihrt man die das dynamische Verhalten der 
Lagerkonstruktion kennzeichnenden komplexen Federkonstanten ¢; ein, so lassen sich die Stiitz- 
krafte durch die VerschiebungsgréBen 9, gemaB 

P,= — 4 (14a) 
ausdriicken. Man kann die auf der rechten Seite der Gleichung (13c) stehenden Integrale wieder 


mit Hilfe einer der bekannten Quadraturformeln auswerten und erhalt mit den durch die Parameter 
des Gebildes festgelegten Abkiirzungen 


L—z L—z2; 
8. = | Ss (2) o(2) dz + YS) == ry(z,) Mj, 
j 


, (14b 
b= | 77) ee) det SYP eyley M, 
J 
und der Setzung 
ry(z;) ¢; = b,; (i= 2,3,..., L—1) (14c) 
die Entwicklungskoeffizienten zu 
s 1 oe i Pm ne a 
esas Tet) (« U+ 8, 1 Pi By Sas | : by; i) 
woraus sich mit den Vektoren 
Gr = (g,,,0,0...08,1), 
6; = (0, b,,) 6, . . . by10) , 
sowie dem Zeilenvektor der komplexen Lager-VerformungsgréBen 
q a CA 92 Qs « : - qi) 
schlieBlich : 
y na 1 “A eS 
[ps —— bh, 14d 
= pagan [2 8+ (e—3 8) at 
-errechnet. Werden dann noch die Matrizen 
Osscalta et (14e) 
Ds) = (6, by. ° . 6) 


eingefiihrt, so 1aBt sich die Gesamtheit der Entwicklungskoeffizienten analog (8e) zu einem Spalten- 


vektor 


f=¢% je a + (6 —7 8) i 3 (148) 


zusammenziehen. : 

Die auf die beiden Stiitzen des zugrundegelegten Ersatzrotors wirkenden Auflagerkrafte be- 
stimmen sich wieder mit den den Eigenformen zugeordneten Auflagerkraften "S, bzw. 'S,, wenn 
man sich dabei auf die angendherte Wiedergabe durch o Kigenformen beschrankt, analog (lla) zu 


Yi ce ae ae ae flea) 


Mit (14f) und der Matrix der Kigenform-Auflagerkrafte des Ersatzrotors 
hag og 
e=(15, 00) 
sowie den zu (13d) analogen Gleichungen 
| A= ij, baw. A= Hi 
lassen sich die beiden Gleichungen (15a) in eine Matrizengleichung 


(2%) =6§ le t+ @ — 8) j 


e111 
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iiberfiihren, die nach einer kleinen Umformung und mit der weiteren Setzung 
ee €; OOS se 
0:0 Oe0'7; 


schlieBlich die Form il 
fe 1 ayy A a mat } | 

é+e3(, 8-6) q-o5zr ht (15) | 

annimmt. 
Die Mittelpunkte der ausgelenkten Lager des verformt umlaufenden Rotors beschreiben Kreis- ||) 


bahnen, die man sich von den als finlanteads Zeiger aufgefaBten VerformungsgréBen G qi erzeugt 


denken kann. Aus diesem Grunde erscheinen die Bewegungen der Lager in einer raumfesten jj) 


Achsialebene als harmonische Schwingungsbewegungen (die im allgemein gegeneinander phasen- |} 


verschoben sein werden), deren komplexe Amplituden nach Betrag und Phasenlage gemessen | | 


werden kénnen. Da der somit als bekannt zu betrachtende Vektor q von der Ordnung I (I< n) 
ist, liegt der gleiche Sachverhalt vor, wie er bereits in Abschnitt 3 beschrieben wurde: Es sind zu 
wenig Glechuseen vorhanden. Die noch fehlenden Gleichungen lassen sich in der erwahnten Weise 
Bévehattens fastlent man den Rotor mehreren MeBlaufen bei entsprechend auszuwahlenden Dreh- 
zahlen unterwirit. 

Die dem Verfahren zugrunde gelegte Isotropie der Lagerung wird in Wirklichkeit nicht erfullt |} 
sein. Bei den iblichen technischen Anwendungen ist aber der dadurch bedingte Fehler meist 
verhaltnismaBig klein, so daB die Ergebnisse nicht nennenswert beeinfluBt werden diirften. 

AbschlieBend sei nochmals ganz besonders darauf verwiesen, daB sich die in den Gleichungen | 
dieses Abschnittes vorkommenden Eigenlésungen auf den zweifach gelagert gedachten Ersatzrotor 
beziehen. Obwohl in den meisten Fallen damit zu rechnen sein wird, daf die auszuwuchtenden 
Objekte ihrer Wichtigkeit entsprechend schon von der Konstruktion her einer schwingungstech- 
nischen Berechnung (neuerdings meist mittels Rechenautomaten bis zu hohen Ordnungszahlen) | 
unterzogen wurden, so ist dies fiir das vorstehend beschriebene Verfahren bedeutungslos, da diese i) 
Rechnung selbstverstandlich fir die wirklich vorhandene Lagerungsart und nicht fiir die hier 
ersatzweise gewahlte zweifache Auflagerung durchgefithrt worden sein wird. In der hier notwendig || 
werdenden nochmaligen Durchrechnung liegt eine wesentliche Erschwernis des Verfahrens, so daB 
man es nur dann anwenden wird, wenn iiber das elastische Verhalten der Lagerkonstruktion von 
vornherein keine verbindlichen Aussagen gemacht werden kénnen. 


5. Experimentelle Lésung. Der nicht unbetrachtliche Aufwand an Berechnungsarbeit zur Be- 
stimmung der Matrix % [(11f) bzw. Abb. 2b] kann durch eine meBtechnische Analyse des aus- 
zuwuchtenden Gebildes umgangen werden, wenn man den dazu erforderlichen, ziemlich umfang- 
reichen mefitechnischen Aufwand in Kauf zu nehmen bereit ist. Voraussetzung ist hierbei, daB |} 
es der konstruktive Aufbau des MeBobjektes gestattet, die erforderliche Anzahl von Tarierunwuch- |} 
ten in den gewahlten MeBebenen anzubringen. 

Das Verfahren eignet sich sowohl fiir starre als auch fiir nachgiebige Lagerung und sei hier fiir 
den Fall starrer Stiitzen kurz angedeutet. Mit dem als I-fach gelagert vorausgesetzten Rotor 
unternimmt man e-MeBlaufe bei den gema8 Abschnitt 3c) auszuwahlenden Drehzahlen. Dabei_ 
mift man die /-Auflagerkrafte nach Betrag und Phase (bezogen auf eine rotorfeste Bezugsebene). 


Damit ist der Spaltenvektor I* gem48 (11b) und (11 e) bestimmt. Nun bringt man Tarier- \ 


unwuchten in einigen der vorher festgelegten n-MeBebenen an, deren Spaltenvektor mit A, fi be- | 
zeichnet sei, der fiir die in Abb. 4b beispielsweise gewahlte Anordnung mit 


uo A 


4, Us = Ty my e% , a My = 14m, e'*, 4, Uy = Tm, e™ 
und it = y=] als imaginare Einheit in transponierter Form durch 
a= (0 4,U, 0 PROP 00 4,U, 0 0) 
gegeben ist. Der mit diesen Tarierunwuchten besetzte Rotor wird nun wieder ¢-MeBlaufen bei | 
den gleichen Drehzahlen wie vorher unterworfen. Die komplexen Amplituden der dabei gemessenen_ | 
Auflagerkrafte faBt man zu einem neuen Spaltenvektor a* gemaB (lle) zusammen. Da fiir beide | 
MeBlaufgruppen Gl. (11g) giiltig ist, 14Bt sich ‘| 
fi=2- 9, (16a) | 


“~ “A am 


+ 4U= 2-1 f* (16b) 
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schreiben. Subtrahiert man diese Gleichungen, so ergibt sich, wenn 


A, Q* = td sre: Ms 
gesetzt wird, 


A pHs tac. A * . 


Der beschriebene Vorgang wird insgesamt n-mal mit abgeanderten Unwuchten A, U;, die jeweils 
so anzuordnen sind, daf sich die den einzelnen MeBlaufgruppen zugeordneten Vektoren ,9{* in 
ihren Elementen méglichst stark voneinander unterscheiden (zur Gewahrleistung der linearen 
Unabhangigkeit eines zu (11g) analog aufzubauenden Gleichungssystems), wiederholt. Die so 
erlangten n-Gleichungen 


“A 


alt = ee (a eg) 


kénnen mit den beiden (n, n)-Matrizen 


A = (4h 45"... 4,50) 
zu emer einzigen Matrizengleichung zusammengezogen werden, aus der sich die gesuchte Matrix 
eee Zu 
“A Aaa aay 
ee ea) (421*) (16c) 


e 


bestimmen 148t, und womit dann die gesuchte Unwuchtverteilung mit Hilfe der Gleichung (16a) 
berechnet werden kann. Man wird, um den Aufwand an Berechnungsarbeit klein zu halten, bedacht 
sein miissen, mit einer méglichst kleinen Zahl von Tarierunwuchten — im Extremfall eine je MeB- 
laufgruppe — auszukommen. 


Bei diesem Verfahren erlangt man durch systematisches ,,Befragen“ Einblick in den struk- 
turellen Aufbau des MeBobjektes, wie dies in grundsatzlich gleicher Weise beim Auswuchten starrer 
Rotoren iblicherweise durchgefiihrt wird. Dieser Vorgang ist, wie man sich leicht ttberzeugen 
‘kann, als Sonderfall in der allgemeinen Lésung gemaé®B (11d) vorhanden, wenn man beachtet, 
daB ein in einer Ebene frei-freies stabférmiges Gebilde den zweifach zu zahlenden Eigenwert null 
‘mit den Higenfunktionen y(z) = 1 baw. y(z) = (z— %)/%, mit z) als Schwerpunktsabszisse des 
Gebildes, besitzt. 


6. Dimensionsbetrachtung. Bei der numerischen Behandlung der Lésungsgleichungen geht man 
zweckmaBigerweise auf dimensionslose Beziehungen iiber. Setzt man die HKigenfunktionen "y(z) 
wie allgemein tiblich dimensionslos an, so erhalten die Elemente von ’© als die den Eigenformen 
-zugeordneten Auflagerkrafte "S, die Dimension einer Kraft je Langeneinheit, wahrend die Elemente 
‘von "¢ und X dimensionsfrei bleiben. 


Mit den fiir die Krafte, Massen und Langen passend zu wahlenden Bezugsgré8en Py bzw. Mo 
baw. |, und den Setzungen 


z= (1/l)z, j = (1/M,) M;, 
W; = (1/b) Wi. e; = (1/h) @;, 
q 0(2) = (%o/Mo) of) » u; = W, Q(z) e(@:) + M; & = (1/My hy) ti; » 
e(z) = (1/l) el) » i, = f 'y*(2) Q(2) dz = (1/Mo) o, » 


Y= (1/P,) X, % = My hs 


6 = (1,/P.) ©» N= (1/M, 4h) U 


(17) 


W= SFLU 


338 E. Hiibner: Auswuchten elastischer Rotoren, ein Problem der Strukturanalyse Ingenieur-Archiv || 


iiberfiihren. Damit wird aber auch die in (11f) vorkommende Produktmatrix, sie werde mit g 
bezeichnet, dimensionsfrei, so daB sich mit der Festlegung 


Oi = (1) Pyare 
die endgiiltige Lésung 
= (Lk) W* 
ergibt. Auf die gleiche Weise gelangt man mit 
it = (1fM, 1,) Ut 
zu 
it = — @F) 1 Ee" (2) 


[siehe (12d)]. SchlieBlich erhalt man mit den weiteren Setzungen 


q= (1h) G- C = (h,/P,)C, 
@ = (1/M,) G, A= (Myl/P,)A. 
B= (lb/P,) B, 


und den Gleichungen (17) die letzte Lésungsgleichung dimensionsfrei: 


E+ S5(78—G)\q=S5eU. 


(Eingegangen am 12. September 1960.) 
Anschrift des Verfassers: Privatdozent Dr. Erhard Hiibner, Aachen, Im Brockenfeld 17. 
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Berechnung der Strémung des instationiren vielstufigen Plattengitters 


Von Ernst-Giinther Feindt 


1. Einleitung. In vielstufigen Axialverdichtern sowie in vielstufigen Gas- und Dampfturbinen 
sind sehr viele Schaufelreihen unmittelbar hintereinander angeordnet, wobei oftmals samtliche 
Schaufeln, Leit- und Laufschaufeln, dasselbe Profil besitzen. Fiir den Konstrukteur ist es wichtig 
zu wissen, wie sich die Geschwindigkeitsverteilung der Schaufel einer mittleren Stufe unterscheidet 
von der Geschwindigkeitsverteilung, die dieselbe Schaufel als Einzelprofil oder als Schaufel eines 
Einzelgitters besitzt. Im folgenden wird dieser Vergleich durchgefiihrt am ungestaffelten Platten- 
gitter, weil dieses mit einfachen mathematischen Mitteln exakt berechnet werden kann. 

In der Praxis steht das Leitrad still, und das Laufrad bewegt sich mit einer Geschwindigkeit 
2 U in Umfangsrichtung. Das strémende Medium rotiert dann im Mittel in der ganzen Maschine 
mit einer Umfangsgeschwindigkeit U, so daf Laufrad und Leitrad eine entgegengesetzt gleiche 
Umfangsgeschwindigkeit U bzw. — U relativ zum Medium haben. Fiir die Rechnung ist es an- 
schaulicher von der Vorstellung auszugehen, da das strémende Medium im Mittel keine Umfangs- 
geschwindigkeit besitzt, und sich dafiir Leit- und Laufrad mit der Umfangsgeschwindigkeit — U 
bzw. + U bewegen. Die Durchstrémung des vielstufigen (genauer: unendlich vielstufigen) Platten- 
gitters wird im folgenden berechnet als Uberlagerung einer Verdrangungsstrémung mit einer Zir- 
kulationsstrémung und einer Parallelstrémung parallel zur €-Achse. Die Parallelstrémung wird 
allerdings im folgenden als trivial nicht weiter beachtet. Die Strémung wird als inkompressibel 
und wirbelfrei vorausgesetzt. 


2. Vorbemerkungen zur konformen Abbildung. Die Berechnung der Strémung des vielstufigen 
Plattengitters geschieht in der Weise, daB ein Periodenrechteck der ¢-Ebene (Abb. 1, 2 und 5) 
auf eine volle x-Ebene (Abb. 3 und 6) konform abgebildet wird. In der x%-Ebene erscheint nur je 
eine Platte des Gitters a und des Gitters b. Die x-Ebene wird dann weiter abgebildet auf ein Recht- 
eck der z-Ebene (Abb. 4 und 7). In der z-Ebene kann die Strémung ohne Schwierigkeit ermittelt 
werden. Die ganze Rechnung muB fiir die Verdrangungsstrémung (Index q) und die Zirkulations- 
strémung (Index J’) getrennt durchgefiihrt werden. Dies liegt daran, daB sich dann in der 


Gitterb Giffer a Gitfer b Gitter a 


i i i i i i i i lattenreihen gezeichnet. 
Abb. 1. Das unendlich vielstufige Plattengitter. Es sind nur vier der unendlich vielen Pla’ 
Fir das Jaklcabéispiel ist = Lt = 0,70711, s = 0,13435. Gitter a: U =—1 (bzw. 0); Gitter b: U=1 (bzw. 2). 


24 
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Gitter a 


1, 


Abb. 2. Skizze der Verdrangungsstrémung. LMTR ist das Periodenrechteck, das konform abgebildet wird auf die volle %q-Ebene. 


| Hy = ST (520506, Oy =73,812°) 


B 
lL. VEN fas M 
Pa Mi (eS EE SEN Lon Schlitz, Ferg 


RETIN TNE 


" 


Abb. 3. Zwischenabbildungen der Verdrangungsstrémung Abb.2. In die x _-Ebene sind die Stromlinien einskizziert. DieStromlinien « 
bis | entsprechen den gleichbezeichneten Stromlinien von Abb. 2. “xT =—*R= Ijkg = 1,04129, %9=—*N = 1, xp = — #4 = 0,998 20, 
“Cc = — xp = 0,21218. 
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Abb. 4. Bildebene der Verdrangungsstrémung. Die eingezeichnete Verteilung der Normal- 
komponente w,, der Geschwindigkeit ist positiv, wenn sie ins Innere des Rechtecks weist. 


d 
Giffer a Gitter b 


ce 


G oe Se % H 


l+s sal 


Abb. 5. Skizze der Zirkulationsstrémung. GHHG ist das Periodenrechteck, das konform abgebildet wird auf die volle xp-Ebene. 
24% 
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ganzen x,- bzw. x,-Ebene eine regulare Strémung ergibt, wenn die Ecken G und H bzw. L, M, L 
und R des abzubildenden Periodenrechtecks in die Staupunkte der Strémung gelegt werden. Die 
Abbildung ist in den Eckpunkten des Periodenrechtecks namlich nicht konform. Dieser Umstand 
stért nicht, wenn die Eckpunkte des Periodenrechtecks in die Staupunkte der Strémung gelegt 
werden. Man sieht dies leicht ein, wenn man die Abbildung genauer betrachtet. Fiir die Ab- 
bildung der ¢-Ebene Abb. 2 auf die x,-Ebene Abb. 3 gilt beispielsweise folgendes: 

Die obere Halfte ONLPMQO des Periodenrechtecks geht in die obere Halfte Im {x,} > 0 
der x,Ebene iiber. Entsprechend geht die untere Halfte ONRSTQO des Periodenrechtecks 
in die untere Halfte Im {rq} < 0 der %- Ebene iiber. In Abb. 2 und Abb. 3 sind einander ent- 
sprechende Punkte mit gleichen Buchstaben bezeichnet. Die Punkte P und S der¢-Ebene gehen in den 
unendlich fernen Punkt %, = 00 der x%4-Ebene iiber. Die x,-Ebene muf man sich langs der negativen 
x-Achse von N bis —oo und langs der positiven Achse von Q bis + 00 aufgeschlitzt denken. Auf 
den oberen Schlitzrandern liegen die Punkte L und M, auf den unteren Schlitzrandern liegen die 
Punkte R und T. Der Streckenzug ONLP der ¢-Ebene geht in die Strecke ON sowie den oberen 
Schlitzrand NL —oo der x,-Ebene iiber. Entsprechend geht der Streckenzug ONRS der 
C-Ebene iiber in die Strecke O N und den unteren Schlitzrand N R—oo der x,-Ebene. Die rechten 
Winkel x ONL, x ONR, x OQM und x OQT sowie die Innenwinkel ~< NLP, <x NRS, 
X QTS und x QMP des Periodenrechtecks werden bei der Abbildung gestreckt von 7/2 auf z, 
und die Winkel 7/4, unter denen die Staupunktstromlinien bei L, N, R, M, Q und T auf die ver- 
tikalen Rechteckseiten auftreffen, werden auf 7/2 aufgebogen. Aus diesen Uberlegungen folgt, 
daB die innerhalb des Periodenrechtecks der ¢-Ebene verlaufenden Stromlinien bei der Abbildung 
auf die x,-Ebene so deformiert werden, da8 sie sich gema8B der Skizze Abb. 3 langs der Schlitz- 
rander ohne Knick schlieBen. 

Entsprechende Betrachtungen kann man anstellen fiir die konforme Abbildung der Zirku- 
lationsstrémung Abb. 5 auf die xp-Ebene Abb. 6. Ein neuer Gesichtspunkt ergibt sich bei dieser 
Abbildung dadurch, daB nur jeweils eine Halfte einer Schaufel des Gitters a und des Gitters } 
konform abgebildet wird, die Schaufeln also bei der Abbildung zerschnitten werden. In der Um- 
gebung der Schnittstelle Punkt E verhalt sich die Abbildungsfunktion! wie die Funktion x, + 1 
= (¢—¢x)*- konst. Entsprechendes gilt fiir die Umgebung von Punkt F’. Als Endergebnis dieser 
Uberlegungen findet man, daB durch die konforme Abbildung die Zirkulationsstrémung des viel- 
stufigen Gitters iibergeht in die Zirkulationsstrémung um zwei einzelne Platten. 


3. Konforme Abbildung fiir die Berechnung der Zirkulationsstrémung. Fiir das Zahlenbcispiel 
wurden die in der Bildunterschrift von Abb. 1 angegebenen Werte gewahlt. Die konforme Ab- 
bildung des Periodenrechtecks GHHG (Abb. 5) auf die xp-Ebene erfolgt* mit Hilfe des ellipti- 
schen Normalintegrals erster Gattung bzw. mit dessen Umkehrfunktion, der Jacobischen 
sn-Funktion. Dem Rechteckverhiltnis (I + s)/t = 0,84146 entspricht** das Modulquadrat 
k7 = 0,31805 sowie die Halbseiten K; = 1,72461, Kp = 2,04955 des Normalrechtecks. Der Zu- 
sammenhang zwischen der ¢- und der xp-Ebene ist also gegeben durch die Gleichung 


“~ °r 


= ERs dx, Sr do, 
2 Kr ya— x3.) a —k?, xh) 2 Kr yl — sin? ce sin? Pr 
0 0 


Fs ; Fp» 8, = 34,3329) 
ee KrF Gr Or = are sin kp) = S fe (1) 
bzw. durch deren Umkehrung 
fm 2Kr os 5 
%*%r = sn ee ig Op ete . (2) 


Diese Beziehungen sind in Tabelle 1 ausgewertet worden durch quadratische Interpolation in den 
Tafeln von Hayashi. Fir die Plattenendpunkte ergeben sich die Werte 0,27083, — 0,27083, 1 und 
—1 (Abb. 6). Im nachsten Schritt der Berechnung wird die ~,-Ebene mit Hilfe der Gleichung 
(rc = 0,27083) 

1 


Tr = pate (3) 


1 A. Betz, Konforme Abbildung, Ziff. 42. Berlin 1948. 
* A. Betz, a. a. O., Ziff. 63. 


** Vel. A. Betz, a. a. O., Abb. 203, oder K. Hayashi, Fiinfstellige Funktionentafeln, Tafel XXX, Berlin 1930. 


XXX. Band 1961 E.-G, Feindt: Berechnung der Stromung des instationaren vielstufigen Plattengitters . 343 


so auf eine tr-Ebene abgebildet, daB die inneren Plattenendpunkte B und C in tp — -+-1 fallen. 
_ Fir die auBeren Plattenpunkte F und E ergibt sich tp = + 3,69237. Die anschlieBende Abbildung 


auf das Rechteck der zp-Ebene (Abb. 7) erfolgt wieder mit Hilfe eines elliptischen Normalintegrals 
erster Gattung* 


ty ; 
= dt. dp : : 
2 | va — t}) (1 — kip th) = —— = F(B, Orr = 74,293°) . (4) 
; 0 


Jn 2 
pee | pom (e800, Y1=34350°) 


~0,07083 G27083 R, tp 


unten P 


wf ERT 
m Up Tr = games 


569237. Aotp 
F° 


369237 Qohen 


Onben 


Cunten 


Abb. 6. Zwischenabbildungen der Zirkulationsstrémung Abb. 5. 


Der Modul kz, der Komodul kp sowie die GréBe 3; werden berechnet aus den Gleichungen 


1 I 


kjrp= 


krr = V1 — kip =0,96263 , 
07 r = are sin ky p = 74,293° . 


GemaB Betz** besteht zwischen sin f (von Betz mit t, bezeichnet) und t, die Beziehung 


— 
ge |r wpa : (5) 


Durch (4) und (5) ist die Zuordnung der Schaufelkonturpunkte der t,- und zp-Ebenen hergestellt. 
Gleichung (4) und (5) sind in den ersten sechs Spalten von Tabelle 2 ausgewertet. Die halbe Breite 
K,r und die halbe Héhe Ky; des Rechtecks der zp-Ebene sind gegeben durch die vollstandigen 
Normalintegrale Krr (krr) = 1,60085 und K7r (krr) = 2,72482. Der AbbildungsmaBstab zwischen 
der ¢- und der zp-Ebene kann nach der Kettenregel 


Ses Sea 
dy, dz, dt, dy, 


(6) 


leicht aus (1), (3) und (4) ermittelt werden. Es ergibt sich die Formel 


\i— (l—hi pt?) 


(l — x3) (lL— kx) |’ 


die in der vorletzten Spalte von Tabelle 2 ausgewertet ist. 


ag 
dy. 


I+ os 
Be oo K 
r 


(7) 


* Vgl. A. Betz, a. a. O., S. 232, Absatz 2 und Gleichung (63, 37). 
** A_a. O., Gleichung (63, 38). 
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4. Berechnung der Zirkulationsstrémung. Die Strémung der tr-Ebene ist eine Zirkulations- 
strémung um zwei einzelne Platten. Die Zirkulationen der beiden Platten sind entgegengesetzt 
gleich gro. Gema® Betz* entspricht dieser Stromung in der z--Ebene die in Abb. 7 dargestellte 
vertikale Parallelstr6mung w=ic. Die gesuchte Geschwindigkeitsverteilung der Platten der 


¢-Ebene ist also 


dy 
wr oe (8) 


Beriicksichtigt man, da8 nur die halben Schaufeln der C-Ebene auf die z--Ebene abgebildet wurden, 
so erhalt man fiir die Schaufelzirkulation J” die Gleichung 


Yj=—l,=4Kire. (9) 
Gleichung (8) ist mit dem in Abschnitt 7 berechneten Wert c = 0,11107 in der letzten Spalte der 
Tabelle 2 ausgewertet. 


2 Tabelle 1. Berechnung der Zirkulations- 
is r strémung. Auswertung der Gleichung (1) 
: F aa | 4 
1 Gl. (1) 
F(pp. Op = 34,332°) 
Aopen Onben Ry sega Sie sa iin 
Sa) willkirlich gewahlt nach Hayashi 
Tab. XXVIL 
—— 
0° 0 

B GP Toes 10° 0,042 65 
20% 0,085 70 
30° 0,12956 
unten i aS 40° 0,174.58 
S 50° 0,22103 
c iy ae 0,244.83 
a: 60° 0,269 05 
65° 0,293 65 
F 70° 0,318.60 
B T5392 0,343 87 
80° 0,369 38 
3 0,379 65 
Rane re 84° 0,389 93 
Abb, 7. Bildebene der Zirkulationsstrémung. Die Zirkulations- 86° 0,400 23 
strémung von Abb. 5 erscheint in der zy-Ebene als Parallel- 88° 0.41055 
strémung mit dem Geschwindigkeitsbetrag |c]. 90° 0,420 87 


5. Konforme Abbildung fiir die Berechnung der Verdrangungsstrémung. Die konforme Abbildung 
der Verdréngungsstrémung unterscheidet sich von der konformen Abbildung der Zirkulations- 
strémung nur dadurch, daf} ein etwas anderes Periodenrechteck der ¢-Ebene abgebildet wird, und 
zwar ist das Rechteckverhaltnis gemafs Abb. 2 gleich 2 (I + s)/t = 1,68292. Nach Hayashis Tabelle 
XXX gehért hierzu ky = 0,96035, On = are sin k, = 73,812°, K, = 2,69727, Ki, = 1,60273 
Damit erhalt man als erste Abbildungsfunktion : nthe 


wn 


ras ee ives . 
C= i ——t _- = F@,, 0, = 13,612- 
K - Op 0 = 13,082 10 
te) CC te oe 
oder als Umkehrfunktion geschrieben 
ss 
tq = sn tee ae 73,812") (11) 
Fiir die anschlieBenden Abbildungen der x,-Ebene auf die Tq-Ebene ergibt sich (%,¢ = 0,21218) 
1 
Vq = €e %q > (12) 


so A, a. OFS. 245, 
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Tabelle 2. Berechnung der Zirkulationsstrémung 


B tp %p or ¢ 1 |\dé/dy p| 3 he 
Gl. (5) Gl. (3) Gl. (4) Gl. (7) Gl. (8) 
xurlich | . 0,27083 : interpolior’ ix] 7 Ae» 8] p = 14,293") sae ee aes 
gewahlt i —k'Prsin?p |” Od eS tated Tabelle I ies poets ; i: a f) i —'i rtp) dé/dy p 
1—x},) (1—kp xp) 

02 1 0,27083 | 15,713° 0,067 2 0 0 oo 
10° 1,014.26 0,27469 | 15,944° 0,0676 0,175 35 0,01133 9,803 2 
20° 1,059 04 0,28682 | 16,668° 0,0708 0,355 81 0,023 34 4,758 8 
30° 1,140 84 0,30897 | 17,997° 0,076 8 0,948 59 0,036 84 3,0149 
40° 1,272.94 0,34475 | 20,167° 0,0860 0,756 90 0,053.05 2.0937 
50° 1,480 52 0,40097 | 23,639° 0,1012 0,995 77 0,074.05 1,4999 
29° 1,62605 0,44038 | 26,128° 0,112 4 1,13104 0,087 45 1,2701 
60° 1,810 72 0,49040 | 29,367° 0,1268 1,28063 0,10378 1,0702 
65° 2,046 11 0,95415 | 33,653° 0,1460 1,431 30 0,124.16 0,8946 
A08 2,345 64 0,63527 | 39,441° 0,1723 1,639 66 0,150.15 0,7397 
Oe 2,799 68 0,758 24 | 49,310° 0,2180 1,908 60 0,19113 0,981 1 
The 2,969 57 0,80425 | 53,540° 0,238 2 2,009 58 0,206 83 0,5370 
80° 3,14208 0,85195 | 58,428° 0,263 8 2,11670 0,224.32 0,495 1 
82° 3,309 04 0,89619 | 63,669° 0,287 2 2,22995 0,24158 0,459 8 
84° 3,461 06 0,93736 | 69,624° 0,3170 2,349 03 0,257 90 0,430 7 
86° 3,983 56 0,97054 | 76,059° 0,348 8 2,473 14 0,271 67 0,4088 
88° 3,664 20 0,99238 | 82,925° 0,384.4 2,60101 0,280.17 0,396 4 
90° 3,692 32 1 90° 0,420 7 2,724.82 0,284.37 0,3906 


und die Abbildung der x,-Ebene auf die z,-Ebene folgt aus den Gleichungen (f ist Parameter) 


Tq 
d, 
r= { aes eas 58 SIN erg (ly (oe: ) 
7 ya 7) (I Ki, A) | Aes 


I 
Vs i — Ky, sin? B 


Wie. in Abschnitt 3 erhalt man den Modul k;, aus der Lage des Plattenendpunktes 1, p. Es ist 
ky, = 1/4,70450 = 0,21256. Damit ergibt sich ky, = 0,97115, 07, = are emi bp = 715732", 
K,, = 1,58900, Kj, = 2,95724. Der AbbildungsmaBstab d&/dy, zwischen der ¢- und der z,-Ebene 
wird wie in Abschnitt 3 berechnet: 


dé 1+ s|]/A—%) (44%) 
—— A ° 14 
dy, qG K, i Gs: ae i) 


Die Abbildungsgleichungen (10) bis (14) werden zweckmaBig nach demselben Schema ausgewertet 
(Tabelle 1 und 2) wie die entsprechenden Gleichungen der Zirkulationsstr6mung. 


6. Berechnung der Verdraingungsstrémung. In der ¢-Ebene sind die Geschwindigkeitskompo- 
nenten normal zu den Schaufeloberflachen gleich U. In der %_- Ebene ist daher die schaufelnormale 


Geschwindigkeitskomponente w,,(2,) 
Wn(Zq) =U. (15) 


Um diese Grenzbedingung in der z,-Ebene zu erfiillen, werden die Rechteckseiten AA und DD mit 
Quellbelegungen q = 2 w,(z,) versehen (Abb. 4) und das Rechteck ADDA durch Spiegeln an den 
vertikalen Rechteckseiten 4A und DD sowie periodisches Fortsetzen iiber die Seiten AD auf die 
volle z,-Ebene erganzt'. Die Tangentialgeschwindigkeit, die von diesen Quellbelegungen induziert 
wird, ergibt sich aus der Gleichung 


Ky i(yg—Y, i(¥q—No 
4 a a 9, (s ) 


w,(Z,) = a W,(2q) g ce me A (ies =) Yq ( ) 
Pah eee ere » || Sa 
eR ie *\ 4 Krq 4 Krq 


1 F.-G. Feindt, Ing.-Arch. 30 (1961), S. 88—95. 
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mit w,(z,) nach (14) und (15). Die gesuchte Tangentialgeschwindigkeit in der (-Ebene ist dann 


w,(0) = 2(z9) SE (17) 


Die 3-Funktionen werden zweckmaBig aus den ersten beiden Gliedern ihrer Reihenentwicklung 
berechnet. Durch Zerlegen in Real- und Imaginarteil ergibt sich* 


ps a seaa i 3 2(ya— ‘) 
“Yq z) . (Yq — Yq) (Ya—Yq 
fie) (y bt SY Do Pl SSS ee et h2 f ea ees 
»,( ‘Rig pa Shag, el ae 
a are aa 4Krq _. ™¥q—Y. _ 32% 9)’ 
a, u(¥q —Yq) 1 Gin die NPE : ) yg Sin yaa rE z ’ 
oe aa ae (18) 
i(¥q— 4) (Sq a) __ 3(¥q — Yq) 
oe hh ee) ae Ve ee 
o:( ee ain Sin 4 Kg + 3h? Gin 4Krq 
(Yq — Yq) a Kye ——— Ya se ee oe 
Ee Gof Sy 
4 Krqg 4 Ktq 4 Krq 


mit h, =e” “1a/C*14) — 0,05375 . 
Der Integrand von (16) ist ersichtlich rein imaginar. Das Integral (16) wird numerisch aus- 


gewertet. 


7. Die Joukowsky-Bedingung. Damit die Plattenhinterkanten nicht umstrémt werden, mu die 
Zirkulationsstrémung zusammen mit der Verdrangungsstrémung an den Plattenhinterkanten die 
Geschwindigkeit Null ergeben, d. h. es mu sein 


(19) 


worin alle GréSen in den Punkten C der z,- bzw. z,-Ebene zu nehmen sind. (Aus Symmetriegriinden 
kénnen die Punkte C statt der Plattenhinterkantenpunkte B genommen werden.) Aus (19) folgt 


a 
dy 
és sm.) SF (20) 
dyy 
Fiir w,(z,) findet man durch Einsetzen von y, = 0 in (16) den Wert w,(z,) = — 0,11252. Der 


zweite l’aktor auf der rechten Seite von (20) hat in C die unbestimmte Form 0/0. Nach (7) und (14) 
kann man fiir den Grenzwert schreiben 


Lae Vz Kg (1 — 2) (1 —k? x8) (1 — 7?) (1 — 43,73) 
GE et P %qo2 Kp (EA) (Lhe 2) (L —x2) (1 — be x8)” 
dyq Trl 
wenn hierin die Werte von K,, K,, ky krq kr, krr sowie %q = %qc = 0,21218 und xp = xo 
= 0,27083 eingesetzt werden. Es ergibt sich 


dé 

iy, _ 4 /t=# 
“= 0198891 lim j= ; 
dyg Feel 


woraus weiter durch Ableiten des Zahlers und des Nenners und Einsetzen von (3) und (12) folgt 


ge dt, dx, 
dy “dt “dl 

hae on 0,98891 O “dig > 0,87530 : “dig ° 
dyq de dt 


* H.-G. Feindt, a.a.O. und A. Betz, a. a. O Gleichungen (78,9) und (78,10). 


bere Maas 
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Die in dieser Gleichung auftretenden Differentialquotienten kénnen aus (2) und (11) berechnet 


werden, wenn man in diesen Gleichungen die sn-Funktionen durch die ersten beiden Glieder der 
Reihenentwicklung ersetzt gemaB der Formeln 


1+ k2 
xp = 4,09909  ——* (4,09909 ¢)8 , 
1+ k2 
Hq = 3,20546 6 — . 1 (3,20546 ¢)° . 
Das Endergebnis der Rechnung lautet 
a 
dy, 
2 = 0,98702 . (21) 
dyq 


Mit dem angegebenen Wert von w,(z,) und (21) folgt aus (20) c = 0,11107. 


8. Rechenergebnis. In Abb. 8 ist fiir die Schaufeloberseite des Gitters b die Geschwindigkeits- 
verteilung der Verdrangungsstrémung aufgetragen. Die Geschwindigkeitsverteilungen der Schaufel- 
unterseite des Gitters b sowie die Geschwindigkeitsverteilung des Gitters a unterscheiden sich hiervon 
gema8 Abb. 2 nur durch das Vorzeichen. Zum Vergleich ist in Abb. 8 mit eingetragen die Geschwin- 
digkeitsverteilung des Einzelgitters, das die gleiche Teilung hat wie das vielstufige Gitter, sowie die 


G8 


OY 


- 04 


=48 


indigkei i a 6 berseiten), U = 1. 
Abb. 8. Geschwindigkeitsverteilung der Verdrangungsstromung (Schaufelo iten), 
a) Gin des S ciseafigen ee und Einzelgitter gleicher Teilung; b) Einzelplatte. \ 
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Abb. 9. Geschwindigkeitsverteilung der Zirkulationsstrémung (Schaufeloberseiten). 
a) Gitter b des vielstufigen Gitters; b) Einzelgitter gleicher Teilung; c) Einzelplatte. 
Die Zirkulation ist in allen drei Fallen gleich, 7” = 1,21060. 


Abb. 10. Geschwindigkeitsverteilung fiir die Strémung mit glattem AbfluB an den Schaufelhinterkanten, U=1. 
a) Gitter a des vielstufigen Gitters; b) Citter b des vielstufigen Gitters; c) Einzelgitter gleicher Teilung. 
Zirkulationsverhiltnis des vielstufigen Plattengitters T/L, = 0,545, des Einzelgitters 7/I, = 0,724. 
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Geschwindigkeitsverteilung der Einzelplatte. Die Vergleichsstrémungen mu8 man sich so ent- 
standen denken, daB sich eine Einzelplatte bzw. ein Kinzelgitter in ruhendem Medium senkrecht 


zu der Platte bzw. parallel zur Gitterfront bewegt mit der Geschwindigkeit U = 1. Die Ver- 


draéngungsstrémung des vielstufigen Plattengitters entsteht dadurch, daB in einem ruhenden 
Medium alle Schaufelreihen b mit der Geschwindigkeit U = 1 in positiver 7-Richtung und alle 
Schaufelreihen a mit der Geschwindigkeit U = —1 in negativer 7-Richtung bewegt werden. 
Uberraschend ist der Umstand, daB die Geschwindigkeitsverteilungskurven des Einzelgitters und 
des vielstufigen Gitters im Rahmen der Rechengenauigkeit (-- 0,5°4) zusammenfallen. Es folgt 
daraus, daB die gegenseitige Beeinflussung nebeneinander liegender Gitter bei der Verdrangungs- 
es sehr gering ist. Kin ahnliches Ergebnis war fiir das Tandemplattengitter gefunden 
worden. 

In Abb. 9 ist die Geschwindigkeitsverteilung der Zirkulationsstrémung fiir eine Platte des 
Gitters b dargestellt. Es sind wieder zum Vergleich eingetragen die entsprechenden Geschwindig- 
keitsverteilungen des Einzelgitters gleicher Teilung und der Einzelplatte. Die Zirkulation ist fiir 
alle drei Kurven gleich, /’= 1,21060. Der Einflu8 der benachbarten Platten bewirkt beim 
Kinzelgitter eine Abnahme der Geschwindigkeit in Schaufelmitte gegeniiber der Einzelplatte. Bei 
dem vielstufigen Gitter wird eine weitere Geschwindigkeitsabnahme in der Schaufelmitte dadurch 
verursacht, daB die Strémung an den Schaufelenden auf den Spalt eingeschniirt wird. Dies fiihrt 
zu einer Geschwindigkeitszunahme an den Schaufelenden und damit wegen der konstanten Zirku- 
lation zu einer Geschwindigkeitsabnahme in Schaufelmitte. 

Als letztes ist in Abb. 10 die Geschwindigkeitsverteilung aufgetragen von je einer Platte der 
Gitter a und der Gitter 6, und zwar fiir den Fall, daB die Joukowsky-Bedingung erfiillt ist. Die 
Geschwindigkeitsverteilung Abb. 10 ist berechnet als Differenz w,,—w,, der Geschwindigkeits- 
verteilungen von Abb. 8 und Abb. 9, die beide an den Plattenhinterkanten unendlich sind. Daher 
ist die Geschwindigkeitsverteilung Abb. 10 in den letzten 10% der Plattentiefe (an den Platten- 
hinterkanten) unsicher. 

Eigentiimlich ist das Uberschneiden der Geschwindigkeitsverteilungskurven in der Nahe der 
Plattenhinterkanten. Diese Erscheinung trat schon beim Tandemgitter auf**. Die allgemeinste 
Geschwindigkeitsverteilung ergibt sich, wenn man der Geschwindigkeitsverteilung Abb. 10 noch 
eine Parallelstrémung parallel zur -Achse tiberlagert. Zum Vergleich ist die Geschwindigkeitsver- 


- teilung des Einzelgitters gleicher Teilung eingetragen, das ebenso wie das vielstufige Gitter mit 


der Geschwindigkeit U in dem ruhenden Medium in Gitterfrontrichtung bewegt wird. Aus den 
Geschwindigkeitsverteilungen geht hervor, da® die Zirkulation des vielstufigen Gitters kleiner 
ist als die des Einzelgitters. Nach (9) ist die Zirkulation einer Platte des vielstufigen Gitters 
TI = 1,21060 (bzw. —1,21060), was einem Verhaltnis J‘/J\, = 0,545 entspricht ([% Zirkulation 
der Einzelplatte). Fir das Einzelgitter ist [’/[) = 0,724. 

Zum SchluB sei noch angemerkt, daB auf ahnliche Weise der Fall berechnet werden kann, da8 
die beiden Gitterreihen a und b des vielstufigen Gitters um 1/2 gegeneinander versetzt sind. 


(Eingegangen am 17. September 1960.) 
Anschrift des Verfassers: Dr. E.-G. Feindt, Hamburg-Garstedt, Breslauer Str. 13. 


* Vel. E.-G. Feindt, a. a. O. 
** F-G. Feindt, a. a. O. 
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Uber Messungen an periodisch instationar angestromten Tragfliigeln* 
Von H. Schwanecke 
Mitteilung Nr. 247 der Versuchsanstalt fiir Wasserbau und Schiffbau, Berlin 


1. Einleitung. Die in dem vorliegenden Bericht mitgeteilten Ergebnisse von Messungen an in- 
stationdr angestrémten Tragfliigeln sind Teil eines Forschungsvorhabens mit dem Ziel der theore- 
tischen und experimentellen Ermittlung der instationaren hydrodynamischen Kraftwirkungen an 
schwingenden Propellern?, Propellern in peripher ungleichférmiger Strémung®*® und an Wasser- 
tragfliigeln unter Seegangshedingungen (Tragfliigelboote). Die experimentellen Untersuchungen, 
die an ebenen rechteckigen Tragfliigeln mit Plattenprofil durchgefiihrt wurden, waren erforderlich, 
um grundsatzliche Fragen hinsichtlich der GréBe und des zeitlichen Verlaufs der Normal- und 
Tangentialkrafte sowie des Strémungsverlaufes in der Umgebung der Tragfliigel zu klaren. Aufer- 
dem war es notwendig, theoretische Untersuchungen, die an Tragfliigeln groBer Spannweite bei 
periodisch instationdrer Anstrémung durchgefiihrt wurden, experimentell nachzupriifen. 

Die Messungen sollten sich in der Hauptsache auf die Vorgainge an einem feststehenden Trag- 
fliigel in periodisch instationarer Strémung erstrecken. Fiir den Vergleich mit den theoretischen 
Rechnungen war es erwiinscht, ein Geschwindigkeitsfeld zu erzeugen, in dem zeitlich veranderliche 
Geschwindigkeitskomponenten nur in der zur Tragfliigeloberflache senkrechten Ebene auftraten, 
die Geschwindigkeitskomponenten in der Ebene des Tragfliigels jedoch konstant blieben. Ferner 
sollte das instationdare Strémungsfeld zweidimensional sein. Wird ein raumfestes rechtwinkliges 
Koordinatensystem (x, y, 2) eingefiihrt, dessen positive x-Achse in Richtung der mittleren Strémungs- 
- geschwindigkeit und dessen positive z-Achse nach oben zeigt, so miissen in einem instationaren 
Geschwindigkeitsfeld mit den erwahnten Eigenschaften die Komponenten des Geschwindigkeits- 
vektors v(t) den Bedingungen geniigen 


0, = U,) = Konst., 
Vy = OF: 
v, = |w| sin (vt —yvt,). 


Wegen vy = konst. kann fiir die Komponente v, gesetzt werden. 


Dias || sin ( fee ; 


Fiir die Bezeichnungen siehe man Abschnitt 5,c). 


2. MeSapparatur und Messungen. a) GroBer Schlagfliigel. Ein instationares Geschwindig- 
keitsfeld, das den gewiinschten Bedingungen entspricht, kann mit einem zwischen zwei senkrechten 
Wanden schwingenden Schlagfliigel, das ist ein Tragfliigel, der translatorische Bewegungen senk- 
recht zu seiner Ebene ausfiihrt, mit ausreichender Genauigkeit erzeugt werden [vgl. Abschnitt 5 (a)]. 
Hierbei ist die Konstruktion des Schlagfliigelantriebes mit groBer Sorgfalt auszufiihren, da eine 
reine Sinus-Bewegung erzielt werden mu. Aus diesem Grunde wurde eine Antriebsanlage mit 
Exzenterscheiben, Gleitsteinen und vertikalen Schubstangen gewahlt, wie aus der Skizze und den 
Fotografien der Abb. 1 zu ersehen ist*. Weitere Antriebselemente waren ein Schneckengetriebe 
(i = 10,3), ein stufenlos regelbares Getriebe (sog. H-Trieb, mit einem Gesamtregelbereich 1:9) und 
ein elektronisch gesteuerter Drehstrommotor (1,5 PS bei 1400 U/min). Die seitlichen parallelen Leit- 
wande, die sich vom Tankboden bis etwa 150 mm iiber die Wasseroberflache erstreckten, befanden 
sich in etwa 2mm Abstand von den Seitenkanten des Schlagfliigels, der eine Breite von b = 1500 mm 
hatte. Bei einer Breite der MeBstrecke des Umlauftanks von ca. 1800 mm verblieben damit an 


* Mit Unterstiitzung der Deutschen Forschungsgemeinschaft durchgefiihrt 
1H. Schwanecke, Schiffstechnik 7 (1960) S. 170" pepe 
2 S. Schuster Z. VDI. 89 (1955) S. 1781. 
y ee Isay, Schiffstechnik 5 (1958) S. 157. 
- Schwanecke, Untersuchung von Tragfliigeln in ungleichformiger Stré als raumlich 
Forschungsvorhaben DFG Schu 24/16 u. 24/21, Bericht Nr (1959) (anverdtfentiiclity: sek Oe Se 
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faktor der Einlaufdiise im Verhaltnis 83.5 
in: der MeBstrecke wurden bei einer derarti 


Abb. 1. Versuchseinrichtung. 


a) Schema. 
b) u. c) Schlagfliigel mit Antrieb. 
d) Tragfliigel. 


Bezeichnungen: a Einlaufdiise des Umlauftanks; b Leitwande; 
¢ Schlagfliigel; d Antriebsmotor; e Getriebe (H-Trieb, stufenlos 
regelbar); f Schneckentrieb; g Antriebswelle; h Exzenterscheiben; 
i Schubstangen; k Tragfliigel; 1 Sechskomp.-Waage. d 


gerufen. In Strémungsrichtung erstreckten sich die Leitwande bis etwa 1000 mm hinter die Aus- 
trittskante des Schlagfliigels. Ein Teil der auf der verglasten Seite der Mefstrecke befindlichen 
Leitwand wurde aus Plexiglas hergestellt, um fotografische Aufnahmen zu erméglichen. 

Die Profillange des aus einer Stahlkonstruktion (Profilstahl mit Mantelblech) bestehenden 
Schlagfliigels betrug 533 mm, die gréBte Dicke 32 mm (d/l = 0,06), als Profil wurde das in Abb. 2a 
angegebene Plattenprofil verwendet. Die Antriebsanlage wurde so ausgelegt, daB bei einer mittleren 
Strémungsgeschwindigkeit in der MeBstrecke von vy = 2 m/s und einer Abwindamplitude von 
|20 = 0,25 vy = 0,5 m/s mit Schlagfrequenzen bis 6 Hz gearbeitet werden konnte. Es konnten 
damit reduzierte Frequenzen fiir den Schlagfliigel von maximal fy = y 1/2 v5 = 5 erreicht werden. 
Um bei allen reduzierten Frequenzen eine Abwindamplitude ungefahr gleicher GréBe zu erhalten, 
war es notwendig, mit steigender reduzierter Frequenz die Schlagamplitude zu verkleinern [vgl. 
Abschnitt 5,a)]. Aus diesem Grunde wurden die Gleitsteine auf den Exzenterscheiben so weit ver- 
schiebbar angeordnet, daB die Schlagamplitude zwischen 0 und 80 mm verandert werden konnte. 

Bei den Versuchen hatte die mittlere Strémungsgeschwindigkeit in der Mefstrecke die Grobe 
Vv) = 1,3 m/s, mit Ausnahme der Messungen bei der kleinsten reduzierten Frequenz; hier wurde sie 
auf v) = 1,4 m/s erhéht. Die Schlagamplituden wurden zwischen 79 mm bei der kleinsten (u = 0,75) 
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und 7,6 mm bei der gréBten reduzierten Frequenz (u = 4,5) variiert. Die Abwindamplituden hatten 
dann im Mittel etwa die GréBe jw) ~ 0,16 vy. Diese Reduktion der Versuchsdaten gegeniiber den 
Daten der Auslegung erwies sich als notwendig, da bei den Messungen an dem feststehenden Trag- 
fliigel die Gefahr des Uberlastens der Sechskomponentenwaage bestand. Die Wassertiefe am Ort 
des Schlagfliigels betrug 1140 mm, die Nullstellung des Fliigels lag etwa 600 mm iiber dem Tank- 
boden. Zur Messung der Schlagfrequenzen diente ein elektrisches Zahlwerk, das von einem an einer 
der vertikalen Schubstangen befestigten Kontakt betatigt wurde. Die mittlere Strémungsge- 
schwindigkeit wurde mit einem etwa 0,8m vor der Schlagfliigeleintrittskante angeordneten 


Prandtlrohr gemessen. 
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a > 
Abb. 2. Plattenprofile. 
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b) Kleiner Schlagfligel. Der groBe Schlagfliigel diente nur zum Erzeugen eines periodisch 
instationéren Strémungsfeldes. Zu den beabsichtigten Messungen des Phasenwinkels zwischen 
Schlagbewegung und Normalkraftverlauf, wozu Dehnungsmefstreifen verwendet werden sollten, 
war er infolge seiner Konstruktion nicht geeignet. Mit diesen Messungen sollten die Ergebnisse der 
Theorie des schwingenden Fliigels im Bereich kleiner reduzierter Frequenzen nachgeprift werden, 
was fiir die Anwendung dieser Theorie auf den Fall eines schwingenden Schiffspropellers von Be- 
deutung ist. Es wurde hierfiir aus einer 10 mm starken Stahlplatte ein weiterer Schlagfliigel mit den 
Abmessungen b = 1500 mm, | = 150 mm, d = 10 mm (d/l = 0,067) angefertigt. Das Profil, eben- 
falls ein Plattenprofil, ergab sich aus dem Profil fiir den groBen Schlagfliigel (Abb. 2a) durch Ver- 
zerren der Ordinaten im Verhaltnis 6,7:6. Um die fiir die in der Fliigelmitte angebrachten Deh- 
nungsmefstreifen notwendigen elastischen Verformungen bei der Schlagbewegung zu erreichen, 
wurde der Tragfliigel an den Schubstangen mit Zapfen gelagert. Auf der einen Seite wurden zwei 
Zapfen angebracht, um das Moment der hydrodynamischen Normalkraft um die waagerechte 
Achse aufzunehmen und auf der anderen Seite nur ein Zapfen in der Nahe des vorderen Neutral- 
punktes des Tragfliigels. Diese Anordnung der Zapfen sollte verhindern, daB bei der Schlaghewegung 
infolge nicht parallelen Laufes der Schubstangen zusatzliche Biege- oder Torsionsspannungen im 
Fligel auftreten, wodurch die Anzeige der Dehnungsmefstreifen verfalscht worden ware. 

Bei den Versuchen wurden die mittlere Strémungsgeschwindigkeit zwischen vy = 0 und v, 
= 3,1 m/s und die Schlagfrequenzen zwischen f = 0,6 Hz und f = 1,5 Hz variiert. Schlagfrequenzen 
und Strémungsgeschwindigkeiten wurden so aufeinander abgestimmt, daB reduzierte Frequenzen 
zwischen 44 = 0,08 und 4 = oo erreicht werden konnten. Dem Bereich kleiner reduzierter Frequen- 
zen wurde besondere Aufmerksamkeit gewidmet. Die Schlagamplitude betrug bei der Mehrzahl 
der Messungen 7,5 mm. Zum Uberpriifen einer méglichen nichtlinearen Abhangigkeit des zeitlichen 
Verlaufes der Normalkraft von der Amplitude wurden vereinzelte Messungen mit Schlagamplituden 
der GréBe 14,7 mm und 30 mm durchgefiihrt. 


Die Dehnungsmefistreifen wurden auf beiden Seiten des Schlagfliigels angeordnet und mit einer 
Tragerfrequenz-MeBbriicke betrieben, die vorher auf mégliche Phasenfehler bei der Anzeige tiber- 
priift wurde. Als Registriergerat diente ein Schleifen-Oszillograph mit Schleifen hoher Eigenfrequenz. 
Die Messungen wurden in Luft und bei den oben angegebenen Variationen von Frequenz, Schlag- 
amplitude und mittlerer Strémungsgeschwindigkeit im Wasser durchgefiihrt. 
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c) Schlagfliigelabwind. Zum Aufmessen der Amplituden des instationaren Strémungsfeldes 
hinter dem groBen Schlagfliigel diente eine Platte mit drei Reihen von je acht Bedensonden, wie sie 
aus der Abb. 3a zu ersehen ist. Die Faden bestanden aus Wolle, waren hinreichend diinn and hatten 
-eine Lange von etwa 50 mm. Thr gegenseitiger Abstand in Strémungsrichtung betrug 55 mm, die 


a 


a) MeBeinrichtung (Fadensonden). 
b) UnbeeinfluBte Strémung. 
ce) Durch den Tragfliigel beeinfluBte Strémung. c 


: Abb. 3. Bestimmung der Amplituden der periodisch instationaren Strémung. 
Reduzierte Frequenz 4 = 1,13; Fortschrittsgeschwindigkeit v) = 1,29 m/sek. Datum: 11. 11.59 u. 20. 1. 60. 


T~500mA 
bei flieBendem 


597 ra 
Sohleife HKKL YI Wasser (= 71sek') 


Emptindlichkeit: 
= 05-0 A/mm Sonde 
Lichtpunkt-Linien-Schreiber +8 


Abb. 4. Hitzfilm-Sonde (VWS). 


a) Ansicht der Sonde. 
b) Stirnfliche der Sonde (3 x 4mm) mit MeBelement (2,5 4 Pt-Folie). 
c) Schaltschema. 


drei Reihen waren je 40 mm voneinander entfernt. Die Platte mit den Fadensonden wurde 250 mm 
neben der Mittelebene des Schlagfliigels angebracht. Sie konnte in Strémungsrichtung so weit 
verschoben werden, daB es méglich war, mit den Fadensonden einen Bereich von 100 mm bis 
830 mm hinter der Austrittskante des groBen Schlagfliigels zu erfassen. Die Platte wurde mit einer 
Kamera mit Teleobjektiv bei Dunkelfeldbeleuchtung und etwa 2—3 s Belichtungszeit von der Seite 
fotografiert. Der groBe Vorteil der Zeitaufnahmen bestand darin, daB sich infolge der Abnahme 
der Geschwindigkeit auf den Wert Null in den Endlagen, diese besonders deutlich markierten. 
Dadurch konnte die Amplitude der periodischen Stérgeschwindigkeiten leicht bestimmt werden. 
AuBerdem wurde durch die Zeitaufnahmen eine Mittelwertbildung erreicht. Mit den Fadensonden 
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Abb. 5. Phasen- u. Nullpunktsdifferenzen der verwendeten MeBwerke (f, ~ 6 Hz) gegeniiber MeBwerk 1(f, > 1 Hz). 
Datum: 25. 1. 60. 
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Abb. 6. Richtungsempfindlichkeit der Hitzfilm-Sonde. Datum: 26. 2. 1960. 


wurden sowohl die Amplituden des unbeeinfluBten periodisch instationaren Strémungsfeldes 
(Abb. 3b) als auch die des durch den festen Tragfliigel gestérten Strémungsfeldes (Abb. 3c) ermittelt. 

Zum Aufmessen des zeitlichen Verlaufes des instationdren Strémungsfeldes wurde eine fiir diesen 
Zweck speziell angefertigte Hitzfilmsonde (Abb. 4) verwendet. Als MeBelement diente eine Platin- 
folie mit den Abmessungen 4 x 0,2 x 2,5- 103mm. Die Sonde wurde in der in Abb. 4c angegebenen 
Weise elektrisch geschaltet. Zur Aufnahme der Signale diente ein Lichtpunktlinienschreiber mit 
hochempfindlichen Galvanometern. Diese MeBwerke hatten infolge ihrer hohen Empfindlichkeit 
nur eine Eigenfrequenz von etwa fy = 6 Hz, so daB es notwendig war, ihre Phasenfehler zu be- 
stimmen. Hierzu wurden Vergleichsmessungen mit einem Galvanometer geringer Empfindlichkeit 
und hoher Eigenfrequenz durchgefiihrt (Abb. 5). 

Eine in der angegebenen Weise betriebene Hitzfilmsonde hat bekanntlich eine Geschwindigkeits- 
Widerstandscharakteristik, die sich mit zunehmender Geschwindigkeit einem Grenzwert nahert. 
Infolgedessen treten im Bereich héherer mittlerer Geschwindigkeiten auch bei relativ starken 
Geschwindigkeitsinderungen keine nennenswerten Widerstandsunterschiede mchr auf. Diese 
Kigenschaft der Sonde konnte zum Bestimmen des zeitlichen Verlaufes der periodischen Zusatz- 
geschwindigkeiten des Strémungsfeldes ausgenutzt werden. Das PlatinmeBelement wurde zu diesem 


XXX. Band 1961 H. Schwanecke: Messungen an periodisch instationaér angestrémten Tragfliigeln 355 


Zweck auf der abgeplatteten Stirnflache des Tragerkérpers befestigt (Abb. 4 b). Bildet die Stirnflache 
mit der Richtung der Anstrémung einen rechten Winkel, so liegt das MeBelement in bzw. in der Nahe 
der Staupunktes. Andert sich nun die Strémungsrichtung um einen gewissen Winkel, so wandert der 
Staupunkt, und das MeBelement gerat mit zunehmendem Winkel in den Bereich héherer Strémungs- 
geschwindigkeiten. Da im Staupunktbereich die értlichen Geschwindigkeiten klein sind, ist die 
Empfindlichkeit der Sonde hier verhaltnismaBig groB. Die Abb. 6 zeigt Messungen, die mit der 
Sonde bei konstanter Stromungsgeschwindigkeit bei verschiedenen Anstellwinkeln durchgefiihrt 
wurden. Man sieht hier deutlich die starke Winkelabhangigkeit bis zu etwa 40°. Da die Sonde nur 
Anderungen der Strémungsgeschwindigkeit anzeigt, gleichgiiltig ob es sich dabei um eine F olge 
von Winkelanderungen oder um Anderungen der mittleren Strémungsgeschwindigkeit handelt, sind 
die Ergebnisse nur qualitativ zu verwerten, was aber fiir den beabsichtigten Zweck ohne Bedeutung 
war. Beim Betrieb wurde die Sonde an einem aus Stahlblech bestehenden Trager befestigt (Abb. 7), 
der in Strémungsrichtung verschoben werden konnte. Es war dadurch méglich, in dem oben 
angegebenen Bereich hinter dem groBen Schlagfliigel in jedem beliebigen Punkt den zeitlichen Verlauf 
der periodischen Zusatzgeschwindigkeiten zu bestimmen. 


Hitz film 


Tragtlige/ 
Schlag tlige/ 

Anstromung 

—— 


a 


Abb. 7. Anordnung der Hitzfilmsonde. 


a) Schema. 
b) Sonde in MeBposition. 
c) Sonde mit Trager. 


d) Feststehender Tragfliigel. Der feststehende Tragfliigel hatte eine Breite b = 1000 mm, 
eine Profillange 1 = 200 mm, eine Dicke d= 10 mm (6/l = 5, d/l = 0,05) und wurde aus einer 
Stahlplatte mit dem in Abb. 2b angegebenen Profil (Plattenprofil) hergestellt. An dem Fliigel wurde 
eine freistehende Mittelstiitze angeordnet, mit der er an dem MeBglied einer Sechskomponenten- 
waage befestigt werden konnte. Um die Stiitze dem Einflu8 der Strémung zu entziehen, wurde sie 
innerhalb einer profilierten Blechverkleidung angeordnet, deren Unterkante sich etwa 5 mm iiber 
dem Tragfliigel befand und die an ihrem oberen Ende mit dem Gehause der Sechskomponentenwaage 
verbunden wurde. Um die Schwingungsfahigkeit des Systems zu vermindern, wurden die Ver- 
kleidung und der Tragfliigel in geeigneter Form mit Drahten verspannt (Abb. ld). = 

Zur Messung der an dem feststehenden Tragfliigel auftretenden Kraftwirkungen diente die in 
der VWS entwickelte kleine Sechskomponentenwaage, iiber die in der Literatur schon verschiedent- 
lich berichtet wurde (vgl. z. B.1). Zur Registrierung der Signale der Sechskomponentenwaage wurde 
der schon erwahnte Schleifen-Oszillograph mit MeBwerken hoher Eigenfrequenzen verwendet. 

Die zu den induktiven Komponentengebern der Sechskomponentenwaage gehorenden Trager- 
frequenz-MeBbriicken waren zunachst auf mégliche Phasenfehler zu untersuchen. Hierzu eel aiaie 
zeitliche Bewegungsverlauf einer der Schubstangen des Schlagfliigels mit Hilfe eines bie oe a 
standes und mit Hilfe eines als Beschleunigungsgeber verwendeten induktiven MeBwert-Gebers, der 
mit einer der MeBbriicken verbunden wurde, mit dem Schleifenoszillographen ae 
Daraus konnte dann der Phasenfehler der MeBbriicke bestimmt werden (Abb. 8). Da der konstruk- 


1 H. Schwanecke, Schiffstechnik 6 (1959) S. 93. is 
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tive Aufbau bei allen Tragerfrequenzbriicken der gleiche war, wurde die Uberpriifung nur eines 
Gerates als ausreichend angesehen. 

Der Abstand zwischen Schlagfliigelaustrittskante und Tragfliigeleintrittskante betrug bei kleinen 
reduzierten Frequenzen 100 mm und bei groBen reduzierten Frequenzen 550 mm. Die groBe Ent- 
fernung muBte gewahlt werden, weil die instationdéren Kraftwirkungen bei dem kleinen Abstand 
zwischen Schlagfliigel und festem Fliigel mit gréBer werdenden reduzierten Frequenzen so groB wurden, 
daB der Tragfliigel einschlieBlich der Sechskomponentenwaage und der Haltebiithne zu schwingen 
begann. Da aber das instationaére Strémungsfeld hinter dem Schlagfliigel in geniigend vielen 
Punkten aufgemessen wurde, hatte die Veranderung des Abstandes zwischen Schlagfliigel und 
festem Fliigel auf die Giite der Messungen keinen EinfluB. 

Im vorliegenden Fall war der Strémungszustand zweidimensional, so daB nur eine Dreikompo- 
nenten-Messung vorgenommen zu werden brauchte. In dem instationdren Geschwindigkeitsfeld 
waren in erster Naiherung nur Zusatzgeschwindigkeiten parallel der 2-Achse mit einem nahezu 
reinen Sinus-Verlauf vorhanden. Die Tragfliigel-Normalkraft hatte daher keinen zeitlichen Mittel- 
wert. Infolgedessen war es méglich, mit Hilfe der Signale von den drei in der Vertikalebene 
liegenden Komponentengebern der Sechskomponentenwaage Normalkraft, Moment der Normal- 
kraft um eine Achse parallel der Fliigeleintrittskante und Tangentialkraft (induzierter Widerstand) 
des Tragfliigels zu messen. Die Abb. 9 zeigt einen Schrieb des Schleifenoszillographen mit den 
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Abb. 8. Phasenfehler der MeBbriicke 1 (Kompon. ,,a“*) beziiglich der Synchronisierungsmarken. Datum: 1. 2. und 8. 2. 60. 


Signalen der Komponentengeber a, b und c. AuBerdem ist ein 100 Hz-Signal als Zeitmarke, mit 
kurzzeitigen Unterbrechungen durch den Drehzahlkontakt des Schlagfliigelantriebes, auf dem Os- 
zillogramm zu erkennen. Die Drehzahlmarke wurde gleichzeitig zur Synchronisierung der Signale 
der Sechskomponentenwaage und der Hitzfilmsonde verwendet. 

Die Abb. 10 zeigt die mit einem Lichtpunktlinienschreiber aufgenommenen Signale der Kom- 
ponente a, der Hitzfilmsonde und des Drehzahlkontaktes bei zwei verschiedenen Schreibegeschwin- 
digkeiten des Registriergerates. Man sieht deutlich die mit der Hitzfilmsonde aufgenommenen 
zeitlichen Schwankungen der Anstrémrichtung. Die zeitliche Lage der Extremwerte lat sich 
daraus leicht bestimmen. Die Stirnflache der Sonde lag mit der Vorderkante des feststehenden 
Fliigels in einer Vertikalebene (Abb. 7a). Die Unsymmetrie des Hitzfilm—Sonden-Signales ist darauf 
zuriickzufiihren, daB die Sonde, je nach der Anstrémrichtung, einmal in der nahezu ungestérten 
Zustrémung und das andere Mal in der durch die Tragfliigelvorderkante beeinfluBten Strémung lag. 


3. Ergebnisse. a) GroBer Schlagfliigel. Mit Hilfe der in Abschnitt 2, c) beschriebenen Faden- 
sonden wurden die dimensionslosen Abwindamplituden w = |w|/v9 des groBen Schlagfliigels fiir 
reduzierte Frequenzen 4“ = 0,75 — 4,5 in 100 mm bis 830 mm Entfernung von der Schlagfliigel- 
austrittskante bestimmt. Eine Darstellung der Abwindamplituden, bezogen auf eine dimensionslose 
Schlagamplitude der GréBe 1, befindet sich in der Abb. 11. Zum Vergleich wurden die mit Hilfe 
der Formeln des Abschnittes 5, a) ermittelten theoretischen Abwindamplituden fiir € — 1,62 und 
€ = co in die Abb. 11 eingetragen. ee 

Es zeigt sich hier eine starke, mit wachsender reduzi i 
der MeBwerte von den theoretisch berechneten Werten, pee ee i enon a She asane 
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Das mit wachsendem Abstand von der Austrittskante durch den turbulenten Strémungszustand 
in der MeBstrecke bedingte Abwandern eines Teiles der in der Wirbelschleppe des Schlagfliigels 
enthaltenen Energie in die iiber und unter ihr befindlichen Strémungsbezirke und der dadurch 
hervorgerufene Abbau der Abwindamplituden auf der Spur des Schlagfliigels macht sich, wie die 
Auswertung der Abwindfotografien zeigte, erst in einer Entfernung von etwa 600 bis 700 mm 
hinter dem Schlagfliigel, d. h. fiir € > 2,4 starker bemerkbar. Die Abweichung der MeBwerte von 
den theoretischen Werten muf daher andere Ursachen haben. Auch scheidet die bei den Versuchen 
nicht eingehaltene Bedingung kleiner Schlagamplituden, die den theoretischen Rechnungen zu- 
grunde liegt, als Fehlerquelle aus, da die starksten Differenzen bei hohen reduzierten Frequenzen 
auftreten, also gerade dort, wo mit geringen Schlagamplituden gearbeitet wurde. Die Abweichungen 
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Abb. 9. Dreikomponentenmessung. Datum: 14. 12. 59. 
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Abb. 10. Kraftverlauf bei periodisch veranderlicher Anstrémung. 


der Abwindamplituden lassen sich daher wohl nur mit einer verminderten Zirkulation am Schlag- 
fliigel selbst erklaren. Die Wassertiefe iiber bzw. unter dem Schlagfliigel (540 mm baw. 600 ee 
war nicht viel gréBer als die Profillange des Schlagfliigels, so dafs mit einem verhaltnismabig starken 
Oberflachen- bzw. Bodeneinflu® auf die Druckverteilung gerechnet werden mu. Offenbar liegt 
hierin der Grund fir die schlechte Ubereinstimmung zwischen Theorie und Experiment. Eine 
nahere Untersuchung dieser Einfliisse wurde nicht vorgenommen, da hier in erster Linie die tat- 
sachlich vorhandenen Abwindamplituden des Schlagfliigels von Interesse waren. 


i iu i i der Phasenwinkel zwischen der 

b) Kleiner Schlagfliigel. Die Ergebnisse der Messungen 
ir ease und hao zeitlichen Verlauf der Normalkraft sind in der Abb. 12 dargestellt. Sie 
wurden mit Hilfe eines harmonischen Analysators aus den Amplituden und Phasendifferenzen der 
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an dem in Luft und im Wasser schwingenden Schlagfliigel gemessenen Kraftverlaufe bestimmt. 
Die Ubereinstimmung zwischen Theorie und Versuch kann als gut bezeichnet werden. Eine Ab- 
hangigkeit von der Schlagamplitude scheint bei geniigend kleinen Amplituden nicht zu bestehen. 
Die Abweichungen bei der groBen Schlagamplitude (z)/1 = 0,20) ist méglicherweise die Folge eines 
Nullpunktfehlers, da sowohl bei der kleinen als auch bei der groBen reduzierten Frequenz die Betrage 
der Abweichungen etwa gleich sind, die Tendenz j edoch mit der Theorie iibereinstimmt. 


c) Feststehender Tragfligel. J. Kraftwirkungen. Zunachst war es erforderlich, die Original- 
aufzeichnungen der Komponentengeber der Sechskomponentenwaage zu linearisieren, da wegen der 
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Abb. 12. Phasenwinkel zwischen Weg- u. Normalkraft-Verlauf, A = oo. vy = 0 + 3,14 m/sek. Datum: 1.—2.2. 60. 


GréBe der Krafte zum Teil nicht mehr im linearen Bereich der Geberelemente gearbeitet werden 
konnte. Im Falle des in Abb. 9 dargestellten Kraftverlaufes war eine derartige Linearisierung bei 
der Komponente a erforderlich. Diese MaBnahme wirkte sich derart aus, daB die positiven Ampli- 
tuden entsprechend den Abweichungen der Gebercharakteristik von einer Geraden vergréBert und 
die negativen Amplituden verkleinert werden muBten. Die linearisierten Kraftverlaufe konnten 
dann mit Hilfe eines harmonischen Analysators nach GréBe und Phasenwinkel analysiert werden 
Mit den Ergebnissen der Analyse erfolgte die Bestimmung der Normalkraftamplituden, der Amaplis 
tuden des Momentes der Normalkraft um die Fliigelvorderkante und der mittleren Tancentialk fale 
(induzierter Widerstand). Sie wurden in der dimensionslosen Form 
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in den Abb. 13, 14 und 15 in Abhangigkeit von der reduzierten Frequenz dargestellt. Die Ergebnisse 
theoretischer Berechnungen [Abschnitt 5, b), Gl. (22), Gl. (23), (23a) und Gl. (25)] wurden eben- 
falls in die Abbildungen eingetragen. 

Die gemessenen Normalkraftamplituden (Abb. 13) liegen in dem ganzen erfaBten Bereich der 
reduzierten Frequenz um einen nahezu konstanten Betrag iiber den fiir groBe Spannweite geltenden 
theoretisch berechneten Werten, obwohl bei dem endlichen Seitenverhialtnis des Tragfliigelmodelles 
eher das gegenteilige Verhalten zu erwarten war. Selbstverstandlich ist hier mit eincr gréBeren 
Ungenauigkeit im Vergleich zu Messungen am stationar angestrémten Profil zu rechnen, schon 
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Abb. 13. Dimensionslose No-malkraftamplituden ae eee 1,29 = 1,44 m/sek. Datum: 14. 12.—21. 12. 59. 
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Abb. 14. Dimensionsloses Moment der Normalkraft, bezogen auf die Eintrittskante es Vo = 1,29 = 1,44 m/sek. 


Datum: 14. 12.—21. 12. 59. 


allein deswegen, weil in der vorliegenden Form die MeBergebnisse als Quotienten von je zwei ge- 
messenen Amplituden (Kraft- und Normalgeschwindigkeitsamplitude) auftreten. Hierdurch lassen 
sich nur starkere Streuungen, keineswegs aber nahezu konstante VergréBerungen der mittleren 
Werte, erklaren. Eine leichte VergréBerung der Normalkraft kann u. U. auf die endliche Profildicke 
zuriickgefithrt werden. Ganzlich ungeklart ist der Einflu8 nichtlinearer Zusammenhinge zwischen 
Zirkulation und induzierter Geschwindigkeit, auf die bei den theoretischen Rechnungen keine 
Riicksicht genommen wurde, die aber bei den verwendeten hydrodynamischen Anstellwinkeln von 
Xp, & + 9° (|w|/v9~ 0,16) schon eine Rolle spielen kénnen. Auch 148t sich der EinfluB der zweifellos 
poiken periodischen Randwirbel auf die Druckverteilung nicht genau voraussagen. ; 

Das dimensionslose Moment (Abb. 14) der Normalkraft um die Eintrittskante des Tragfliigels 
konnte infolge meBtechnischer Schwierigkeiten nur fiir kleine reduzierte Frequenzen ermittelt 
werden und ist seinem Betrag nach gréfer als es die theoretischen Rechnungen erwarten a 
Diese Differenz kann jedoch mit der ebenfalls nicht mit der Theorie itibereinstimmenden Normal- 
kraft begriindet werden. 
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In der Abb. 15 ist die mittlere Tangentialkraft (zeitlicher Mittelwert des induzierten Wider- 
standes) als Funktion der reduzierten Frequenz dargestellt. Im Gegensatz zu den Ergebnissen 
der Normalkraft- und der Momentenmessung kann die Ubereinstimmung der mittleren Tangential- 
kraft mit den theoretisch berechneten Werten als gut bezeichnet werden. Allerdings ist hier die 
Auswertung einfacher, da es sich, wie schon gesagt, um einen zeitlichen Mittelwert handelt. Die 
Streuungen der MeBwerte diirften eine Folge der quadratischen Abhangigkeit von den Amplituden 
des instationdéren Strémungsfeldes sein. 

Das Verhaltnis von Tangential- und Normalkraft, die sogenannte Gleitzahl des Tragfliigels, die 
im Gegensatz zum stationaren Fall bei periodisch instationaren Vorgangen negative Werte hat, ist 
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Abb. 15. Dimensionslose Tangentialkraft (zeitlicher Mittelwert) - Vy = 1,29~+ 1,44 m/sek. 


Datum: 14, 12.—21, 12. 59. 
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Datum: 14, 12.—21. 12. 59. 


in der Abb. 16 dargestellt. Die starke Streuung der MeBwerte ist darauf zuruckzufiihren, daB die 
Tangentialkraft vom Quadrat, die Normalkraft dagegen nur von der ersten Potenz des hydro- 
dynamischen Anstellwinkels abhangig ist. Aus versuchstechnischen Griinden war es jedoch nicht 
moglich, immer mit den gleichen Werten w zu arbeiten. 


Die auf den hydrodynamischen Anstellwinkel bezogene instationare Gleitzahl, die in Abb. 17 
als Funktion der reduzierten Frequenz angegeben ist, weicht im Bereich kleiner ju stark von den 
theoretisch berechneten Werten [Gl. (28) mit Gl. (22) und (25)] ab, was auf den groBen Unter- 


schied zwischen den theoretisch berechneten und den gemessenen Normalkraftamplituden zuruck- 
zufiihren ist. : 
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2. Phasenwinkel. In der Abb. 18 sind die Phasenwiakel der mit der an der Tragfligeleintritts- 
kante angebrachten Hitzfilmsonde bestimmten Extremwerte der periodisch instationdéren An- 
strémung beziiglich des als Synchronisierung dienenden Drehzahlkontaktes des Schlagfliigelantriebes 
angegeben. Die Abstande x, sind die bei den Versuchen eingestellten Entfernungen zwischen 
Schlagfliigelaustritts- und Tragfliigeleintrittskante. 

Mit Hilfe der Ergebnisse der harmonischen Analyse des Kraftverlaufes am Tragfliigel wurden 
dann die Phasenwinkel der Kraftextremwerte beziiglich der Synchronisierungsmarken bestimmt. 
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Abb. 17. Instationire Gleitzahl finstl® (bezogen auf die Normalgeschwindigkeitsamplitude %). v = 1,29 + 1,44 m/sek. 
Datum: 14. 12.—21. 12. 59. 
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Abb. 18. Phasenwinkel des Q-Extremwerte der veranderlichen 
Anstrémung beziigilch der Synchronisierungsmarken. 
Datum: 6. 1. 1960. 


Sie sind in der Abb. 19 dargestellt, Aus den zeitlichen Kraftverlaufen und dem zeitlichen Verlauf 
der instationaren Anstr6mung konnten nunmehr die Phasenwinkel zwischen Anstrémungs- und 
Kraftverlauf ermittelt werden. Sie sind, getrennt fiir Normalkraft und Tangentialkraft, in den 
Abb. 20 und 21 dargestellt. Um einen Vergleich mit den theoretisch berechneten Werten zu ermbg- 
lichen [Gl. (22), (26)] war es notwendig, die in Abb. 18 angegebenen Phasenwinkel 0, die die zeit- 
lichen Anderungen der Anstrémgeschwindigkeit an der Tragfliigeleintrittskante darstellen, umzu- 
rechnen auf die Verhaltnisse am Ort der Profilmitte. Mit Ausnahme des Bereiches kleiner -Werte 
ist die Ubereinstimmung von Theorie und Experiment zufriedenstellend, 


362 H. Schwanecke: Messungen an periodisch instationaér angestromten Tragfliigeln —_—sIngenieur-Archiv 


3. Riickwirkung des Tragfliigels auf das Strémungsfeld. An der Oberflache des Tragfliigels ver- 
schwinden infolge seiner Undurchdringlichkeit die Normalkomponenten der Geschwindigkeit. In 
einem instationdiren Strémungsfeld mit periodischen Geschwindigkeitskomponenten normal zur 
Tragfliigeloberflache bedeutet das, da die Amplituden dieser Normalkomponenten mit zunehmen- 
dem Abstand von der Eintrittskante des Tragfliigels mehr oder weniger vermindert werden. Die 
Abbn. 3b u. 3c zeigten diesen Vorgang deutlich. Es tritt hierbei eine Energietransformation ein, derart, 
daB von dem Tragfliigel dem instationaren Strémungsfeld Energie entzogen wird, die sich in einem 
negativen Widerstand des Tragfliigels bemerkbar macht, d. h. in Vortriebsenergie umgesetzt wird. 
Ein MaB fiir diese Energieiibertragung ist das Verhaltnis der Amplituden des instationaren Stré- 
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Abb. 20. Phasenwinkel der Normalkraft Pu): vy = 1,29 + 1,44 m/sek. 
Datum: 14, 12.—21. 12. 59. 
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Abb. 21. Phascnwinkel der Tangentialkraft P,(u): V) = 1,29 = 1,44 m/sek. 
Datum: 14, 12,—21. 12. 59. 


mungsfeldes vor und hinter dem Tragfliigel. Man kann hier von einem instationdren Giitegrad des 
Profiles sprechen, der den Wert 1 erreicht, wenn dic Amplituden der Geschwindigkeitsnormalkom- 
ponenten vollstandig abgebaut werden und ihn in der Form 

Ninse = Li = 
definieren [Gl. (27)], wenn die Amplituden vor dem Tragfliigel mit w und die Amplituden hinter 
dem Tragfliigel mit w, bezeichnet werden. Die Abb. 22 zeigt die durch Messung bestimmten und 
die theoretisch berechneten Giitegrade. Wie aus dem Diagramm zu ersehen ist, liegen die durch 


a ll 
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Messung bestimmten Giitegrade bei kleinen reduzierten Frequenzen héher als die theoretisch 
berechneten Werte. Das deutet darauf hin, daB die Geschwindigkeitsnormalkomponenten bei 


kleinen yw-Werten starker abgebaut werden, als es mit den linearen Ansitzen der Theorie erfat 
wird. 


4. Zusammenfassung. Es wurde iiber Messungen an zwei schwingenden Tragfliigeln (Schlag- 
fliigel) in gleichférmiger Grundstrémung und an einem feststehenden Tragfliigel in periodisch in- 
stationdrer Strémung berichtet. Diese Messungen sollten einige grundsatzliche Fragen klaren, die 
im Zusammenhang mit der Ermittlung der an schwingenden Schiffspropellern wirksamen hydro- 
dynamischen Massen und Dampfungen, der hydrodynamischen Erregung von Schiffswellenanlagen 
infolge ungleichférmigen Zustromes zum Propeller und der hydrodynamischen Kraftwirkungen an 
den Tragwerken von Tragfliigelbooten bei Fahrt im Seegang auftreten. 

Die an einem zwischen zwei senkrechten Wanden schwingenden Tragfliigel durchgefiihrten 
Messungen des Phasenwinkels zwischen Bewegungs- und Kraftverlauf zeigten eine gute Uber- 
einstimmung mit den fiir diesen Fall theoretisch berechneten Werten. Dieses Ergebnis ist im Hin- 
blick auf die an schwingenden Schiffspropellern auftretenden hydrodynamischen Kraftwirkungen 
von Bedeutung, weil die GréBe des Phasenwinkels bei der Berechnung der hydrodynamischen 
Dampfungen und der hydrodynamischen Massen eine wesentliche Rolle spielt *. 

Messungen der GréBe und des zeitlichen Verlaufes der an einem ruhenden Tragfliigel mit dem 
Seitenverhaltnis A =5 bei periodisch instationdrer Anstrémung auftretenden Kraftwirkungen 
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Abb. 22. Instationarer Giitegrad. vy = 1,29 + 1,44 m/sek. Datum: 23. 9.—23. 11.59 und 19. 1.—21. 1. 60. 
ergaben hinsichtlich ihrer GréBe und ihres Verlaufes in Abhangigkeit von der reduzierten Frequenz 


im allgemeinen eine zufriedenstellende Ubereinstimmung mit theoretischen Rechnungen, die fiir 
groBe Spannweite des Tragfliigels (ebener Fall) durchgefiihrt wurden. Merkliche Differenzen 


- zwischen Theorie und Experiment zeigten sich bei allen reduzierten Frequenzen im Fall der Normal- 


kraftamplituden und im Bereich kleiner reduzierter Frequenzen, im Fall der Tangentialkraft sowie 
bei den Phasenwinkeln der Normal- und Tangentialkraft beziiglich der Anderungen der An- 
strémung. Fiir die, verglichen mit der Theorie, zu groBen Normalkraftamplituden kann keine 
befriedigende Erklarung gegeben werden, dagegen diirften die bei Kleinen Werten der reduzierten 
Frequenz festgestellten Abweichungen der tibrigen durch Messung bestimmten GréBen von den 
theoretischen Werten eine Folge des endlichen Seitenverhdltnisses des Tragfliigels sein. Dieses 
macht sich um so starker bemerkbar, je kleiner die reduzierte Frequenz ist?. ; se 
Die urspriinglich vorgesehenen Messungen mit systematisch verkleinertem Seitenverhaltnis 
des Tragfliigels konnten im Rahmen des im Abschnitt 2 beschriebenen Versuchsprogrammes nicht 
mehr vorgenommen werden. Sie werden jedoch, da fiir die Ubertragung der Ergebnisse der Ver- 
suche an einfachen Tragfliigeln auf die Fliigel eines Schraubenpropellers insbesondere kleine Seiten- 
verhaltnisse interessant sind, zu gegebener Zeit mit verbesserten MeBmethoden nachgeholt werden. 


* Vel. te 1 Seite 350. 
# ee Beet nuns des induzierten Geschwindigkeitsfeldes eines gebundenen Stabwirbels endlicher 


Lange mit zeitlich verdnderlicher Wirbelstarke bei gleichformiger Queranstromung, VWS-Bericht Nr. 140 
(1959) (unveréffentlicht). 
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Die experimentelle Bestimmung des instationéren Giitegrades des Tragfliigels, der ein Mab 
fiir die dem instationaren Strémungsfeld entnommene Energie ist, zeigte, daB schon bei verhaltnis- 
mabig Kleinen reduzierten Frequenzen, d.h. bei verhaltnismaBig groBen Wellenlangen der normal 
zum Tragfliigel gerichteten periodischen Zusatzgeschwindigkeiten, ein merklicher Energieentzug 
eintritt. Die experimentell bestimmten Werte liegen bei kleinen reduzierten Frequenzen zum Teil 
erheblich itber den mit Hilfe eines einfachen linearen Abklingungsgesetzes ermittelten Werten der 
Theorie. 

AbschlieBend méchte der Verfasser nicht versiumen, Herrn Prof. Dr. Schuster fiir seine ver- 
standnisvolle Férderung der Arbeiten, sowie Herrn Dipl.-Ing. Dernedde und Herrn Strube fir ihre 
Hilfe bei der Durchfiihrung und der Auswertung der Messungen zu danken. 


5. Anhang. a) Ergebnisse der ebenen Theorie des Schlagfligels in gleichférmiger 
Strémung. Die an einem Tragfliigel mit groBer Spannweite bei schwingender Bewegung auf- 
tretenden hydrodynamischen Kraftwirkungen kénnen mit Hilfe der bekannten Theorien von 
Theodorsen}, v. Karmén und Sears®, KiiBner? und KiiBner und Schwarz* unter den iblichen Vor- 
aussetzungen berechnet werden. Es werden hier nur die Ergebnisse und auch nur soweit sie fiir 
die im Hauptteil der vorliegenden Arbeit mitgeteilten experimentellen Untersuchungen von Be- 
deutung sind, wiedergegeben. 


Lautet das Bewegungsgesetz der Profilsehne des Schlagfligels 
z(t) =z, exp (re) 


so hat der Auftrieb cines aus einem Fligel groBer Spannweite herausgeschnittenen Streifens mit 
der Breite b in der vom Verfasser bevorzugten Kiissnerschen Darstellung die GréBe® 


P, = qb x20 u{T"(u) + w—i [lL + T’(w)]} exp (iv 8) (1) 
oder, in etwas verallgemeinerter Form, 
P,= qb % 2m F(u)| exp ( O(u) + 17 1) (1a) 
mit 
|F()| =" VL + + |T(w)P +2 (Tu) + wT’) 
a By) = — are tg ey ris y 


Das Moment des Auftriebes, bezogen auf den Sehnenmittelpunkt, hat die GréBe 
M,=q1b Zot {T’’ (wu) —ifl+ T’(u)]} exp (iv #). (3) 
Der induzierte Widerstand ergibt sich zu 


2 
P, = Pry — qb 0 2 {T"%u) —[1 + Tw)P—i2 Tu) 1+ TR exp(i2r 
mit dem zeitlichen Mittelwert 
-2 
Pom =— Gb | aye [l + 2 Tu) + | Ty)? I. (5) 
Der Vollstandigkeit halber wird auch der Ausdruck fiir die aufgenommene Leistung angegeben: 


L=L,, + qb? 20p? {1+ Tu) + ilu + T’(u)]} exp (i2 v2) (6) 


1 Th. Theodorsen, A general theory of aerodynamic instabilit d i 
Ratio. Gos. y y instability and the mechanism of flutter, NACA-Rep.- 


® Th. v. Karman u. W. R. Sears, Journ. Aero. Se. 5 (1938) S. 379. 


3 H.G. Kii eos i 2 spe es : a 
See asta ee Bericht tiber den instationaren Auftrieb von Fliigeln, Luftf.-For- 


4 H.G. KiiBner u. L. Schwarz, Der schwingende Fliigel mit i eli 
Forschg. 17 (1040) $337, gende Fliigel mit aerodynamisch ausgeglichenem Ruder, Luftf.- 


° Beziiglich der gewahlten Bezeichnungen vgl. Abschnitt 5c). 
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mit dem zeitlichen Mittelwert 
Be bY 7 
Lim = 4%) b = 20 [1 + T’(u)). 
‘Der Abwind des Schlagfliigels ergibt sich in bekannter Weise zu 
+1 ay oo 


S$ 


| [ 2h da—ie [ sede Gun) fC feds, exp (int) (8) 


x,=—1 x,=—l LX, 


2(x, = {T"'(u) — LT (tt )}} = aes a2 j1— x ; 


Das uneigentliche Interal in (8) hat die Lésung 


—E—x 


[ oak ee 2) (p= exp (—1 1 €) eet dt 
t 


= exp (—ipé) (6 [uv (E—m)] feat S.[u(§—x)]). 


Clue a = egg tL lea) | 


S,[u(E — x) = — St lee — a) 


gesetzt, so folgt aus (10) 


ioe) 


= [ SESE de temp (i fe — PE) i 6 weal +S le ai] 


E—x 


u(E — x) 


_ Damit lautet die Beziehung fiir den Abwind [vgl. (8)] 
+1 
— = — 3 py? exp (ivt—ipé) | 8(x1) exp (i ux) dx, 
0 


x, =—1 


41 
ss a (x,) exp (4 x) {f, [wu (& —*,)] —?f, [4 E — %) i a 


flw(é—a)) = Silu(E— =) —F + SEY 


flulé—s)1 = Glas) 


~ Fiir «(€ — x,) > 1 verschwinden die Ausdriicke (14) und der Abwind ergibt sich zu 
+1 
© — *0 y2 exp (iv t—i pe) | ee exp (i %) dx, « 
0 


pl 


| Mit g(x,) nach vt lautet die Lésung des Integralausdruckes in (15) 


Ms 8(%) exp (iu x,)dx,=% (Jo(4) Ip ye J 67) [E— T’(u)| 
a i eee te) T200)})- 


secant fn) \ i 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(16) 
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Nimmt die reduzierte Frequenz u groBe Werte an, so vereinfacht sich (16) wegen 


Tu), Tw)>?, 


cos («4 cos (“4 *) 
J(u) & —— —9 J) 2 fir u > 1 
% % 
Vez nz 
zu 
si os 
| g(x) exp (iu x,) dx, & 2 [— cos w + sin w —i (cos u + sin y)] . (17) 


Fir w > 1 und € 1 hat dann der Abwind die GréBe 


= we — > Vue x [— cos w + sin pp — i (cos w + sin y)] exp (ivt—ipé) 8) | 


ch) 
mit der Amplitude 
Lae a2 a we. (19) if 


U9 


b) Ergebnisse der ebenen Theorie des ruhenden oder gleichférmig bewegten | | 
Tragfliigels in periodisch veranderlicher Strémung. Die hydrodynamischen Kraft- | 


wirkungen, die an einem ruhenden oder sich gleichférmig bewegenden Tragfliigel groBer Spann- 
weite bei periodischer Anderung des hydrodynamischen Anstellwinkels auftreten, kénnen mit | 
Hilfe einer vom Verfasser entwickelten und auf der bereits erwahnten Arbeit* aufbauenden 
Theorie! berechnet werden. Die Ergebnisse dieser theoretischen Untersuchungen, die fiir die Aus- 
wertung und die physikalische Deutung der im Abschnitt 3 mitgeteilten MeBergebnisse wichtig 
sind, werden hier in Kiirze mitgeteilt. 


Der Anstellwinkel des als nicht gewélbt angenommenen Tragfliigels ist gegeben durch die | 
Beziehung 


a(x, t)=a)+—~ 3 |w,| exp (inyt—inps), (20) 
2) n=1 


wobei der Summenausdruck der durch die periodischen Anderungen der Strémungsrichtung her- 
vorgerufene zusatzliche hydrodynamisch bedingte Anstellwinkel ist. Da die im Abschnitt 2 be- 
schriebenen Messungen an einem Tragfliigel mit dem mittleren Anstellwinkel Null in einem in- |} 
stationdaren Strémungsfeld mit nahezu rein sinusférmigem Verlauf der normal zum Tragfligel |} 
gerichteten Stérgeschwindigkeiten durchgefiihrt wurden, vereinfacht sich der Ausdruck fiir den |} 
Anstellwinkel zu 

|r| 


a(x, t) = — exp (ivt—ipx). (20a) i} 


Vo 


Mit einem Abklingungsfaktor o, der in der zitierten Arbeit naher beschrieben wird, ergibt sich 
fir emen hydrodynamischen Anstellwinkel entsprechend (20a) die Normalkraft des Tragfliigels zu 


Py = qb Um (Ja(u) UL + TH) —o)] + Ju 


OM Get | 


+3 [Ja(u) TCH) 0) + J(u) 


iM (Ou cea eee pe °|| exp(ivt) (21) 


bzw. 
ewe a“ lep,| exp [? Pp(u) + ig] (21a) | 


* Siehe FuBnote 3 von Seite 364. 
1 H. Schwanecke, Schiffstechnik 5 (1958) S. 97. 
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mit 
Jol? (1 + PH) Gol + TW) CL — oy} + JH) [rw a—o 
_ 3—T ; 
sae of |") (o—1)+ ae o|| 
+2 Solu) Jntu) {(1 + Tq) Uo) | 7”) (10) 37] 
cs A =e (22) 
+ T’(u) (l—o)|T(u) (0-1) +14 syed 
und 


(Jou) 7”) —o) + Hw] Toye—W) 41454 | 
Op (u) = arc tg 


Te ss 
Fw) 1 + TH) A —o)) + Jue) r" Gyo cra | 
Das Moment der Normalkraft um die Tragfliigeleintrittskante hat die GréBe 


M,,, =—qbP ms i ( o(/4) [1 + T’(u) (1 —o)] + J,(n) ED) (l 0) 2 “| 


+ Lut (la) —2 TW) (eI) +14(5P + 2) ol} exp (iv) (23) 
oder, analog (21a), 


+ Ji(1) 


My, =—4b 1 = ley, | oP [tomy (u) + ive], (23a) 
Die GréBe der Tangentialkraft ergibt sich aus der Beziehung 


P= Paw —9b1(2) 5 (ou) + TH) a + h(a) [TG @—a) —* GT oll 
| Folu) Tu) A —0) + Au) | Tw) (1) + 1+ GPa} 
+42 {Iolu) +7) Ao] + A(w) |T"w) 1-0) — Gah x 


x | Jou) T’’(u) (L—o) + J(u) rw @—1) +1458 o|}) exp (i291) (24) 
mit 


Pew =—9b1 (ES (13) [TP —20 + 04) + Tw) @—20) +] 
ee ery teat +4 in| + T(u) (2 Gao galt T"(u) = + ian 


+ Jo(u) Ja(u) [ITUNES 6 + TW) T+ Tw) (440) 7), (25) 


Werden auch hier Beiwerte eingefiihrt, so ae aus (24) 


P,= qb (RY (cee \er,| exp (12 Dp (u) +i2yt)] (24a) 

"4 —— (Pem nach (25)) 
cee gl eal 

und oe 

( (ney 0 + ray a—on + 409] 7°09 — 1 dx 

o nar t 2 TR a ee I T’ 2 ese 
oe ors [Jo [1+ TW ol + J) rw —a)—— me) a va 
iB 5) dl 


x [Jo T"(u) (1 —0) + x(n) rw (— 


wt “a (26) 
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Der instationare Giitegrad, der ein MaB fiir den Abbau der normal zum Tra 
Geschwindigkeitskomponenten des Strémungsfeldes ist, wurde in Gleichung (61) 
i Arbeit* in der Form 


definiert. L, und L, sind die in dem instationaren Strémungsfeld enthaltenen sekundlichen Energie- 
betrige vor und hinter dem Tragfliigel. Sie sind dem Quadrat der Amplituden der normal zum 
Tragfliigel gerichteten Geschwindigkeitskomponenten proportional. Werden die ungestorte Nor- 
malgeschwindigkeitsamplitude mit |w| und die durch die Tragfliigel beeinflu8te Amplitude mit 
\w,| bezeichnet, so kann der instationare Giitegrad in der Form 


Neacas! aa i (27) 


|2e0| 


angegeben werden. 
Sinngema8B wie im stationaren Fall kann auch bei instationarer Strémung eine Gleitzahl 


Pm Pam 
Einst = IP, -— || (28) 


Z |ep,| 


definiert werden. 
Teg re Die instationire Gleitzahl ist dann von Bedeutung, wenn man den Tragfliigel zur Vortriebs 
erzeugung benutzen will. Sie gibt an, wie groB die mit einem geforderten Vortrieb zwangslaufig 
auftretenden Normalkraftamplituden sind. 


c) Bezeichnungen. 


I Profillange, 

d Profildicke, 

b Spannweite (Fligelbreite), ; 

Bo Schlagamplitude des Schlagfliigels, 

t Zeit, 

Y Kreisfrequenz, 

f Schlagfrequenz des Schlagfliigels, 

Vo mittlere Anstrémgeschwindigkeit, konstant, 
w(t) Zusatzgeschwindigkeit senkrecht zur Richtung der mittleren Anstrémgeschwindigkeit, 

|u| - zugehérige Amplitude, 


i de desgl. dimensionslos, 
: 0 
— q Staudruck, 
% Rechtwinkliges Koordinatensystem, die x-Achse zeigt in die Richtung der mittleren 
iy Anstrémgeschwindigkeit, die z-Achse senkrecht nach oben. Die Bezugslange fiir die — 
z dimensionslosen Koordinaten ist die halbe Profillange, -) 
By Integrationsveranderliche am Ort des Profiles, a 
é _ Aufpunkt im Strémungsfeld auf der x-Achse, : 
= 1 = 
= es reduzierte Frequenz, ; | 
T(u) Komplexe Funktion**, : 


J(u), J,(u) Besselfunktionen, 
a C,(u), S,(u) Integralcosinus bzw. -sinus, : 
—o(u) __ Abklingungsfaktor fiir die Geschwindigkeitskomponenten senkrecht zur Tragfliigelober- 


eres flache, Z 
Cp,Cy _—— Krraft- bzw. Momentenbeiwert, <a 
a ‘(|4 | Amplitude einer komplexen GréBe Z = Z’ + iZ’. 


T 


(Eingegangen am 28. September 1960.) 


Anschrift des Verfassers: Dr.-Ing. Helmut Schwanecke, Berlin NW 21, Schleuseninsel im Tiexparcent 


(ethane aig tetas Fate —_ ways nx iees 


* Vgl. FuBnote 1 von Seite 366. ~~ 
** Vel. FuBnoten 3, 4 von Seite 364 und FuBnote 1 von Seite 366. 
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Spannungsoptik 
Kin Lehr- und Nachschlagebuch fiir Forschung, Technik und Unterricht 


Von Dr.-Ing. Helmut Wolf, 


Leiter des spannungsoptischen Laboratoriums der Siemens-Schuckertwerke AG, Miulheim/Rubr 


Mit 311 zum Teil farbigen Abbildungen im Text und auf einer Tafel, 3 Tabellen und 1119 Literatur- 


hinweisen von 1814 bis 1959. XX, 582 Seiten Gr.-8°. 1961. Ganzleinen DM 66,— 


ZUR INFORMATION 


Im vorliegenden Buch werden, aufer den elastizitatstheoretischen und optischen Grundlagen der Spannungs- 
optik, die wichtigsten spannungsoptischen Verfahren zur Untersuchung ebener und raumlicher Spannungs- 
Rasta ausfiihrlich behandelt. Verschiedene Kapitel sind den besonders fiir die Praxis wichtigen Fragen 
des Modellmaterials, der Modellherstellung und der Versuchstechnik gewidmet. Auch spezielle Anwen- 
dungen des spannungsoptischen Verfahrens, wie z.B. die Untersuchung von Rotationsproblemen und 
elastoplastischen Spannungszustanden, werden beschrieben. In einem groBeren Literaturverzeichnis sind 
alle wichtigen spannungsoptischen Veréffentlichungen des In- und Auslandes, beginnend mit den ersten 
Anfangen der Spannungsoptik, zusammengestellt. Das Buch wendet sich, auBer an Studierende des Ma- 
schinenbaus, des Bauingenieurwesens, der Physik und Mathematik, an alle diejenigen Ingenieure, Kon- 
strukteure, Physiker und Mathematiker, die sich in Lehre, Forschung oder Praxis mit Spannungs- bzw. 


Festigkeitsproblemen beschaftigen. 


INHALTSUBERSICHT 


Die elastizitatstheoretischen und optischen Grundlagen - Die spannungsoptische Grundausriistung - Das 
spannungsoptische Grundgesetz in zwei Dimensionen - Isochromaten - Isoklinen und Hauptspannungslinien - 
Die wichtigsten Verfahren zur Bestimmung ebener Spannungszustande - Ubertragung der Modellversuche - 
Modellmaterial - Zur Modeltherstellung - Versuchstechnik und photographische Ausriistung - Die wichtig- 
sten Verfahren zur Bestimmung raumlicher Spannungszustande: Zur Technik dreidimensionaler spannungs- 
optischer Untersuchungen - Spezielle Anwendungen spannungsoptischer Methoden - Literaturhinweise fiir 


die einzelnen Abschnitte - Literatur-, Namen- und Sachverzeichnis. 
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Einflu8felder orthogonal anisotroper Platten 


Influence surfaces of orthogonal anisotropic plates 


Von Dr.-Ing. Siegfried Krug und Dr.-Ing. Peter Stein, 


Technische Universitat Berlin 


In deutscher und englischer Sprache 


Englische Ubersetzung durch Dipl.-Ing. HEINZ JUHL, Harrisburg, USA 


Mit 31 Abbildungen im Text und 193 Tafeln, VIII, 32 Seiten. Quer 4°. 1961. Ganzlemen DM 67,50 


ZUR INFORMATION 


Dieses Buch wendet sich an den praktisch tatigen 
Ingenieur und erleichtert ihm die Berechnung ortho- 
gonal anisotroper Rechteckplatten des Leicht-, 
Stahl- und Stahlbetonbaues. 


Die MomenteneinfluBfelder sind nach den Rand- 
bedingungen des freien, des frei drehbaren und des 
eingespannten Randes sowie deren Kombinationen 
geordnet. Die verschiedensten Seitenverhaltnisse, 
die unterschiedlichen Biegesteifigkeiten in der 
Langs- und Querrichtung sowie die Torsionssteifig- 
keit werden durch zwei Parameter ¢ und y erfaBt. 
Es wird gezeigt, daS durch die ausgewablten Werte 
€ = 2,031,253; 1,0; 0,8; 6,5 und y = 0,8; 0,4; 0 prak- 
tisch alle Seiten- und Steifigkeitsverhaltmisse er- 
fat werden kénnen, da fiir zwischenliegende Werte 
einfache Interpolationsformeln angegeben sind. 


Auf die Beriicksichtigung einer beliebigen Quer- 
kontraktionszahl sowie auf die Verwendung der 
EinfluBfelder zur Berechnung durchlaufender sym- 
metrischer Platten ist hingewiesen. 


Fir gleichmaBig verteilte Vollbelastung sind die 
Biegemomente tabellarisch angegeben. 


Die Anwendung der Tafeln, insbesonderé die Aus- 
wertung sowie die Interpolation, werden durch Bei- 
spiele in einfacher Weise. dargeleet. 


INHALTSVERZEICHNIS 


Einleitung. Transformationen zur Berechnung einer 
heliebigen orthotropen Platte. Formeln zur Be- 
rechnung der Momente. Auswertung der Einfluf- 
felder, Interpolation und Extrapolation. beliebiger 
Platten. Allgemeine Betrachtungen iiber die An- 
wendung dér HinfluBfelder. Literaturverzeichnis. 
Tafelverzeichnis. Momente aus gleichmafig ver- 
teilter Vollbelastung. EinfluBfelder Nr.1 bis 
Nr, 193. 


FOR INFORMATION 


This book appeals to the practising engineer. It 
simplifies the analysis of rectangular plates, ortho- 
gonally anisotropic, of the light construction and of 
steel and reinforced concrete structures. 


The moment influence surfaces are arranged ac- 
cording to edge conditions, as free, hinged and fixed 
edges and to their combinations. The different ben- 
ding stiffnesses in longitudinal and transverse 
direction, and the torsional stiffness are being con- 
sidered by two parameters, ¢ and y. Itis shown that 
practically all side length and stiffness ratios can 
be considered with the chosen values of ¢ = 2,0; 
1,25; 1,0; 0,8; 0,5 and 7 = 0,8; 0,4; 0, because 
simple formulas of interpolation are given for values 
lying between. 

The general consideration of any Poisson’s ratio, 
and the application of the influence surfaces for the 
analysis of continuous, symmetrical plates is 
mentioned. 

The evaluation for the bending moments for uni- 
formly distributed load has been done. 


The application of the charts, particularly the 
evaluation and the interpolation, are explained in 
an easily understandable way by examples. 


CONTENTS 


Introduction. Transformations for the analysis of 
the general orthotropic plate. Formulae for the 
computation of the moments. Evaluation of the 
influence surfaces. Interpolation and extrapolation 
of any general plate. General considerations on the 
application of the influence surfaces. Index of re- 
ferences. Chart index, Moments for uniformly 
distributed load over the entire plate. Influence 
surfaces No. 1 to No. 193, 
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